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Alternative Theories of Gravity 

Timothy Clifton 
SUMMARY 

This work investigates alternative theories of gravity, the solutions to their field 
equations and the constraints that can be imposed upon them from observation and 
experiment. Specifically, we consider the cosmologies and spherically symmetric so- 
lutions that can be expected to result from scalar-tensor and fourth-order theories 
of gravity. We find exact cosmological solutions of various different kinds; isotropic 
and anisotropic, homogeneous and inhomogeneous. These solutions are used to in- 
vestigate the behaviour of the Universes at both late and early times, to investigate 
the effects of corrections to general relativity on approach to an initial singularity 
and to look for effects which may be observable in the present day Universe. We 
use physical processes, such as the primordial nucleosynthesis of the light elements, 
to impose constraints upon any deviations from the standard model. Furthermore, 
we investigate the vacuum spherically symmetric solutions of these theories. This 
environment is of particular interest for considerations of the local effects of gravity, 
where the most accurate experiments and observations of gravitational phenomena 
can be performed. Exact solutions are obtained for this situation and their stabihty 
analysed. It is found that a variety of new behaviour is obtainable in these theo- 
ries that was not previously possible in the standard model. This new behaviour 
allows us an extended framework in which to consider gravitational physics, and its 
cosmological consequences. 
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CHAPTER 1 



Introduction 



The standard model of gravitation and cosmology is based on the theory of general 
relativity (GR). GR is one of the cornerstones of modern theoretical physics and 
has been shown to be enormously successful not just in its conceptual ingenuity and 
mathematical elegance, but also in its ability to explain real physical phenomena. 
This is particularly astonishing considering that the vast majority of relativistic 
gravitational effects were completely unknown at the time of its conception. 

Whilst in most branches of physics it has been the case that experiments have 
been performed and theory developed to fit these observations, this has not been the 
case for the development of relativity theory. Instead, due to the extreme weakness 
of gravity, it has been the situation that the theory of relativity was developed with 
only a minimal input from experiment. The theory then had only to satisfy the few 
experiments that were available at the time. Today this is no longer the case; an ever 
increasing body of experimental and observational evidence is being accumulated. 
This wealth of data has brought gravitational physics up to the status of a real 
experimental science, and so it must be confirmed that the theory which has been 
developed does in fact adequately model the observations. 

In making any comparison between theory and observation it is useful to have 
a framework in which to consider the theory. Such frameworks have been con- 
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structed for gravitational physics by a number of authors for a variety of different 
circumstances. For example, the parameterised post-Newtonian approach has been 
developed to provide a framework in which weak field tests of gravity can be inter- 
preted. Once the parameters in this very general approach have been constrained 
by the appropriate observations, these constraints can then be compared to the pre- 
dictions of the relativistic theory. However, if we have only a single theory it is 
difficult to interpret exactly what these observations mean; we can only say if they 
are consistent with that one theory or not. In order to make a proper evaluation 
of just how influential these constraints are it is necessary to have a variety of the- 
ories, each differing from the others in their own particular modiflcations. It is the 
development of just such alternative theories that is the subject of this work. 

As well as providing a foil against which to test GR, alternative theories of grav- 
ity are also interesting to study in their own right. The very existence of logically 
consistent alternatives to GR demands investigation to establish their potential com- 
patibility with experiment, how they manifest the usual relativistic phenomena and 
if any new phenomena can be identifled that are not present in the standard model. 
This avenue of research becomes particularly pressing if one wishes to attempt to ex- 
plain aspects of gravitational physics which appear unsatisfactorily accounted for in 
the standard model. A number of these apparent shortcomings present themselves 
both from observations of physical phenomena as well from more philosophical con- 
cerns regarding the theory itself. Examples which are often cited in the literature 
are the apparent need for more than 95% of our Universe to be made from dark 
matter and dark energy (neither of which is satisfactorily understood), the anoma- 
lous acceleration of Pioneer, the incompatibility of GR with quantum fleld theory 
and the lack of application of Mach's principle. Any one of these provides sufficient 
motivation to investigate potential alternatives to GR. 

In considering alternative theories of gravity it is necessary to impose some a 
priori conditions, to limit the number of candidate theories. We proceed as Dicke 
[HT], and choose that 

1. Space-time is a differentiable four-manifold with a metric and a connection. 

2. The theory must be relativistic. 
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3. General coordinate covariance should be manifest. 

4. There should exist no a priori geometric structure. 

5. The field equations should be derivable from an action principle. 

6. The theory should be simple. 

The first condition is the statement that we are restricting ourselves to what are 
often called 'metric theories of gravity'. Associating a metric with the curvature 
of space-time is a powerful idea, leading directly to the equivalence of inertial and 
gravitational mass as neutral test-particles (those with negligible self-gravity) follow 
geodesies of the space-time. This equivalence of gravitational and inertial masses has 
been shown to very high accuracy by Eotvos type experiments. Furthermore, this 
condition limits our considerations to four dimensional theories with Riemannian 
geometry. 

The second condition can be well justified from both experimental and philosoph- 
ical considerations and is closely related to the first. It is known that at any point 
on a smooth manifold we can make a local coordinate transformation to normal co- 
ordinates. In such a coordinate system the effects of gravity should be transformed 
away (up to tidal forces) and we should recover special relativity. This is a mani- 
festation of the weak equivalence principle which states: ' a freely falling frame in a 
gravitational fi,eld is equivalent to an inertial frame in the absence of gravity\ 

The condition of general coordinate covariance implies that we should be free to 
choose any coordinates to mark the positions of physical events on our space-time 
manifold. It is the position of the events themselves and the geometry of space-time 
that is important, not the arbitrary set of coordinates we put on the manifold. This 
condition requires that the theory should be formulated in terms of tensors. 

Condition four is the idea that there should be no absolute space. It is the 
requirement that the geometry of space-time should be entirely determined by its 
matter content. This condition is in keeping with Mach's principle which states that 
the inertial forces acting on a body can only be defined with respect to other bodies, 
and not with respect to abstract concepts such as absolute space-time. 

The condition that the field equations of the theory should be derivable from an 
action principle gives us covariantly conserved quantities, such as energy-momentum. 
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via Nother's theorem and allows us to use Lagrangian mechanics. Further justifi- 
cation for this condition may be found in quantum mechanics; if a satisfactory 
quantum theory of gravity is to be found this strongly suggests that the theory 
should be derivable from an action principle. 

The last of the above conditions is perhaps the most tenuous and is the one which 
could be most readily relaxed. We choose to interpret this condition as meaning 
that the gravitational Lagrangian should be a function of simple contractions of the 
Riemann tensor and one extra scalar field only. This still leaves us with a great 
number of alternative theories, including the most general formulation of the scalar- 
tensor theories as well as the fourth-order theories of gravity. Both of these classes 
of theory contain GR as a special (or limiting) case and it is these theories which 
this work will be concerned with. 

It is clear that GR satisfies all of the above conditions. However, there are a 
number of additional conditions that GR also obeys which will not necessarily be 
obeyed by the theories we will be considering. These conditions are: 

1. The strong equivalence principle. 

2. The space and time independence of Newton's constant. 

3. The linearity of the action in second derivatives of the metric. 

The first of these conditions, the strong equivalence principle, is the statement 
that: 'massive self- gravitating objects should follow geodesies of the space-time'' . 
Whilst the metric formahsm above assures us of this property for neutral test- 
particles, we are no longer assured that it should be true in general for extended 
massive bodies. These objects distort the background space-time in which they exist 
and whilst it is true that in GR they still follow geodesies, this is no longer true in 
general for the alternative theories we will be considering. This type of violation of 
the strong equivalence principle is often referred to as the Nordvedt effect. 

The second condition that is met by GR is that at any point in space-time one 
can take the weak-field and slow motion limit to recover Newtonian gravity, with the 
same value of Newton's constant as is measured on Earth. This condition is explicitly 
violated by some of the theories we will consider. In fact, the Brans-Dicke theory 
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was devised specifically in order to create a theory of gravitation in which space- 
time variations of Newton's constant could occur. Constraining the parameters of 
these theories is then equivalent to constraining the potential variations of Newton's 
constant which they model. 

The violation of the last of these three conditions is perhaps the most worrying. 
It is the hnearity in second derivatives of the Einstein-Hilbert action that ensures 
that the field equations of GR are of no higher than second order in derivatives of the 
metric. Whilst this is also true of scalar-tensor theories, it is no longer true of theories 
which are constructed from contractions of the Riemann tensor only (except for 
the Einstein-Hilbert action with a cosmological constant). Generating Lagrangians 
which are non-linear in second derivatives yield field equations which are fourth-order 
in derivatives of the metric (with the exception of additional terms in the Lagrangian 
which are total diverges, such as the Gauss-Bonnet term). Also, variation of the 
action with respect to the metric is then no longer equivalent to Palatini approach 
of varying with respect to the metric and an independent connection. In fourth- 
order theories these two approaches yield two different sets of field equations. In 
this work we will only be concerned with the metric variation approach. 

The investigation of alternative theories of gravity is hampered by the additional 
complexities that modifications to GR introduce. It seems that GR is unique not 
only in satisfying all of the conditions listed above, but also in being the simplest 
relativistic metric theory of gravitation that can be conceived of. Any modification 
to GR invariably results in a set of field equations that is considerably more com- 
plicated than Einstein's equations. In dealing with this complexity there appear in 
the literature two difi^erent ways of proceeding. The first is to look for approximate 
solutions to a specific theory of interest. This approach is particularly useful in 
the weak-field and low-velocity regime where frameworks such as the parameterised 
post-Newtonian approach can often be employed. It is the approach that is most 
frequently taken up when investigating theories that are motivated by a desire to 
overcome the perceived shortcomings of the standard theory (examples of which 
were given above). The second approach is to look for exact solutions of the sim- 
plest possible modifications to GR. The idea behind this approach is to understand 
as well as possible the effect of modifying the standard theory. Once this behaviour 
is well understood it can then be used as an approximation to more complicated 
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theories of specific interest, as well for considerations of that particular modification. 
It is the second approach that will be most closely adhered to in this work. 

The program for this study will be to investigate the form of the solutions to 
the simplest alternative theories that meet the demands listed above; these will be 
the scalar-tensor theories and the fourth-order theories, which are often collectively 
referred to as the extended theories of gravity. Particular attention will be focussed 
on highly symmetric situations. High symmetry space-times are the most readily 
solved for and are often the ones of most physical interest. These solutions will then 
be used to model physical processes that occur in the Universe. Comparing these 
models with observation will allow constraints to be imposed upon the theory, which 
will limit the allowed deviations of this kind from GR. In performing this analysis 
we will be able to note new behaviour that was not previously obtainable in the 
general relativistic case. 

We begin in chapter 2 by formally introducing the theories that will be consid- 
ered. The generating Lagrangians will be explicitly given and the field equations 
derived from them. The form of these field equations will immediately be seen to 
be more compficated than the general relativistic case. We will then give a brief 
discussion of how these theories can be conformally transformed to look like general 
relativity minimally coupled to a scalar field. This conformal equivalence will be 
useful in finding solutions to the field equations. 

In chapter 3 the homogeneous and isotropic solutions to the scalar-tensor field 
equations will be given. For the simplest case of a spatially fiat Brans-Dicke universe 
containing a single perfect fiuid the general solutions to the problem can be found. 
These general solutions have a late time power law attractor for the evolution of 
the scale factor and at early times are generally dominated for a period by the free 
component of the gravitational scalar itself. For this reason the general solution 
for the radiation dominated universe is presented explicitly, as it is this component 
of a realistic universe which is expected to dominate at early times, and the late 
time attractors arc given for dust and vacuum energy densities. Vacuum solutions 
with non-zero spatial curvature are also given, and will be used in chapter 7 for a 
discussion of inhomogeneous cosmologies. The spatially-fiat perfect fiuid solutions 
will be used in chapter 9 to model the expansion of the Universe during primordial 
nucleosynthesis. The homogeneous and isotropic cosmological solutions for more 
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general scalar-tensor theories are also given in this chapter. 

In chapter 4 the homogeneous and isotropic spatially flat cosmologies of a simple 
class of fourth-order gravity theories are investigated. A phase plane analysis is 
performed of the solutions and it is found that the late time attr actor solutions for 
both vacuum cosmologies and universes filled with perfect fluids are simple power- 
law functions of time. These solutions are used in chapter 9 to model the expansion 
of the Universe during both primordial nucleosynthesis and microwave background 
formation. 

In chapter 5 a variety of other homogeneous fourth-order cosmologies are pre- 
sented. These include solutions of the Godel type, as well as Einstein static, de-Sitter 
and Bianchi type I Kasner-like solutions. For the case of the Godel solutions the 
existence of closed time-like curves, and hence the theoretical possibility of time 
travel, is investigated. The Kasner-like solutions are used to perform an analysis of 
the approach towards an anisotropic cosmological singularity in these theories. It is 
found that, unlike in the general relativistic case, these cosmologies do not undergo 
an inflnite sequence of chaotic oscillations as the initial singularity is approached. 
It is also conflrmed that matter sources do not generally dominate the evolution of 
these universes on approach to the singularity. 

In chapter 6 static and spherically symmetric vacuum solutions are presented for 
both the scalar-tensor and fourth-order cases. For the case of the simple class of 
fourth-order theories that is studied it is found, using a dynamical systems approach, 
that the general solution to the problem is not asymptotically flat. An explicit r 
dependent expression is found for the asymptotic attractor and an exact solution 
is presented which has the same asymptotic behaviour. A perturbative analysis is 
performed which conflrms that the exact solution corresponds to the mode with the 
appropriate Newtonian limit and also indicates the existence of extra modes which 
are damped oscillatory waves. This solution is used in chapter 10 to model solar 
system gravitational experiments in such a space-time. 

Chapter 7 investigates the inhomogeneous cosmologies of these alternative the- 
ories. In the scalar-tensor case exact solutions are presented for inhomogeneous 
universes which asymptotically approach the homogeneous and isotropic solutions 
of chapter 3. A detailed investigation of the spherical collapse model in these the- 
ories is also performed: firstly for exact fiat and closed vacuum solutions which 
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are matched at a boundary and secondly for flat and closed universes containing a 
perfect fluid. For the case of spherical collapse with a perfect fluid numerical simula- 
tions are performed, as no exact solutions are know for closed perfect fluid universes. 
The evolution of the scalar field, and hence Newton's constant, are followed both 
inside and outside the collapsing region and variations are tracked. For the case of 
fourth-order theories an exact spherically symmetric vacuum solution is presented 
which asymptotically approaches the homogeneous and isotropic vacuum solution 
of chapter 4. This solution is used in chapter 10, as an alternative to the static 
solution found in chapter 6, to model gravitational experiments that are carried out 
in the solar system. A linear perturbative analysis is carried out on the background 
of this solution and it is found that the exact solution does indeed reduce to the ap- 
propriate Newtonian limit and that the extra modes present in the general solution 
are manifest as modes growing as in the linearisation. 

In chapter 8 an analysis is made of scalar-tensor cosmologies where energy is 
allowed to be exchanged between the gravitational scalar field and a perfect fiuid. 
In such a theory all of the conditions required of the theory above are still met. 
However, some of the meaning of condition 1 is lost as test-particles no longer fol- 
low exact geodesies of the space-time geometry. This is due to the interaction with 
the gravitational scalar, which acts as an external force on the test-particle. The 
exchange of energy and entropy between the scalar field and the perfect fiuid is 
discussed and the evolution of homogeneous and isotropic universes is then consid- 
ered. For a general form of the parameter controlling the energy exchange the field 
equations for the scale factor and scalar field can be decoupled and integrated to 
two first order ordinary differential equations. From this point some simple spe- 
ciflc functions of the energy exchange parameter are proposed and these equations 
are then integrated again to give explicit expressions for the evolutions of the scale 
factor and scalar fleld. 

In chapter 9 cosmological processes are modelled in the solutions that were pre- 
sented in previous chapters. Speciflcally, primordial nucleosynthesis is modelled in 
the general solutions for radiation dominated spatially flat scalar-tensor universes. 
Previous analyses have been performed on primordial nucleosynthesis in these the- 
ories, but have always taken the late time power law solutions to describe the evo- 
lution of the Universe. In this study we consider the general form of the solution 
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and hence allow for Newton's constant to be varying whilst nucleosynthesis occurs. 
This allows, for the first time, observational constraints to be imposed upon the 
free scalar-dominated period of such a universe's expansion. In addition, we study 
primordial nucleosynthesis and microwave background formation in the power law 
solutions of the fourth-order theories found in chapter 4. These observable processes 
allow the theory to be constrained from observation. 

Chapter 10 provides an analysis of the weak field constraints available on the 
theory. A brief explanation of the parameterised post-Newtonian approach, and 
how it can be used to place constraints on the scalar-tensor theories of gravity, is 
given. The spherically symmetric vacuum solutions to the fourth order theories 
found in chapters 6 and 7 are then analysed. It is shown that these theories do 
not fit so readily into the usual formalism, and that in these theories one should 
expect extra gravitational forces in the weak field limit which are not present in 
GR. For the static solution this extra force takes the form of a force which drops off 
as and for the inhomogeneous cosmological solution this force looks hke a friction 
force, for a slowly moving test-particle. The r"^ force in the static solution has a 
dramatic impact on the perihelion precession of closed orbits in such a space-time 
and, applying observational constraints from the orbit of Mercury, it is possible to 
constrain the theory dramatically if this solution is used to model the solar system. 
The inhomogeneous cosmological solution fits better into the parameterised post- 
Newtonian formalism and, if it is used to model the solar system, can be used to 
place tight constraints on this theory from the observations of the Shapiro time delay 
of radio signals from the Cassini space craft as it passes behind the Sun. The extra 
gravitational forces present in these solutions, and their lack of asymptotic flatness, 
show explicitly new behaviour in these theories which is not present in GR. 

In chapter 11 we summarise our results and provide concluding remarks. 
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Conventions 

The Riemann curvature tensor is 

^ bed — ^ bd,c — i 6c,(i + i ce^ bd " i de^ c& 

wliere Latin indices refer to space-time components and we use Greek indices to 
refer to space components. The signature of the metric is taken to be (— , +, +, +) 
and we choose units in which c = 1. 



10 



CHAPTER 2 



Extended Theories of Gravity 



'Extended theories of gravity' is the collective term usually apphed to the scalar- 
tensor and fourth-order theories of gravity. In this chapter we will introduce these 
theories, deriving their field equations and showing their conformal equivalence to 
GR. 

2.1 Actions and Field Equations 

We present in this section the generating Lagrangians for the theories that will be 
considered in the remainder of this work. The actions obtained from integrating 
these Lagrangians over all space will be extremised to obtain the relevant field 
equations. 

Considering gravitational theories derived from Lagrangians has the advantage 
of ensuring coordinate independence of the derived field equations, as well as the 
covariant conservation of quantities via Nother's constraint. 
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2.1.1 Scalar-Tensor Theories 

A general form of the scalar-tensor theory 1143t I187j can be derived from the 
Langrangian density 

>C = T^v^ ifi<P)R - 9i^)da^dy - 2A(0)) + C^"^, h{(P)gab) (2.1.1) 

where f,g,h and A are arbitrary functions of the scalar field and Cm is the La- 
grangian density of the matter fields \1'. The function /;.(0) can be absorbed into the 
metric by a conformal transformation of the form jHO] 

h{(p)gab gab- 

The conformal frame picked out by this choice is the one in which there is no direct 
interaction between the scalar field and matter fields and is usually referred to as 
the Jordan frame (see section 12.21 for a discussion of conformal transformations in 
these theories). Test-particles in this conformal frame follow geodesies of the metric 
and the weak equivalence principle is obeyed. The effect of this transformation on 
the rest of the Lagrangian can then be absorbed into redefinitions of /, g and A. 

By a redefinition of the scalar field we can now set /(0) 0, without loss of 
generality. The Lagrangian density ()2.1.H) can then be written as 

^ = T^v^ (<PR - ^9,0a'^0 - 2A(0)^ + £„(^, gab) (2.1.2) 



IGvr 

where uj{(f)) is an arbitrary function which has been chosen in analogy with the Brans- 
Dicke theory and A is a 0-dependent generalisation of the cosmological constant. 
This theory reduces to the Brans-Dicke theory in the limit u —inconstant and A ^ 0, 
and approaches general relativity in the limit u ^ oo, u'/uj"^ — and A — 0. 
The variation of the action derived from integrating ()2.1.2|) over all space is 



51 = -— 

lOTT 



^g^'dgabC + <P5R + ^d^<pd'<p5gab + ^'^T'^'^gab 



(2.1.3) 



where use has been made of 



Sg""' = -g^^g'^Sgcd, 5^g = \g^'5gab (2.1.4) 
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and 

rpab _ ^ 

Now using the relation [^Oj 

^ -cdt 



we can write 



where 



^Rab = -g'' {SgacM + Sgbd;ac " Sgab;cd " '^fl'cdjab) (2.1.5) 

- 1^ / ^^v^^-*^"' + 8^T-')^9ab (2.1.6) 



gab ^ ^^j^ab _ l^a5^) _ ^^aCgbd _ gabgCd^^^^^ _ ^^gac^bd _ Ig'^'^g^^)^^^^^^ + ^-^A 

and ~ means equal up to total divergences. Such terms are irrelevant here as 
they can be transformed via Gauss's theorem to terms on the boundary which are 
assumed to vanish. 

The first variation of the action corresponding to ()2.1.2|) . with respect to the 
scalar field 0, is now 

51- f dQ^/^ (r + 2uj— - ^0 + —0 c0' - 2A' I S(f) (2.1.7) 



where primes denote differentiation with respect to 0. Eliminating R with the trace 
of ()2.1.6|) and looking for stationary points of the action, by setting the first variation 
to zero, then gives the field equations 

cf^iRab-^gabR) = 87rT,6+(5^5,^-(7a6/')0;cd+^(5:5,'-^(?„6/')0,c0,d-^7a6A (2.1.8) 

and 

{2u + 3)00 = 8nT - cu'0,a0° - 4A + 20A'. (2.1.9) 
2.1.2 Fourth- Order Theories 

The fourth-order theories we shall consider are derived from functions of the three 
possible linear and quadratic contractions of the Riemann curvature tensor: R, 
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RabR"^^ and RabcdR"'^'^'^ ^ ■ The relevant weight-zero scalar density for this general 
class of theories is then given by 



C = x"'^f{X,Y,Z) 



(2.1.10) 



where /(X, F, Z) is an arbitrary function of X, Y and Z which are defined by 
X = R, Y = RabR"^^ and Z = Rabcd.R"''^'^'^ , and x is an arbitrary constant which 
can be determined from the appropriate Newtonian limit. The action is obtained, 
as usual, by integrating this density, together with that of the matter fields, over 
all space. The addition of supplementary terms to the density ()2.1.1()j) in order to 
cancel total divergences which can be transformed to integrals on the boundary can 
be problematic (see e.g. |121j ) and so, for simplicity, they will all be assumed to 
vanish. 

The variation of the action derived from integrating the density ()2.1.1()j) over all 
space is 



SI = 



X 



1 



fg^'Sgab + fx6X + fY5Y + fz6Z 



^fg^'Sgab - fxiR^'Sgab - g^'SRab) 



- 2fY{R'^''R'\5gab - R'^'SRab) 

— '^f z{Rcde^^ R'^^'^'^^'^gab — Ra'^^^R^'bcd) 

where use has again been made of ()2.1.4|) . Using the relations ()2.1.5|) and 



f2.1.111 



5R\cd - -g^^'iSged-M + 5g^b-4c - hdb.ec - 5g.c-M ' ^geb-M + 5g^b;ed). 



We can then write 



fxg^'^Rab ^ -fx-Ag^'g"' - g^'g'')5g 



ab 



2fYR''^5R, 



UifyR^') + {fyR 



cd\ 



cdg 



ab 



2{fyR 



Sg, 



ab 



2fzRa'"'5R%cd ^ ^fzR'^'''^'');cdSgab. 



^As pointed out by Dc Witt IW there is a fourth possibility, namely £'^'"^'^RefabR'^^cd- However, 
this contraction is of no physical interest due to parity considerations. 
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Substituting these expressions back into ()2.1.1H) then gives 



where 



P'^' = -^fg'^' + fxR"' + 2fYR'^'^R''\ + 2fzR"'''^''R'\de + fxAf'g"" - f'd"") 
+ □(/yi?'^') + g''\fYK%,d - 2{fyR<-)% - 4(/zi?'^('^^)^);,,. (2.1.12) 

The notation denotes partial differentiation of / with respect to A^. Looking for 
a stationary point of the action by setting the first variation to zero gives the field 
equations 

Pab=\Tab-gabk (2.1.13) 

when matter fields and a cosmological constant are included. Here, A is the cosmo- 
logical constant (defined independent of /(X, F, Z)) and is the energy- momentum 
tensor of the matter. These field equations are generally of fourth-order, with the 
exception of the cases in which the function / is linear in the second derivatives of 
the metric [7^ . 

We will often be considering gravitational theories derived from the choice 

f = f{R) = R'+' (2.1.14) 

where 5 7^ is a real number. The limit 5 — gives us the familiar Einstein- 
Hilbert Lagrangian of GR. We take the quantity to be the positive real root of 
R throughout this work. This choice of / gives 

Pa, = 5{1 - 6')R'^^ - 5{1 + 5)R'?f + (1 + 5)R'Rab - \gabRR' 

- gab6{l - 5^)R'?^ + 5(1 + S)gabR'^. (2.1.15) 

These equations have the useful property that in vacuum they are conformally equiv- 
alent to Einstein's equations with a scalar field in an exponential potential, as shown 
below. 
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The Newtonian Limit 



By comparing the geodesic equation to Newton's gravitational force law it can be 
seen that, as usual, 

C = (2.1.16) 

where $ is the Newtonian gravitational potential. All the other Christoffel symbols 
have r°j^ = 0, to the required order of accuracy. 

We now seek an approximation to the field equations ()2.1.15p that is of the form 
of Poisson's equation; this will allow us to fix the constant x- Constructing the 
components of the Riemann tensor from ()2.1.1(ij) we obtain the standard results 

^"'^' = ^ih ^oo = V^$. (2.1.17) 

The 00 component of the field equations ()2.1.15p can now be written 

{l + 5)R,,-]^g,,R=^^ (2.1.18) 

where terms containing derivatives of R have been discarded as they will contain 
third and fourth derivatives of $, which will have no counterparts in Poisson's 
equation. Subtracting the trace of equation ()2.1.18j) gives 

(1 + '^)^oo = ^ (too - (^-^-^^^ 

where T is the trace of the stress-energy tensor. Assuming a perfect-fluid form for 
T we should have, to first-order. 

Too ~ p and T ~ 3j9 - p ~ -p. (2.1.20) 

Substituting (|2.1.2()jl and (!2.1.17jl into (!2.1.19j) gives 

4(l-52)i?5- 
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Comparison of this expression with Poisson's equation allows one to read off 

X = 16^G^|r3§^o (2.1.21) 

where Rq is the value of the Ricci tensor at the time G is measured. It can be 
seen that the Newtonian limit of the field equations ()2.1.15|) does not reduce to the 
usual relation oc p, but instead contains an extra factor of . This can be 
interpreted as being the space-time dependence of Newton's constant, in this theory. 
Such a dependence should be expected as the Lagrangian ()2.1.14j) can be shown to be 
equivalent to a scalar-tensor theory, after an appropriate Legendre transformation^ 
(see e.g, |124j ). This type of Newtonian gravity theory admits a range of simple 
exact solutions in the case where the effective value of G is a power-law in time [20] . 



2.2 Conformal Transformations 

A conformal transformation of the metric gah into cjab can be written 

gab = e'^^^^gab (2.2.1) 



where r(a;) is an arbitrary function of the coordinates The line-element is then 
transformed according to 

7 9 Wt-r^ 



and the square root of the determinant of the metric as 



— 4r / — - 
-g = e V-^, 



in four dimensions. Conformal transformations of this kind are local, isotropic trans- 
formations of the standard of 'size', as defined by the line-element. These transfor- 



■^This equivalence to scalar-tensor theories should not be taken to imply that bounds on the 
Brans-Dicke parameter uj are immediately applicable to this theory. It can be shown that a 
potential for the scalar-field can have a non-trivial effect on the resulting phenomenology of the 
theory 147 . Furthermore, the form of the perturbation to general relativity that we are considering 
does not allow an expansion of the corresponding scalar field of the form (f)o+4>i where (f)o is constant 
and 1^1 1 << 1 00 1, so that any constraints obtained in a weak- field expansion of this sort cannot be 
applied to this situation. 
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mations are known as conformal as they leave the angle between any two vectors in 
the space-time unaltered. 

After performing such a transformation we use the term 'conformal frame' to 
distinguish the new, rescaled metric from the original. Among the infinite possible 
conformal frames there are two which are most commonly used and have specific 
interpretations: the Jordan frame and the Einstein frame. The Jordan frame is the 
one in which the energy-momentum tensor is covariantly conserved and in which 
test-particles follow geodesies of the space-time metric. The Brans-Dicke theory 
[52] . for instance, is most usually formulated in the Jordan frame. The Einstein 
frame is the conformal frame in which the field equations of the theory take the 
form of the Einstein equations (unlike the Jordan frame, the Einstein frame can 
only be defined for some theories). In the Einstein frame the field equations are 
second order but the energy-momentum tensor of the matter fields is not always 
covariantly conserved and test-particles do not necessarily follow geodesies of the 
space-time metric. The Einstein frame, therefore, is particularly useful for finding 
vacuum solutions, but less useful for finding solutions with matter fields present. 

Under the transformation 1)2.2.11) it can be shown that the Ricci tensor and scalar 
transform as (see e.g. |198j ) 

Rab = Rab - 2r.a, + 2r,„r,fe - 2^,b^^''r,,r,rf - gabOr (2.2.2) 
e^^R = R- 6^''^r,„r,6 - ear (2.2.3) 

and the d'Alembertian transforms as 

e^rn0 = Q0 + 2r'r>,fe (2.2.4) 

where overbars over operators or indices denote that they are defined using the 
metric cjab- We will now use these transformations to show how the scalar-tensor 
and some fourth-order theories can be transformed from the Jordan frames to the 
Einstein frame. 
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2.2.1 Scalar-Tensor Theories 

All of the scalar-tensor theories defined by the Lagrangian ()2.1.2|) . in the Jordan 
frame, can be transformed to an equivalent in the Einstein frame. To see this 
we will consider the various terms in the generating Lagrangian ()2.1.2j) seperately. 
Firstly, we consider the term containing the Ricci scalar 

lOTT 

which under the conformal transformation ()2.2.1|) becomes 

1 



= -i-V^<pe'^{R - Qr''^,ar,b - 6Qr). (2.2.5) 

The non-minimal coupling to the Ricci scalar can now be removed by making the 
choice of conformal factor 

e^^ = 0-^ such that gat = (2.2.6) 

This choice of F defines the conformal transformation between the Jordan and Ein- 
stein frames in the scalar-tensor theories. Making this choice of F it can be seen 
that the last term on the right-hand side of ()2.2.5|1 is a total divergence, which can 
be transformed to an integral over a surface at infinity and therefore ignored. 
Considering now the second term in ()2.1.2|) we get 

iOTT (p 47r 

under the transformation ()2. 2.111 and using the definition of F above. The third term 
then gives 

£3 = -J-^A = _±^e^^A 
so that the whole Lagrangian can be written 

^ = T^v^ - 2(3 + 2a;)^'^^F,,F,, - 2e^^A) + £^(^, e^^^,,). (2.2.7) 

lOTT ^ 
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Now, by making the definitions 



^'^ ^ and 8nV{^) = e^^A (2.2.8) 



(3 + 2uj) dip 

for tlie scalar ip and tlie function V{ip), we can write the transformed Lagrangian 
(ITT7|l as 

C = ^V^R + (-^rV,aV^,6 - Viij)^ + £^(vl>, e'^U)- (2.2.9) 

In the absense of any matter fields the scalar-tensor theories can now clearly be seen 
to be conformally related to Einstein's theory in the presence of a scalar field in a 
potential (that potential disappearing in the absense of A). 

In the Brans-Dicke theory [22] the coupling constant is a constant quantity 
and the scalar fields and ip are related by 



For more general theories with u = uj{(j)) the definition of ip, ()2.2.8p . must be inte- 
grated to obtain a relation between and ip. 

By extremising the action ()2.2.9|1 with respect to cjab we get the Einstein frame 
field equations 

Rab = 8n(Tab - -Qabf + 1p,a'^P,b) 

and by extremizing with respect to ip we get the Einstein frame propagation equation 

\I}ip = —y/AnaT, 

where a"^ = 3 + 2uj and we have defined the energy-momentum tensor Tab with 
respect to cjab so that T""^ = e^^T"'''. Correspondingly, whilst the Jordan frame 
energy-momentum tensor is always covariantly conserved, T^^.f, = 0, its counterpart 
in the Einstein frame is not, T!^'^ = \/A^aTip\ The interaction between the scalar 
field %p and the matter fields described by Tab can now be seen explicitly. 
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2.2.2 Fourth- Order Theories 



We will now show the conformal transformation between the Jordan frame and the 
Einstein frame for theories derivable from an abritrary function of the Ricci scalar 
|25| I122j . Extremising the action 

S = J dn^^f{R) (2.2.10) 

with respect to the metric gab gives the field equations, according to ()2.1.12j) . 

f'Rab - If gab - f[ab + QabUf = |T„, - (^^.A (2.2.11) 

where primes denote differentiation with respect to R. Under the conformal trans- 
formation 

Qab = f'gab, 

and making the definition 

X 

the field equations 1)2.2.111) then become 

Rab - \gabR = I (^4>,a4>,b - ^gabrU,c4>4 - gabV^ + (2.2.12) 

where 

(Rf-f) 2A 

2r Xf 

It can now be seen that, in the absense of any matter fields, theories derived from 
an action of the form ()2.2.10|) are conformally equivalent to GR in the presence of 
a scalar field in a potential. 

We will now consider the more specific case of theories derived from the La- 
grangian ()2.1.14|) . In this case the relevant conformal transformation can be written 
as (jab = ^oR^gab, whcrc Qq is a constant. The scalar field can then be defined as 



levrG 

and the field equations, in the absense of matter fields and any cosmological constant. 
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can be rewritten as 



Gab = SttG - ^^ab(rt,c0,<i + 21/(0))^ (2.2.13) 



and 

where is given by 



5sign(i?) 167rG(l-5) ^ . ^ . .^ 

The magnitude of the quantity Qq is not physically important and simply corre- 
sponds to the rescaling of the metric by a constant quantity, which can be absorbed 
by an appropriate rescaling of units. It is, however, important to ensure that fig > 
in order to maintain the +2 signature of the metric. 
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CHAPTER 3 



Homogeneous and Isotropic Scalar- Tensor 

Cosmologies 



We will consider in this section the homogeneous and isotropic cosmological solutions 
to the scalar-tensor fields equations ()2.1.8p and ()2.1.9p . These are space-times with 
a six- dimensional group of motions, Gq, acting transitively on a three-dimensional 
space-like subspace. 

3.1 Friedmann- Robertson- Walker Cosmologies 

Spatially homogeneous and isotropic space-times are described by the Friedmann- 
Robertson- Walker (FRW) line-element 



where a{t) is the scale factor and k is the spatial curvature. Using the FRW metric 
in the scalar-tensor field equations ()2.1.8|) and ()2.1.9|) . and setting A = 0, gives 



ds 



'2 = -dt"^ + a^(^t) 



dr"^ 



+ r'^{dd'^ + sm^ ddcj)'^) 



1 — kr'^ 




(3.1.1) 
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0^87r (p-3p) dj 

(2a; + 3) (2cu + 3)0' ^ ' ' 

and 

where H = a/a is the Hubble rate, an overdot denotes differentiation with respect 
to comoving proper time, t, p is the matter density, and p is the pressure. Each 
non-interacting ffuid source p{p) separately satisfies a conservation equation: 

p + 3H{p + p) = 0. (3.1.4) 

Substituting equations ()3.1.2p and ()3.1.3|) into equation ()3.1.ip gives 



H + H'^ - + = 8^ + 3p + pu) ^ 1 



UJ 



3 02 30 (2a; + 3) 2 (2a; + 3) 

It is the solutions of these equations that concern us in this chapter. 



(3.1.5) 



3.2 Brans-Dicke Theory 

For the Brans-Dicke theory (a; =constant) there are a variety of known exact ho- 
mogeneous and isotropic solutions (see e.g. Wl \ ITHH| I144j ). We will present the 
solutions here that are of particular significance, and that will be of use later on. 

In these theories a;(0) = a; is a constant. The three essential field equations for 
the evolution of the scalar field 0(t) and the expansion scale factor a{t) in a Brans- 
Dicke universe are ()3.1.2j) . ()3.1.3|1 . and ()3.1.4j) . Now a; is a constant parameter and 
the theory reduces to GR in the limit a; ^ oo where = G^^ — > constant. The gen- 
eral solutions to the FRW Brans-Dicke theory field equations are fully understood. 
The vacuum solutions are the t — attractors for the perfect-fluid solutions and 
the general solutions with equation of state 

P={l-l)p (3.2.1) 

and /c = can all be found. At early times they approach the vacuum solutions but 



24 



at late time they approach particular power-law exact solutions 




(3.2.2) 



0(t) = 0^i[2(4-37)/[4+3..7{2~7)] 



(3.2.3) 



These particular exact power-law solutions for a{t) and (/)(t) are 'Machian' in the 
sense that the cosmological evolution is driven by the matter content rather than 
by the kinetic energy of the free field. 

The power-law solutions ()3.2.2|1 and ()3.2.3|) are valid for any fluid with equation of 
state 1)3.2.11) . We will now consider seperately the solutions for universes dominated 
by radiation, matter and vacuum energy. The general solutions for A; = are known 
for any fluid that satisfies ()3.2.1|) |97j . Here we will only present the general solution 
for a radiation-dominated universe as the late-time power-law attractors for matter 
and vacuum dominated eras are expected to be sufficient for any realistic universe 
(see below). 

3.2.1 Radiation Domination 

We consider first the case of a fiat radiation-dominated universe. Assuming the 
equation of state p = \p and defining the conformal time variable, rj, by adrj = dt, 
the equations ()3.1.1|1 . ()3.1.2|1 and ()3.1.3|) integrate to give jHZI 



V + V2 



(3.2.4) 



if + (2?72 + 3?7i)?7 + rjl + 3r]ir]2 - \^'>il 
3?7i 



(3.2.5) 



0(r/) 



r]^ + (2?72 + 3?7i)?7 + r]l + 'ir]i'r]2 - ^ur]f 



where rji and ri2 are integration constants. 

For uj > —3/2 the solutions to these equations are 



a{r]) = ai{r] + r]+)' [r] + r]_ 




(3.2.6) 



+ 



(3.2.7) 
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where ?7± = ''72 + f??! ± 1^71 y ^ + f"^' '^1 ^^'^ '^i integration constants, and 
S7ipro/S(j)ia1 = 1 (subscript indicates a quantity measured at the present day 
and we rescale so that = 1, throughout). 
For uj < —3/2 we find 




a(r/) = aiJ{r] + r]_)'^ + r]lexp tan"^ ^ ^ , (3.2.8 




0(77) = 01 , tan-i , (3.2.9) 



where r]^ = f^iy ||^^| — V- = V2 + ^Vi^ Svrpro/S^iai = 1 and constants have been 
absorbed into ai and 0i. 

In these solutions rji determines the evolution of the scalar field during radiation 
domination, and is a physically interesting quantity; 172 sets the origin of the con- 
formal time coordinate. The evolution of a and (p, for uj > —3/2 and uj < —3/2, 
is shown in figures 13.11 and 13.21 respectively. A value of u = 40000 is chosen, in 
agreement with recent solar system observations 

For UJ > —3/2 {u! < —3/2) we see that the scale factor here undergoes an 
initial period of rapid (slow) expansion and at late times is attracted towards the 
solution a{'r]) oc 77, or a(t) oc . Similarly, can be seen to be changing rapidly at 
early times and slowly at late times. We attribute these two different behaviours, 
at early and late times, to periods of free scalar-field domination and radiation 
domination, respectively. In fact, setting r/i = in ()3.2.4j] and ()3.2.5|1 . we remove the 
scalar-dominated period and gain the power-law exact solutions 0(77) oc r] and = 
constant. These are the "Machian" solutions usually considered in the literature. 
Allowing rji to be nonzero we will have a nonconstant 0, and hence G, during 
radiation domination. If p^o = is chosen then these solutions become vacuum 
ones that are driven by the field. For k = these solutions do not have GR 
counterparts. 

We see from ()3.2.9|) . and figure |SI2l that for < — 3/2 the initial singularity is 
avoided; for a more detailed discussion of this effect, see p?T| . 
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Figure 3.1: The evolution of a and as functions of rj for uj > — | with \rii\ = 
10~^'^GeV~^ , UJ = 40000, and rj^ set so that a(0) = 0. The solid line corresponds to 
771 > and the dashed line to rji < 0. 
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Figure 3.2: The evolution of a and as functions of rj for uj < — | with |?7i| 



UJ 



40000, and 7]2 set so that a'(0) = 0. The solid line corresponds to 



771 > and the dashed line to 771 < 0. 
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3.2.2 Matter Domination 

When considering a matter-dominated universe we could proceed as above and 
determine a set of general solutions to ()3.1.1|) . ()3.1.2p and ()3.1.3|) that at early times 
are described by free scalar-field domination and at late times by matter domination 
(see jnZI), but for our purposes this is unnecessary. For a realistic universe we require 
a period of radiation domination during which primordial nucleosynthesis can occur. 
If the scalar-field-dominated period of the Universe's history were to impinge upon 
the usual matter-dominated period then we would effectively lose the radiation- 
dominated era. For this reason it is sufficient to ignore the free scalar component of 
the general solution and consider only the Machian component. This is equivalent 
to imposing the condition ^ as a ^ 0. With this additional constraint the 
solutions to ()3.1.1|) . ()3.1.2|) and ()3.1.3p . for = and p = 0, are given by |135j 

, > 2+2u , , 2 

a{t)=aj^ and ^(t) = (pj^ (3.2.10) 

where 

Svr pmo _ 2(3 + 2tu) 
~ 3(4 + 3w) 

and a* and 0^, are constants. These solutions can be seen to approach their GR 
counterparts as cx; ^ oo. 



3.2.3 Vacuum Energy Domination 

Similarly, for a vacuum {p = —p) dominated period of expansion we can impose the 
condition ipa? ^ as a ^ to get the power law exact solutions 



a{t) = a^t^2+'^ and 0(t) = 0^^^ (3.2.11) 

where 

^PAo = ^(3 + 2a;)(5 + 6c^) 

and and (p^ are constants. Note that this vacuum stress does not produce a de 
Sitter metric, as in GR. 
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3.2.4 Vacuum Solutions 



In the Brans-Dicke theory it is possible to have spatially flat and positively curved 
exact vacuum solutions, unlike in the general relativistic cosmologies. 

Spatially flat solution 

Assuming solutions of the form (j) cc and a oc and setting a(0) = gives on 
substitution into (jSHHj), ()H.1.2|1 and the k = Brans-Dicke vacuum solutions 



a 



(t) = tla+2(i-v/5G^)-) ^(t)=ct>m(-] . (3.2.12) 



to 



Closed solution 



For the closed region we now follow the method given in ^Hj to flnd expressions for 
a(r) and 0(t). We start by introducing conformal time, rj, deflned by adrj = dr; 
then equation ()3.1.2|) becomes 



a 

This integrates directly to yield 



2 _ 



2 



V3A(2cj + 3)-^/^ (3.2.13) 



where A is a constant. We now introduce the variable y = (ppO? to write ()3.1.3|1 as 

= -4ky^ + ^0,2 a\2uj + 3). (3.2.14) 
Equations (13.2. 13jl and (|3.2.14jl then give 

^ = V3Ay-\2uj + 3)-^/^ and y,l= -Aky^ + A". (3.2.15) 



The solutions of equations ()3.2.15|) . when A; > 0, are given by 

y(r]) = sm{2Vk{r] + 5), (3.2.16) 
2\/k 
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and 

(j){ri) = C tan^ {Vk{ri + B)) (3.2.17) 

where B and C are arbitrary constants. We now fix the conformal time origin by 
setting S = 0, so that y = (ppO^ gives 



sm 



1/2 



air]) oc " ■ (3.2.18) 

tanV5^2^(v^r/) 

3.3 More General Theories 

In order to evaluate more general scalar-tensor theories it is convenient to work in 
the Einstein frame. In this frame, substituting the FRW line-element into the field 
equations gives 

a\ ^ Svr , 1 -9, , , 

a J 3 2 

^ = --{p + 3p + 2r) (3.3.2) 
a 3 

^ + 3-iP = -^a(p - 3p) (3.3.3) 
a 

where over-dots here denote differentiation with respect to i and d{t) = e~^a{t) 
is the scale-factor in the Einstein frame. Defining N = ln(a/ao), Damour and 
Nordtvedt [ZHI write dSZH), and (jTT^ as 

^ ^ -ij" + (2 - 7 - -e)tP' = -(4 - 37)-^, (3.3.4) 



and 



(3-47r^'2)^ ' ' 3 ' 

where e = Sk/Sirpd'^, p = {'j — l)p and ' denotes differentiation with respect to A^. 
Now we consider a coupling parameter of the form 

r(v>) = r(V'oo) + aoim - v-oo) + \p{m - V-oo)' 

so that 

a{t) = ^ + -y={m - (3.3.5) 
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This is of the same form as chosen by Santiago, Kalhgas and Wagoner |167j and 
Damour and Pichon ||22j- Santiago, Kalhgas and Wagoner arrive at this form of a by 
assuming the evolution of the Universe has been close to the GR solutions throughout 
the period from primordial nucleosynthesis to the present. They therefore consider 
themselves justified in performing a Taylor expansion about the asymptotic value of 
t/' and discarding terms of second order or higher. Damour and Pichon consider F 
to be a potential down which t/^ runs. They assume that a particle near a minimum 
of a potential experiences a generically parabolic form for that potential. This leads 
them to consider a quadratic form for F which gives, on differentiation, a as above. 
These two lines of reasoning are, of course, equivalent as the parabolic form of a 
particle near its minimum of potential can be found using a Taylor series. Theories 
with this form of a belong to the attractor class which approach GR at late times, 
if we impose the additional condition = 0. 

3.3.1 Radiation Domination 

For the case of a radiation-dominated flat universe, the general solutions of ()3.3.1|) . 
fnnr^ and (nnn?|l with the choice (EH), are, for u > -3/2: 



167r Ir^il \il -r]2 + \r}i\ 



and, for u < —3/2, 



^ - ^1 = -^/^;^ + tan- ^'(^ - 



16 Ir^il V 47r |?7i| \ \rii 



Here, r/i determines the evolution of the scalar field, ipi is the value it approaches 
asymptotically and r]2 sets the origin of the conformal time coordinate defined by 
a{rj)dri = dt. For uj > —3/2, rji ElZ such that E 71 whilst for u < —3/2, r^i G X 
such that E I and uj eTZ. 
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The corresponding solutions in the Jordan frame are 

0(r7)=exp(-/?(^(r7)-V'oo)') (3.3.6) 

and 




(3.3.7) 



They exhibit the same features as their Brans-Dicke counterparts: at early times 
there is a period of free-scalar-field domination, and at late times they approach 
a (X and =constant. 



3.3.2 Matter Domination 

Solutions during the matter-dominated era are difficult to find because the energy- 
momentum tensor for the matter field is not conserved. However, it is possible to 
obtain an evolution equation for ip. Using ()3.3.4|) . we get 

~ -r + (1 - U^' = -#(V^ - ^oo). (3.3.8) 



(3 - 47rV^'2) ' 3 Att 

For a fiat universe e = and we can solve ()3.3.8j) by making the simplifying assump- 
tion Anip''^ ^ 3. This gives the solution 



^(iV) - = Ae-t(i+v^)^ + Be'ld-v^)^ (3.3.9) 
where A and B are constants of integration. 



3.3.3 Vacuum Domination 

Similarly, for the case of a vacuum-dominated universe the evolution of can be 
approximated using ()3.3.4|) to obtain 

^(AT) - = Ce-i(i+v^)^ + De-i(i-v^)^ (3.3.10) 
where C and D are constants of integration. 
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CHAPTER 4 



Homogeneous and Isotropic Fourth-Order 

Cosmologies 



This chapter is based on the the work of Chfton and Barrow |7()j . 

4.1 Freidmann- Robertson- Walker Cosmologies 

In this chapter we will again be concerned with the idealised homogeneous and 
isotropic space-times described by the FRW metric: 

ds" = -de + a\t) ( — + r^de' + sin^ OdA . (4.1.1) 
\ (1 — kr^) J 

Substituting this metric ansatz into the C = R^^^ field equations ()2.1.15|) . and 
assuming the Universe to be filled with a perfect fiuid of pressure p and density p, 
gives the generalised version of the Friedmann equations 

^1 - 6)R^+^ + 35(1 + 6)R^ I ^ + 3-^ ] - 35(1 - 6^)R^^ = ^{p - 3p) (4.1.2) 



R aR ' ' i?2 2 

. . 5)i?^ + ^ + 35(1 + 5)^9. 
a I a R 



-3-(l + 5)R' + ^ + 35(1 + 5)-^R' = (4.1.3) 
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where, as usual, 

a k , ^ 

i? = 6- + 6— + 6— . 4.1.4 

a 

It can be seen that in the hmit 5 — > these equations reduce to the standard 
Friedmann equations of GR. A study of the vacuum solutions to these equations for 
all k has been made by Schmidt (see the review |171j ) and a quahtative study of 
the perfect-fluid evolution for all k has been made by Carloni et al [58/. Various 
conclusions are also immediate from the general analysis of f{R) Lagrangians made 
in reference [31] by specialising them to the case / = R^~^^. In what follows we shall 
be interested in extracting the physically relevant aspects of the general evolution 
so that observational bounds can be placed on the allowed values of 6. 

Assuming a perfect-fiuid equation of state of the form p = (7 — l)p gives the 
usual conservation equation p oc a~^'^. Substituting this into equations ()4.1.2|1 and 
()4.1.3|) . with k = 0, gives the power-law exact Friedmann solution for 7 7^ 

a{t) = (4.1.5) 

where 

(1 - 25) (2 - 357 - 25^(1 + 37)) = 12ttG{1 - 6)^'^pc (4.1.6) 

and pc is the critical density of the universe. 

Alternatively, if 7 = 0, there exists the de Sitter solution 

a{t) = e"* 

where 



3(1 - 26)n^ = 8nG{l - 6)p,. 

The critical density ()4.1.fjj) is shown graphically, in figure 14.11 in terms of the 
density parameter fio = ^3^2" as a function of S for pressureless dust (7 = 1) and 
black-body radiation (7 = 4/3). It can be seen from the graph that the density of 
matter required for a fiat universe is dramatically reduced for positive 6, or large 
negative 6. In order for the critical density to correspond to a positive matter density 
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Figure 4.1: The density parameter, Qq, as a function of 5. Solid line corresponds to 
pressureless dust and dashed line to black-body radiation. 

we require 6 to lie in the range 



4.2 The Dynamical Systems Approach 

The system of field equations ()4.1.2p and ()4.1.3p have been studied previously using a 
dynamical systems approach by Carloni, Dunsby, Capozziello and Troisi for general 
k We elaborate on their work by studying in detail the spatially fiat, k = 0, 
subspace of solutions. This allows us to draw conclusions about the asymptotic 
solutions of ()4.1.2p and ()4.1.3p when k = and so investigate the stability of the 
power-law exact solution ()4.1.5p and the extent to which it can be considered an 
attractor solution. By restricting to /c = we avoid 'instabilities' associated with 
the curvature which are already present in general relativistic cosmologies. 

In performing this analysis we choose to work in the conformal time coordinate 



V73 + QQu + 9uj^ + 37 
4(1 + 37) 



<6 < 



V73 + QQu + 9uj^ - 37 



(4.1.7) 
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Making the definitions 



R' , a' 

X = — and y = — 
R a 



where a prime indicates differentiation with respect to r, the field equations ()4.1.2|1 
and ()4.1.3j) can be written as the autonomous set of first order equations 

, _ 2 - ^(37 - 2) _S^_ (4-^(37-2))xt/ _ 2(1-%^ 
" 6^(1 + 6) 2 26 6^ ^ ' 

.--^.^(2 + 3.).. + i^±^i|p^. (4.2.3) 

These coordinate definitions are closely related to those chosen by Holden and 
Wands |104j for their phase-plane analysis of Brans-Dicke cosmologies and allow us 
to proceed in a similar fashion. 

4.2.1 Locating the Critical Points 

The critical points at finite distances in the system of equations ()4.2.2|) and ()4.2.3|) 
are located at 

4-37 1 

Xi2 = ± — and yi2 = ± , (4.2.4) 

V(l + 5)(5-37) ' ^(l + 5)(5-37) 

and at 



and 



The critical points {xi,yi) and {x2,y2) do not exist for 7 = 5/3. 

The exact form of a{t) at these critical points, and the stability of these solutions, 
can be easily deduced. At the critical point (xj, yi) the forms of a(r) and R{t) are 
given by 

a(r) = aoe^'" and R{t) = R^e^^^ (4.2.7) 
where ao and -Rq are constants of integration. In terms of r the perfect-fiuid conser- 



X3 4 = T , ^ = (4.2.5) 

v/(l + 5)(2-3(57-52(2 + 67)) 



y34 = ± — . 4.2.6) 

'2 - 357 - 52(2 + 67) ^ ' 
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vation equation can be integrated to give 

p = poe-^^y'\ 

where po is another positive constant. Substituting into the definition of r now gives 

dr cx e-5^^'"-i^'>"c^t 

or, integrating, 

t-tooc 3 L^et^^»"+^'^ (4.2.8) 

It can now be seen that if S'yyi + 6xi > then t — > cx) as r — > 00 and t — > to 
as r — ^ —00. Conversely, if 3'yyi + 6xi < then t ^ to as t ^ 00 and t — ^ —00 as 
r —00. 

In terms of t time the sokitions corresponding to the critical points at finite 
distances can now be written as 

a(t) oc (t - to)'^"»+*"« and oc (t - to)^^^^- 

The critical points 1 and 2 can now been seen to correspond to a oc and the 
points 3 and 4 correspond to ()4.1.5j) . 

In order to analyse the behaviour of the solutions as they approach infinity it is 
convenient to transform to the polar coordinates 

X = fcoscj) 
y = f sin0. 

The infinite phase plane can then be compacted into a finite size by introducing the 
coordinate 

r 

r = . 

1 + f 
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The equations ()4.2.2j) and ()4.2.3j) then become 
-1 



4(1 - 2r)(5(l + 6) sin0 - (2 - 5(37 - 2)) cos< 



452(1 + s) 

- r\{Q + 86 + 35^ + 6'^- 126^) cos0 + (1 + 5) (2 - 25 - 5^ - 25^) cos ; 

- 25(3 - 5(37 - 2) + 3 cos 20) sin (p)] (4.2.9) 

and 



-1 



(25(1 + 5) COS0 + 2(2 - 5(37 - 2)) sin0)(l - 2r) 



252(1 + 5)(1 -r)r 

- (^5(1 + 5) cos 0(1 - 3 cos 20) - 4 sin + 4(1 - 5)^ sin^ 

+ 25(37- 2 + 5(l + 5)(l + 25)cos2 0)sin0)r2 1. (4.2.10) 



In the hmit r — > 1 {f oo) it can be seen that critical points at infinity satisfy 

sin0j(5cos0j + sin 0j) (5(1 + 25) cos0i + 2(1 — 5) sin0j) = 
and so are located at 

</'5,(6) = (+7r) (4.2.11) 
07,(8) =tan-i(-5) (+7r) (4.2.12) 

09,(10) = tan-1 (-f^fZ^) (+^)- (4-2-13) 
The form of a{t) can now be calculated for each of these critical points by proceeding 
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as Holden and Wands j|in4j . Firstly, as r — 1 equation ()4.2.9|1 approaches 
^' ^ TT^ I -^(l + 25(37 - 2)) sin (pi - 35 sin 30^ 

- (2 - (5(2 + 5)) cos 0i + (2 - 5(2 + 5 + 25^)) cos 3( 

which allows the integral 

r-1 = /(</.,) (r -To) 

where the constant of integration, tq, has been set so that r — >• 1 as r — > Tq. Now 
the definition of x allows us to write 

R' r /(0j)(r - To) + 1 cos0i 



cos 0i = r COS v^i 



R (l-r) /(0,)(r-ro) /(0,)(r-ro) 

as r — > tq. Integrating this it can be seen that 



-Roc|r — ro| ^'-'^^^ as r ^ 1. 

Similarly, 

sin (p^ 

a oc |r — tqI •''('^i) as r — 1. 
The definition of r ()4.2.ip now gives 

37 sin (tn |5 cos 0^ 

c/r OC |r - To I =^ fi-f.)^^7m dt 

which integrates to 

t-tooc-^^|r-ro| W (4.2.14) 

where 

37sin0i + 5cos0i - 2/(0i) 



2 

The location of critical points at infinity can now be written in terms of t as the 
power-law solutions 

cos (j)^ sin 

R{t) oc (t - to)^ and a{t) oc (t - to)^- (4.2.15) 
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Direct substitution of the critical points (|4.2.njl . (j4.2.T4 and (|4.2.i:-ijl into 
(|4.2.15jl gives 



05,6 (^) constant 
a7,8(t) Vt - to 

6(1+2S) 

agM^) (t - to) 

as r ^ 1. Moreover, it can be seen from ()4.2.14j) that as r — > 1 and t ^ tq so t ^ 
as long as F {(pi) / f < 0, as is the case for the stationary points considered here 
(as long as the value of 6 lies within the range given by ()4.1.7p ). 

The exact forms of a{t) at all the critical points are summarised in the table I^^Tl 



Critical point 


a(t) 


1, 2, 7 and 8 

3 and 4 
5 and 6 

9 and 10 


i 

2(1+5) 
t 37 

constant 

5(1+2*) 
t (l-*) 



Table 4.1: The form of a{t) at critical points of the phase plane. 
4.2.2 Stability of the Critical Points 

The stability of the critical points at finite distances can be established by perturbing 
X and y as 

x{r) = Xi + u{r) and y{r) = yi + v{r) (4.2.16) 
and checking the sign of the eigenvalues, A,, of the linearised equations 

u' = XiU and v' = XiV. 

Substituting ()4.2.16|) into equations ()4.2.2p and ()4.2.3|) and linearising in u and 
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V gives 



2 V 2 5 

The eigenvalues Aj are therefore the roots of the quadratic equation 

+ + C = 

where 

S = -^(2 + 5)x,-^(37-2)y, 

C = -^(2 + 35)x^ - (2 + 5(1 + ?>5))xiyi + -^(8 - 85 - 67(1 + 35) + 957')2/- • 
2 2o 

If i? > and C > then both values of Aj are negative, and we have a stable critical 
point. If i? < and C > both values of Aj are positive, and the critical point is 
unstable to perturbations. C < gives a saddle-point. 

For points 1 (upper branch) and 2 (lower branch) this gives 

*^(l+i)(5-37) 6(1 + S) 

and for points 3 (upper branch) and 4 (lower branch) 

" ■ , 3(2 + 2.-7(1 + 2.)) c=<l±i^. 

v/2(l + 5)(2 - 357 - 252(1 + 37)) 5(1 + 5) 

The stability of the critical points at finite distances for a universe filled with 
pressureless dust, for various different values of 5, are given in table 14.21 and for 
black-body radiation are given in table 14.31 

Values of 5 < _ V^'^+^^t^+^t' g^^id 5 > "^^^^^^^3^^ ^"^ have not been considered 
here as they lead to negative values of po for the solution ()4.1.5|1 . (The reader may 
note the difference here between the range of 5 for which point 3 is a stable attractor 
compared with the analysis of Carloni et. al.). 
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Critical 
point 


B 


c 


V^+3 ^ 1 
16 ^ " ^ 4 


-l<5<0 


0<6< 

16 


1 


(1+26) 


(l+4d-) 


Saddle 


Stable 


Saddle 


5^2(1+5) 


5(1+5) 


2 


(1+25) 


(1+45) 


Saddle 


Unstable 


Saddle 


5y/2(l+5) 
3 


5(1+5) 


3 


(1+45) 


Stable 


Saddle 


Stable 


a/2(1+5)(2-35-8<52) 


5(1+5) 


4 


3 


(1+45) 


Unstable 


Saddle 


Unstable 


A/2(l+<5)(2-3i-852) 


5(1+5) 



Table 4.2: The stability of critical points at finite distances for a universe filled with 
pressureless dust. 



Critical point 


B 


c 


v/6+1 < 5 < Ve-i 

5 5 


1 
2 
3 

4 


3 


4 


Saddle 
Saddle 
Stable 

Unstable 


3 


(l+<5) 
4 


VT+5 
(1-5) 


(1+5) 
4 


a/(1+5)(1-25-552) 
(1-5) 


(1+5) 
4 


v/(l+5)(l-25-552) 


(1+5) 



Table 4.3: The stability of critical points at finite distances for a universe filled with 
black-body radiation 

Point 3 lies in the y > region and so corresponds to the expanding power-law 
solution ()4.1.5j) . It can be seen from the table above that this solution is stable 
for certain ranges of 6 and a saddle-point for others. In contrast, point 4, the 
contracting power-law solution ()4.1.5|) . is unstable or a saddle-point. The nature of 
the stability of these points and the trajectories which are attracted towards them 
will be explained further in the next subsection. 

A similar analysis can be performed for the critical points at infinity. This time 
only the variable will be perturbed as 

= 0, + g(t). (4.2.17) 

The conditions for stability of the critical points are now that r' > and the 
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eigenvalue n of the linearised equation q' = fiq satisfies /i < in the hmit r — 1. 
If both of these conditions are satisfied then the point is a stable attractor; if only 
one is satisfied the point is a saddle; and if neither are satisfied then the point is 
repulsive. 

Substituting ()4.2.17j) into ()4.2.1()j) and linearising in q(t) gives, in the limit r —>■ 1, 



452(1 -r) 



(6(1 -5) + 5^(1 + 25)) cos 0i 

- 3(2(1 - 5) - 5^(1 + 25)) cos 30^ - 35 sin (pi + 95 sin 30 )g 



= fiq. 



The sign of r' in the limit r — > 1 can be read off from ()4.2.9|1 . The stability properties 

of each of the stationary points at infinity can now be summarised in the table 03 

A/18+387+972-37 



where A^i 



A/18+387+972+37 
4(1+37) 



and N2 



4(1+37) 



Critical point 


-A^i < 5 < 


-i < 5 < 


0<5<i 


i<S<N2 


5 


Stable 


Saddle 


Unstable 


Unstable 


6 


Unstable 


Saddle 


Stable 


Stable 


7 


Unstable 


Unstable 


Saddle 


Stable 


8 


Stable 


Stable 


Saddle 


Unstable 


9 


Saddle 


Stable 


Stable 


Saddle 


10 


Saddle 


Unstable 


Unstable 


Saddle 



Table 4.4: The stability of critical points at infinity. 



4.2.3 Illustration of the Phase Plane 

Some representative illustrations of the phase plane are now presented. Firstly, the 
compactified phase plane for a universe filled with pressureless dust, 7 = 1, and a 
value of 5 = 0.1 is shown in figure 

Figure 14.21 is seen to be split into three separate regions labelled I, II and III. 
The boundaries between these regions are the sub- manifolds R = 0. As pointed out 
by Carloni et. al. the plane i? = is an invariant sub-manifold of the phase space 
through which trajectories cannot pass. 
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Figure 4.2: Phase plane of cosmological solutions for 7 = 1 and 5 = 0.1. r is used 
as the radial coordinate and (j) is measured anti-clockwise from y = 0, a; > 0. 
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The equation for i? in a FRW universe, ()4.1.4j) . can now be rewritten as 

= ^((1 + 5)y' + 5{l + 5)xy - 1). (4.2.18) 

This shows that the boundary i? = is given in terms of x and ?/ by (1 + 5)?/^ + 
(5(1 + (5) xy — 1 = and that in region I the sign of R must be opposite to the sign 
of 5 in order to have a positive p. Similarly, in regions II and III, R must have the 
same sign as 5 in order to ensure a positive p. 

It can be seen that regions II and III are symmetric under a rotation of vr and a 
reversal of the direction of the trajectories. As region II is exclusively in the semi- 
circle ?/ ^ all trajectories confined to this region correspond to eternally expanding 
(or expanding and asymptotically static) universes. Similarly, region III is confined 
to the semi-circle y ^ and so all trajectories confined to this region correspond to 
eternally contracting (or contracting and asymptotically static) universes. Region 
I, however, spans the y = Q plane and so can have trajectories which correspond to 
universes with both expanding and contracting phases. In fact, it can be seen from 
figure W?^ that, for 5 = 0.1 all trajectories in region I are initially expanding and 
eventually contracting. 

It can be seen from figure 14.21 that in region I the only stable attractors are, at 
early times, the expanding point 10 and at late times the contracting point 9. (By 
'attractors at early times' we mean the critical points which are approached if the 
trajectories are followed backwards in time). Both of these points correspond to the 
solution 

c(l+2S) 

a(xt (1-*) 

which describes a slow evolution independent of the matter content of the universe. 
Notably, region I only has stable attractor points, at both early and late times, which 
are reached in a finite time as t —inconstant. In region II the only stable attractors 
can be seen to be the static point 5 at some early finite time, to, and the expanding 
matter-driven expansion described by point 3 as t — 5> oo. Conversely, in region III 
the only stable attractors are the contracting point 4 for t — > — oo and the static 
point 6 for t ^ to- 

Figure 14.31 shows the compactified phase plane for a universe containing pres- 



46 



sureless dust and having 6 = —0.1. Figure lOl is split into three separate regions 
in a similar way to figure with the boundary between the regions again corre- 
sponding to -R = and is given in terms of x and y by ()4.2.18|) . Regions II and 
III again correspond to expanding and contacting solutions, respectively. Region I, 
still has point 10 as the early-time attractor and point 9 as the late-time attractor, 
but now has all trajectories initially contracting and eventually expanding. Stable 
attractors in Region I are still reached in a finite time. Region II now has point 7 
as an early-time stable attractor solution and point 1 as a late-time stable attractor 
solution, corresponding to a ^ as t ^ oo. Point 3, which was the stable attractor 
at late times when 6 = 0.1, is now no longer located in Region II and can instead be 
located in region I where it is now a saddle-point in the phase plane. Interestingly, 
the value of 6 for which point 3 ceases to behave as a stable attractor {6 = 0) is 
exactly the same value of 6 at which the point moves from region II into region I; 
so as long as point 3 can be located in region I, it is the late-time stable attractor 
solution and as soon as it moves into region I it becomes a saddle-point. At this 
same value of 6, point 1 ceases to be a saddle-point and becomes the late-time stable 
attractor for region II, so that region II always has a stable late-time attractor as 
t oo. Region III behaves in a similar way to the description given for region II 
above, under a rotation of vr and with the directions of the trajectories reversed. 

Phase planes diagrams for 7 = 1 with values of S other than 0.1 and —0.1, but 
still within the range 

v'ls + 387 + 972 + 37 v^is + 387 + 972 - 37 

4(1 + 37) ^ ^ 4(1 + 37) ' 

look qualitatively similar to those above with some of the attractor properties of the 
critical points being exchanged as they pass each other. In particular, for 5 < — ^ 
point 3 returns to region II and once again becomes the stable late-time attractor 
for trajectories in that region. The points that are the stable attractors for any 
particular value of 6 can be read off from the tables in the previous section. 

Universes filled with perfect fluid black-body radiation also retain qualitatively 
similar phase-plane diagrams to the ones above; with the notable difference that the 
point 3 is always located in region II and is always the late-time stable attractor of 
that region. This can be seen directly from the Ricci scalar for the solution ()4.1.5|) 
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Figure 4.3: Phase plane of cosmological solutions for 7 = 1 and S = —0.1. r is used 
as the radial coordinate and (j) is measured anti-clockwise from y = 0, a; > 0. 
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which is given by 

3^(1 + 5) 

and can be seen to have the same sign as 6, for 6 > —1, and so is always found in 
region II. 

For a spatially-flat, expanding FRW universe containing black-body radiation 
we therefore have that ()4.1.5p is the generic attractor as t ^ oo. Similarly, for a 
spatially-flat, expanding, matter-dominated FRW universe ()4.1.5j) is the attractor 
solution as t — > oo; except when < 5 < 0, in which case it is point 1 (a oc t^). 

2 

If we require a stable period of matter domination, during which a{t) ~ ta, we 
therefore have the theoretical constraint 6 > (or 5 < —\)- Such a period is neces- 
sary for structure to form through gravitational collapse in the post-recombination 
era of the Universe's expansion. 

The effect of a non-zero curvature. A; 7^ 0, on the cosmological dynamics is 
similar to the general relativistic case. The role of negative curvature {k = —1) can 
be deduced by noting that its effect is similar to that of a fluid with 7 = 2/3. The 
solution ()4.1.5j) is unstable to any perturbation away from flatness and will diverge 
away from A; = as t — 00. This is usually referred to as the 'flatness problem' and 
can be seen to exist in this theory from the analysis of Carloni et. al. |58j . 
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CHAPTER 5 



Other Homogeneous Fourth-Order Cosmologies 



We will consider in this chapter the conditions required for gravity theories that are 
derived from Lagrangians that are functions of the scalar curvature and Ricci and 
Riemann curvature invariants to possess solutions which are homogeneous space- 
times other than the standard FRW solutions. We will firstly consider space-times 
of the Godel jHS [UHl E], Einstein static [21], and de Sitter [102 types. The 
Godel solution admits a G5 group of motions, the Einstein static solution has a G7 
and the de Sitter solution has the maximal Giq. In the Godel case we determine 
the conditions for the existence or non-existence of closed time-like curves. These 
three space times share the common property that the three curvature invariants 
X = R, Y = R^^Rab and Z = R°'^^'^Rahcd are constant for each of them. This gives 
the fourth-order field equations a more simple form, allowing solutions to be found 
more readily. 

We will also consider the vacuum Kasner solutions ^111] and their fluid-filled 
counterparts that form type I of the Bianchi classification of three-dimensional ho- 
mogeneous spaces. These universes have a G3 group of motions acting transitively on 
a three dimensional subspace and are geometrically special. In vacuum (or perfect- 
fluid) cases they are defined by just one (or two) free constant(s) compared to the 
four (eight) that specify the most general spatially homogeneous vacuum (or per- 
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feet fluid) solutions. However, they have proved to provide an excellent dynamical 
description of the evolution of the most general models over finite time intervals. In 
GR the chaotic vacuum Mixmaster universe of Bianchi Type IX undergoes a infinite 
sequence of chaotic space-time oscillations on approach to its initial or final singular- 
ities which is well approximated by a sequence of different Kasner epochs which form 
a Poincare return mapping for the chaotic dynamical system [IHl 1126^ |H1 El EH I159j . 

This chapter is based on the work of Clifton and Barrow [71], |72] and Barrow 
and Clifton 

5.1 Godel, Einstein and de-Sitter Universes 

These three homogeneous space-times, and the investigations of their stability, have 
played a central role in our understanding of the dynamics of GR and in the pos- 
sible astrophysical consequences of general relativistic cosmologies. Higher-order 
modifications of GR are of importance in assessing the corrections that might be 
introduced in high- curvature environments and can also be of use in explaining the 
late-time acceleration of the Universe. Furthermore, investigations of these theories 
allows for an evaluation of the special nature of GR itself. Here we study the condi- 
tions under which the Godel, Einstein static, and de Sitter universes exist in a wide 
class of non-linear gravity theories. 

Considerable simplification of the fourth-order field equations ()2.1.12j) occurs if 
the three curvature scalars X, Y and Z are constant. In this case the fourth-order 
field equations, ()2.1.12|1 . reduce to 



P""' = -\fg'''+R''\fx-2Rfz)-2R'''''RcdUY+^fz)+\g''\X^-AY+Z)fz (5.1.1) 



where use has been made of the identities (7^1 



R 



:abcd 



-DR"" + R 



lac d 



■,bc — 



c 



R 



lac d 



1 c 




idc 



1 
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and 



2RR 



ab 



1 



AR 



acdb 



1 



Red + AR\R'' + -g'^'R'. 



Equation ()5.1.1|) is only of second order and is therefore a significant simplification 
of the original system of equations. We will investigate the Godel, Einstein static 
and de Sitter universes solutions of these field equations. 



5.1.1 Godel Universes 

The Godel universe 11481 E] is a rotating homogeneous space-time given by the 
line-element 



ds^ = ~{dt + C(r)#)2 + D'^{r)di/j'^ + dr^ + dz^ 
= -dt^ - 2C{r)dtd^ + G{r)dip'^ + dr^ + dz^ 



(5.1.2) 



where 



C{r) 

D{r) 
G{r) 



4Q , 2 

sinn I 1 



\ 2 J 



— sinh (mr) 
m 

D\r)-C\r) 
sinh^ I 



4 . , 2 /mr\ 



\ 2 J 



1 + ( 1 _ 1^ ) sinh" ( 



2 fmr 



V 2 



and m and VL are constants, and Vl controls the rate of rotation. The existence of 
Godel universes in one particular /(X, F) theory has been previously studied by 
Accioly jl] where a solution was found in vacuum. We extend this analysis to the 
more general class of theories above and to universes filled with matter fields. 

The Godel universe is of particular theoretical interest as it allows the possibility 
of closed time-like curves, and hence time-travel The condition required to 

avoid the existence of closed time-like curves is |158| EEI 



G{r) > for r > 
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or 

m2^4^]l (5.1.3) 

By investigating the existence of Godel universes in this general class of gravity 
theories we will also be able to determine those theories in which the condition 
()5.1.3|1 is satisfied. Consequently, we will be able to determine those theories in 
which time travel is a theoretical possibility. 

It is convenient to work in the non-holonomic basis defined by 

e^o = e(2)2 = 6(^)3 = 1 

= D{r) 
e(o)i = C(r). 

The Godel line-element fj5.1.2j) then becomes 

where the one-forms 6"^^^ are given by = e''^\dx"' (capital Latin letters denote 
tetrad indices and lower case Latin letters denote space-time indices). The inverses 
of e^^\ can be calculated from the relations ~ ^'^^\b) and the non-zero 

elements of the Riemann tensor in this basis are then given by 

^(o){i)(o)(i) = Rmmm = 

-R(i){2)(i){2) = 3ri^ — m^. 

The perfect-fiuid energy-momentum tensor is defined in the usual way with re- 
spect to the comoving 4-velocity U"" = (—1,0,0,0) and its covariant counterpart 
Ua = (1, 0, C(r), 0) such that its non-zero components in the non-holonomic basis 
are given by 

T{o){Q) = p and T(^i){i) = ^(2)(2) = ^(3)(3) = P- 
The field equations, ()2.1.13j) . for this space-time can then be manipulated into 
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the form 

A-fp=^/ (5.1.4) 

= {2VL^ - m^){fx - 4[^]2 - m^]fz) + 2(/y + Afz){8Q^ - 5Q^m^ + m^) (5.1.5) 

|(p + J9) = 2n\fx - - ^Vz) + 4(/y + 4fz)n\2n' - m"). (5.1.6) 

Solving the field equations has now been reduced to solving these three algebraic 
relations for some specified /(X, Y, Z). 

1 f = f{X) 

For p + p we see from ()5.1.6p that fx 7^ 0. From equation ()5.1.5|) it can then be 
seen that = 2fi^, as in general relativity. Therefore for any theory of the type 
/ = f{X) the inequality ()5.1.3|1 is not satisfied and closed time- like curves exist, 
when p + p ^ 0. 

It now remains to investigate the case p + p = 0. It can immediately be seen 
from ()5.1.6p that we must have /x = in order for a solution to exist. This sets the 
relation between m and Q and automatically satisfies equation ()5.1.5|) . The required 
value of A can then be read off from equation ()5.1.4|1 . The condition /x = will 
now be investigated for a variety of specific theories. 

la f = f(X)=X + aX^ 

Theories of this kind have been much studied [3^1 1113j as they have a number of 
interesting properties, not least of which is that they display divergences which are 
normalisable at the one loop level |181j . 

In these theories, the condition /x = is equivalent to 

ni^ = n^ + — . 

4a 

The condition under which this theory then satisfies the inequality ()5.1.3p . and hence 
does not admit closed time-like curves, is ^ (12a)~^. Therefore, for any given 
theory of this kind, defined only by a choice of the constant a, there is a range of 
values of fl for which closed time- like curves do not exist, when a > and p + p = 0. 
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However, when a < this condition is never satisfied and closed time-hke curves 
exist for all values of Q. 

lb f = f(X) = X + ^ 

Theories of this type have generated considerable interest as they introduce cosmo- 
logical effects at late times, when R is small, which may be able to mimic the effects 
of dark energy on the Hubble flow j^Hl I14()| lllfflj . The square in the factor is 
introduced here as these theories require a positive value for this coefficient in order 
for the field equations to have a solution. 

When /x = we have the two possible relations 

2 

where a is the positive real root of a^. The upper branch of this solution then allows 
the condition ()5.1.3|) to be satisfied if QQ"^ ^ a; or, for the lower branch, if ^ —a. 
For any particular value of a > 0, the upper branch always admits a range of fl 
for which closed time-like curves do not exist. For the lower branch, however, the 
inequality is never satisfied and closed time-like curves are permitted for any 

value of Q. 

Ic f = f(X) = \X\^+' 

This scale-invariant class of theories is of interest as its particularly simple form 
allows a number of physically relevant exact solutions to be found [ZHl IIHI EUl IHE] ■ 
In order for solutions of the Godel type to exist in these theories we must impose 
upon 6 the constraint 6^0. 

The condition /x = now gives 

Evidently, the condition ()5.1.3p is never satisfied in this case: Godel solutions always 
exist and closed time-like curves are permitted for any value of Q. 
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2 f = fiX,Y,Z) 

In this general case, equations ()5.1.5|) and ()5.1.6p can be manipulated into the form 
/x-4(n^-,n.^)/.= ' ^ X(P + P) (5.1.7) 

when 7^ 4^2^. The special case = 4^2^ gives 

p + p = (5.1.9) 
/x = 41]2(/y + /z). (5.1.10) 

First we consider p+p > 0. It is clear that = 4^2^ is not a solution in this case. 
Equations ()5.1.7|) and ()5.1.8|) show that in order to satisfy the inequality ()5.1.3p . and 
avoid the existence of closed time-like curves, the following two conditions must be 
satisfied simultaneously 

fx + 4{m^-n^)fz<0 (5.1.11) 
/y + 4/z<0. (5.1.12) 

For p + p < the inequalities must be reversed in these two equations. It is now 
clear that if it is possible to construct a theory which has a solution of the Godel 
type without closed time-like curves for a non-zero p + p of a given sign, then this 
theory cannot ensure the non-existence of closed time-like curves for the opposite 
sign of p + p. 

It remains to investigate the case p + p = 0. In order to have a solution in this 
case we require that the two equations 

fx-^n^-m^)fz = (5.1.13) 
fy + 4fz = (5.1.14) 

are simultaneously satisfied, for 7^ 4^2^, or 

fx = 4n\fY + fz) (5.1.15) 
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for = 4f2^. Any theory which satisfies therefore, does not allow the 

existence of closed time-like curves when p + p = 0. 
We now illustrate these considerations by example. 

2a f = X + aX^ + /3Y + -fZ 

Using the identity If 

- — (y^{X^ — AY + Z)\ = pure divergence 

it is possible to rewrite / as / = X + aX^ + (3Y where d = a — 7 and (3 = (3 + 47. 
This is the theory considered by Accioly jT]. 

For the case p + p > 0, the inequalities ()5.1.11|) and ()5.1.12p then become 

/3<0 and (5.1.16) 

The first of these inequalities can be satisfied trivially. The second can be only be 
satisfied for all if has some non-zero minimum value. From the field equations 
fl5.1.5|) and ()5.1.6|) we can obtain 

(2^2 _ ^2^^^ ^ 4^[^2 _ ^2]) ^ -2/3(80^ - 5nW + m^) 

which gives two solutions for m as a function of fi, a and (3. Substituting either of 
these values of m into R = 2(f2^ — m^) then gives that R —>■ as ^ It can 
now be seen that the second equality in ()5.1.16|1 cannot be satisfied for all if d is 
finite. It is, therefore, not possible to construct a theory of this type which excludes 
the possibility of closed time-like curves for all Q, when p + p ^ 0. 

For p + p < the inequalities in ()5.1.16|) must be reversed. In this case f3 > 
and d > allows a range of Q in which closed time-like curves are not permitted 
and P > and d < does not allow closed time-like curves for any values of Q. 

The case p + p = was studied by Accioly P, who found that equation ()5.1.15|1 
is satisfied if 

4(]2 ^ ^2 ^ 

(3d + p) 

Hence closed time-like curves do not exist in Godel universes for these theories, when 
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p + p = 0. 



2b f = ciX + I3Y 

This class of theories is introduced as an example which excludes the possibility of 
closed time- like curves for p + p > 0. The inequalities ()5.1.11|) and ()5.1.12j) in this 
case are 

a < and /5 < 0, 

which can be trivially satisfied. For p + p = the only non-trivial solution is given 
by = 4f2^, the value of Q then being given in terms of a and (3 by ()5.1.15|) . 

This example shows explicitly that it is possible to construct a theory in which 
closed time-like curves do not occur in Godel universes when p + p ^ (though 
we do not consider it as physically viable as a is required to have the 'wrong' sign 

Table 15.11 summarises the results found in this section. 

5.1.2 Einstein Static Universes 

The Einstein static universe is a homogeneous and isotropic space-time with line- 
element 

ds^ = -dt" + " , + r\d9^ + sin^ edcf)^). 
(1 — kr^) 

Here, k parametrizes the curvature of the space-like slices orthogonal to t and the 
scale-factor has been rescaled to 1. For a universe containing pressureless dust the 
field equations ()2.1.13p can now be written as 

XP = ^kfx + lQk\fy + fz) (5.1.17) 
A = i/ - 2kfx - Sk'ify + fz). (5.1.18) 

It can now be seen immediately that solutions exist for an Einstein static universe for 
any f{X, Y, Z) that is differentiable in all its arguments. The corresponding values 
of p and A are simply read off from equations ()5.1.17|) and ()5.1.18|) . It remains to be 
studied under what circumstances these solutions are stable. The investigation by 
Barrow, Ellis, Maartens, and Tsagas j22j shows that this is an issue that depends 
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Table 5.1: A summary of the conditions under which closed time-like curves can exist 
in Godel universes, for various different gravitational theories, defined by f{X, Y, Z): 
•/ denotes their existence for all values of and x denotes that they are not allowed 
for any value of 1). The symbol // x means that closed time-like curves are allowed 
to exist for some restricted range of only; the ranges are given in the main text. 

upon the material content and the equation of state of matter in a delicate fashion. In 
GR, there is first-order stability against density perturbations when the sound speed 
exceeds a critical value (l/\/5 of the speed of light) because the Jeans length exceeds 
the size of the universe j9H |2Zj • However, there is instability against homogeneous 
gravitational- wave modes of Bianchi IX type [2Zj • In general, we expect a universe 
with compact space sections and Killing vectors to display linearisation instability 
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5.1.3 de Sitter Universes 



The line element for the maximally symmetric de Sitter vacuum universe can be 
written as 

ds' = -dt" + ^^\dr^ + r\de^ + sin^ ed<P^)) 

where A is the cosmological constant in GR. For de Sitter space-time all components 
of the Riemann and Ricci tensors can be written in terms of the Ricci scalar using 
the equations 

Rabcd = —-^R{9adgbc — QacQbd) 
Rab = -^Rgab 

where the Ricci scalar is i? = 4A. The field equations ()2.1.13|) can now be reduced 
to the single equation 

^/ - A = A/x + 2AVy + ^A'fz- (5.1.19) 

This equation must be satisfied by / if the de Sitter universe is to be a solution in 
any particular gravitational theory. This problem reduces to that studied by Barrow 
and Ottewill [36J in the case / = f{X). This result establishes the situations where 
inflation of de Sitter sort can arise from higher-order corrections to the gravitational 
Lagrangian. In the case where / = f{X) alone it is appreciated that the resulting 
theory is conformally equivalent to general relativity plus a scalar field P3lll22j with 
an asymmetric exponential potential and so either de Sitter or power-law inflation 
is possible. Our results establish when de Sitter inflation is possible in situations 
where the other invariants, Y and Z contribute to the Lagrangian, and the conformal 
equivalence with GR is broken. The stability of the de Sitter solutions in these 
theories will be studied elsewhere. 

5.2 Anisotropic Cosmologies 

The way in which the Kasner metric has played a central role in the elucidation of 
the existence and structure of anisotropic cosmological models and their singularities 
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in GR makes it the obvious starting point for an extension of that understanding 
to cosmological solutions of higher-order gravity theories. We will determine the 
conditions for the existence of Kasner-like solutions and find their exact forms. In 
some cases these solutions are required to be isotropic and correspond to exact 
FRW vacuum solutions with zero spatial curvature. By studying gravity theories 
whose Lagrangians are derived from arbitrary powers of the curvature invariants we 
are able to find simple exact solutions. Past studies have usually focussed on the 
addition of higher-order curvature terms to the Einstein-Hilbert Lagrangian. This 
results in enormous algebraic complexity and exact solutions cannot be found. The 
results presented here provide a tractable route into understanding the behaviours 
of anisotropic cosmological models in situations where the higher-order curvature 
corrections are expected to dominate. The solutions we present are vacuum solu- 
tions but we provide a simple analysis which determines when the introduction of 
perfect fiuids with non-comoving velocities has a neghgible effect on the cosmological 
evolution at early times. 

5.2.1 Field Equations 

We will be looking for spatially homogeneous, vacuum solutions of Bianchi type I, 
described by the line-element 



where pi, p2 and ps are constants to be determined. For the special case Pi = P2 = Ps 
these solutions correspond to spatially fiat FRW metrics, which will be found to exist 
for various higher-order theories both in vacuum (where none exist in GR) and in 
the presence of a perfect fiuid. 

We can determine the number of independently arbitrary functions of three space 
variables that will characterise the general vacuum solution of the field equations on 
a Gauchy surface of constant time in these higher-order theories. The field equations 
are in general fourth order in time; so if we choose a synchronous reference system 
then we need 6 functions each for the symmetric 3x3 tensors Qap, Qap, 9 a/3 a-nd 9 ap- 
This gives 24 functions, but they may be reduced to 20 by using the 4 constraint 



3 




(5.2.1) 



i=l 
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equations, and finally again to 16 by using the 4 coordinate covariances. If a perfect 
fluid were included as a matter source the final number would rise to 20 due to 
the inclusion of the density and 3 non-comoving velocity components in the initial 
data count. We note that in GR the general vacuum solution is prescribed by 3 
arbitrary functions of 3 space variables, and that in GR the Kasner vacuum solution 
is prescribed by one free constant, as the three pi satisfy two algebraic constraints. 

5.2.2 Exact Solutions 

For a realistic theory we should expect the dominant term of the analytic function 
f{X,Y,Z) to be of the Einstein-Hilbert form in the Newtonian limit. However, 
there is no reason to expect such a term to dominate in the high curvature limit. In 
fact, this is the limit in which quantum corrections should become dominant. We 
therefore allow the dominant term in a power series expansion of f{X, Y, Z) to be 
of the form i?", {RabR'^^Y {RabcdR"''^'"'')^ approach to the singularity, where n 
is a constant. These three different cases will be investigated separately below. 

f = R^ 

Substituting / = f{X) = R" into the fourth-order field equations ()2.1.12j) together 
with the line-element ()5.2.H1 gives the two independent algebraic constraints, 

(2^2 - 4n + 3)P +{n- 2)Q - 3{n - l){2n - 1) = 

and 

2{n^ - 1)P + (n - 2)(P2 - Q) - 3{n - l)\2n - 1) = 0, 
where we have defined 

P = J2pi and Q = J2pI (5-2-2) 

i i 

The constraint equations have two classes of solution. The first is given by 

^_ 3(n-l)(2n-l) 
{2-n) 
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which is only solved by the isotropic solution pi = P2 = Ps = -P/3. This is the 
zero-curvature isotropic vacuum universe found by Bleyer and Schmidt |17H IHH] . It 
can be seen that these isotropic cosmologies are valid for all n 7^ 2 and correspond 
to expanding universes for n<l/2orl<n<2 and to collapsing universes for 
1/2 < n < 1 or 2 < n. 

The second class of solutions to ()5.2.2|1 is given by 

P = 2n-1 (5.2.3) 
Q = (2n - 3)(1 - 2n). 

This class represents a generalisation of the Kasner metric of GR. It can be shown 
that it only exists for n in the range 1/2 ^ n ^ 5/4 and that the values of the 
constants Pi must then lie within the ranges 

2n - 1 - 2v/(2n- l)(5-4n) ^ 3j9i ^ 2n - 1 - v^(2n - l)(5-4n), 
2n-l- v/(2n-l)(5-4n) ^ 3p2^2n-l + ^ {2n - l)(5-4n), 
2n - 1 + v/(2n-l)(5-4n) ^ 3j93 ^ 2n - 1 + 2v/(2n - 1)(5 - 4n), 

where it has been assumed without loss of generality that p\ ^ p2 ^ Ps- These ranges 
are shown on figure 5.1 and can be read off by drawing a horizontal line of constant n 
on the plot. The four points at which the horizontal line crosses the two closed curves 
gives the four boundary values for the allowed ranges of the constants Pi. It can 
be seen that the points at which the curves cross the abscissa, which corresponds 
to = 1, give the boundary values — 1/3 ^ pi ^ ^ p2 ^ 2/3 ^ J93 ^ 1, in 
agreement with the Kasner solution of GR. For n > 1/2 these solutions correspond 
to expanding universes with a curvature singularity at t = 0. For n = 1/2 the only 
solution is Minkowski space. 

For a universe filled with a perfect fluid having an equation of state p = (7 — l)p, 
7 constant, relating the fluid pressure p to its density p, the field equations ()2.1.12|1 
have the isotropic solution 

2n 

Pi=P2=P3 = — , (5.2.4) 
37 

for 7 7^ 0. This reduces to the spatially-flat FRW solution of GR in the limit n ^ 1. 
For n > 0, these isotropic cosmologies are expanding, and for n < they are con- 
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n 




Figure 5.1: For the solution i5. 2. S^) . the intervals in which the Kasner indices pi 
are allowed to lie can be read off' from this graph. For any value of n in the range 
^ ^ n ^ ^ a horizontal line is drawn; the boundaries of the intervals in which the 
Pi lie are then given by the four points at which the horizontal line crosses the two 
closed curves. For n = 1 these boundaries can be seen to be — |, 0, | and 1, as 
expected for the Kasner solution of general relativity. 

tracting, with n = giving Minkowski space. The stabihty and observational con- 
sequences of cosmologies of this type were investigated in [73 EH] , where primordial 
nucleosynthesis and the microwave background were used to impose observational 
constraints on the admissible values of n. 
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Substituting f = f{Y) = {RabR'^Y into the fourth-order field equations (|2. 1.121) 
along with the metric ansatz ()5.2.H) gives the two independent equations, 

yn-l j^(p2 ^ g _ ^pg ^ p2g ^ g2)_ 

2(3P2 + p3 ^ 3g _ gpg ^ 2Q'^)n + 2(4P2 + 4Q - 8Pg)) = (5.2.5) 

and 

yn-i (^^2AP - 2P2 - 2P3 - 30g + 10Pg)n + (8P2 - 32P + 24g)n2 

+ (p2 - 4P + 2p3 + gg - lOpg - p^g + 3g2)) = o, (5.2.6) 

where P and Q are defined as before. Equations ()5.2.5|1 and ()5.2.6|1 have five classes 
of solutions. 

The first of these classes is given by 

p = g = 0. 

This is only satisfied for Pi = P2 = Ps = 0, which is Minkowski space. It can be seen 
from ()5.2.5j) and ()5.2.6j) that as this solution corresponds to F = it only exists 
for n ^ (for n < the premultiplicative factor in ()5.2.5j) and ()5.2.6|1 causes the 
left-hand side of those equations to diverge). 
The second class of solutions is given by 

p = g = 1. 

This is just the Kasner solution of GR, for which the values of the constants Pi are 
constrained to lie within the ranges — l/3^pi ^0^p2^2/3^p3^1. Again, 
this solution corresponds to F = and so is only valid for n ^ 0. 
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The third class of solutions is given by 

^ _ 3(1 - 3n + 4n2) ± + lOn - - W + 4871^) 

" 2(1 -n) 

^ 3 

The only solution belonging to this class corresponds to an isotropic and spatially 
flat vacuum FRW cosmology. The values of the constants Pi are all equal to P/ 3 in 
this case and the solution is valid for all n 7^ 1. 
The fourth class of solutions is given by 

P={l-2nf (5.2.7) 
Q = l-Sn+ 16n^ - 8n^. 

The solutions belonging to this class are anisotropic cosmologies which reduce to the 
standard Kasner form, P — Q — 1, in the limit n — > 1. For this class of solutions, 
the values of the constants pi are in general constrained to lie within the ranges 

(1 - 2nf -2A ^ 3pi ^ (1 - 2nf - A 
(l-2nf-A ^ 3p2^(l-2nf + A 
{l-2nf + A ^ 3p3 ^ {1 - 2nf + 2A 

where A = a/(1 — 2n)(l — 6n + An^) and the pi have again been ordered so that 
Pi ^ P2 ^ Ps- These boundaries arc shown in figure 5.2 which can be read in the 
same way as figure 5.1, by taking a horizontal line of constant n and noting the 
four points at which that line intersects the curves. These intercepts determine the 
allowed intervals for the values of the pi. 

For real- valued pi, the value of n must lie either in the range ni ^ n ^ n2 or in 
the range 1/2 ^ n ^ 713, where ni, n2 and 77,3 are the roots of the cubic polynomial 
1 — 6n + An^ — and are chosen such that rii < n2 < n^. These generahsations of 
the Kasner metric always correspond to expanding cosmologies, independent of the 
value of n. 



66 




Figure 5.2: For the solution \5.2. 7| ) the intervals in which the Kasner indices pi are 
constrained to lie can be read off' this graph, as with figure 5.1. A horizontal line of 
constant n is drawn; the boundaries of the intervals in which the pi lie are then given 
by the four points at which the horizontal line crosses the curves. For n = 1 these 
boundaries can be seen to be — |, 0, | and 1, as expected for the Kasner solution of 
general relativity. 

The fifth class of solutions to ()5.2.5|) and ()5.2.6|) is given by 
P = An-l 

Q = -3 + 12n- ± 2(1 - 2n)^/2{l - 4n + 2n'^). 

This class describes complex- valued pi for all values of n (except for n = 1/4, for 
which this class of solutions reduces to the first class) and is therefore of limited 
interest. 

The isotropic solution for a universe filled with a fluid with equation of state 
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p = (7 — l)p is given by the choice 

An 

Pl=P2=P3 = —■ 

37 

This solution reduces to the spatially-flat FRW cosmology of GR in the limit n 
1/2 and corresponds to an expanding universe for n > and to a contracting universe 
for n < 0. 

/ r? ryabcd\n 

— yi^abcdi^ ) 

Substituting f = f{Z) = {RabcdR^^'^'^T and (IKTTI into the field equations ^ll.Vl\ 
gives the two independent equations 



+ (y + + 2Pg - 2P^Q + 3Q^ - 105 + -^)n = (5.2. 



and 



^n-i fP' p 13Q 7PQ 3P'Q 

^ l^lT"^~^~T + ^ + ^ 2- + ^"^^ 

+ (6P + _ I6Q - 2PQ + 8S)n + 8(Q - P)n^^ = (5.2.9) 
where P and Q are defined as before and we have also now defined 

i 

Equations ()5.2.8p and ()5.2.9|) have four different classes of solution. 
The first class of solutions is given by 

P = Q = S = 0. 

The only solution that belongs to this class is Minkowski space. As RabcdR"'^'^'^ = 
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for Minkowski space it can only be a solution for n ^ 0, due to the premultiplier in 
equations ^^TT^ and ^^2M- 

The second class of solutions is given by 

P = Q = S = 1. 

The only solution that belongs to this class corresponds to pi = p2 = and = 1. 
Making the coordinate transformations z = t sinh z and t = t cosh z |117j allows the 
line-element ()5.2.H1 to then be written in the form 

ds"^ = -dp + dx^ + dy"^ + dz"^ 

which is clearly Minkowski space again. This is only a solution for tt, ^ 0, as with 
the first class of solutions. 

The third class of solutions to ()5.2.8|) and ()5.2.9p is given by 



3(1 - 2n + ± V-1 + lOn - 16^2 + 16^^ 
~ 2(1 -n) ' 

p2 
p3 

which only has the isotropic and spatially flat vacuum FRW cosmology as a solution, 
where pi = P2 = Ps = P/^- This is a solution for all 7^ 1. 
The fourth and last class of solutions is given by 

P = 4n - 1 

Q = l{16n2 - 8n - 1 ± 4:^/2[n{l - 2n) + S{1 - n)]}. 
3 

This class of solutions corresponds to anisotropic metrics with Z = and is unusual 
in that it cannot be expressed in the form P = constant and Q = constant. This 
feature means that the standard picture of a plane intersecting an ellipsoid is no 
longer valid for this class. Solutions in this class do not appear to have any range of 
n for which the constants pi take real values, and so are of limited physical interest. 
The isotropic solution for a universe filled with a fluid with equation of state 



69 



p = (7 — l)p is given by 

An 

Pl=P2=P3 = ^• 

37 

As before, this reduces to the spatially flat FRW solution of GR in the limit n —* 1/2 
and corresponds to an expanding universe for n > and to a contracting universe 
for n < 0. 

5.2.3 Investigation of the Effects of Matter 

We have identified above two generalisations of the Kasner metric of GR which are 
solutions to the scale-invariant theories of gravity we are investigating. The first 
of these is a solution to the theory defined by £ = i?"; the second is a solution 
to the theory C = {RabR""^)^- These are the first exact anisotropic solutions to be 
found for higher-order gravity theories. In this section we will investigate some of 
the properties of these cosmologies. 

The behaviour of these solutions is of particular interest when considering the 
Bianchi type VIII or type IX 'mixmaster' cosmologies. The field equations for 
these cosmologies can be cast into the form of the equations of motion of a particle 
moving inside an exponentially steep triangular potential well with open channels 
in the corners jl26j . The three steep-sided walls are created by the 3-curvature 
anisotropics. In the region where the potential is negligible (far from the walls) 
the behaviour of the solutions approaches that of the Kasner metric. As the expo- 
nentially steep potential wall is approached the universe 'particle' is reflected and 
re-enters a Kasner-like regime with the Kasner indices pi systematically changed to 
some new values by the rule of reflection from the potential wall. In general rela- 
tivistic cosmologies of Bianchi types VIII and IX this process is repeated an infinite 
number of times as the singularity is approached [IHl |Sl El EH H'^^j so long as 
matter obeys an equation of state with p < p. After reflection from the potential 
wall the Kasner index that was previously negative is permuted to a positive value 
and the lowest-valued positive Kasner index is permuted to a negative value. This 
is repeated ad infinitum in the general relativistic cosmology as one of the Kasner 
indices must be kept negative while the other two are positive. In the generalisations 
of the Kasner metric presented above this is no longer the case; for some solutions 
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it is possible for all of the Kasner indices to be positive. If the Kasner indices are 
being permuted in a chaotic fashion for long enough then, eventually, they will end 
up in such a configuration. Once this occurs the oscillations will end, as all spatial 
dimensions will be contracting, and the singularity will be reached monotonically 
without further oscillation of the scale factors lllj. 

From figure 5.1 it can be seen that for solution ()5.2.H|1 all Kasner indices can be 
made positive for values ofl<n<5/4in the L = i?" theories. Similarly, from 
figure 5.2 it can be seen that for the solution ()5.2.7|) . with rii < n < n2 ot 1 < n < n^, 
all the Kasner indices can be made positive in the theory C = [RabR"'^)"' (where the 
rii are defined as before). For both of these solutions, with these ranges of n, there 
is therefore no infinite series of chaotic mixmaster oscillations as the singularity is 
approached. However, the solutions to both of these theories are still expected to 
exhibit an infinite number of chaotic oscillations if 1/2 < n < 1, as at least one of the 
Kasner indices must be kept negative whilst another is kept positive. A separate 
detailed analysis is require to determine where there are differences between this 
behaviour and the chaotic oscillations found in GR. 

It now remains to show that the vacuum solutions found above are the asymp- 
totic attractor solutions in the presence of non-comoving matter motions as t ^ 0. 
This analysis will follow closely that of |2S1 and |117j . We aim to show that the 
fluid stresses diverge more slowly than the vacuum terms as t — and so produce 
negligible metric perturbations to an anisotropic Kasner like universe. If this is the 
case, then the vacuum solutions above can indeed be considered as the asymptotic 
attractor solutions even in the presence of matter. Matter will just be carried along 
by the expansion and behave like a test fluid. 

We now consider a perfect fluid with equation of state p = (7 — l)p where 
1 ^ 7 < 2 and the non-zero 4- velocity components, W, are normalised so that 
WUi = —1. The conservation equations T"*;fc = on the metric background ()5. 2.111 
can then be written in the form |117| ITU IHj 




d 



ip, - (7 - mu^p^k. 
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Neglecting spatial derivatives with respect to time derivatives, these equations inte- 
grate to 

t^Uopi = constant 

(7-1) 

Uap ^ = constant. 

The neglecting of spatial derivatives means that these equations are valid on scales 
larger than the particle horizon as in the velocity-dominated approximation (al- 
though we have not so far restricted the fluid motions to be non-relativistic) |117t 

EaiiHii- 

From the second of the integrals above, it can be seen that all the covariant 
components of the spatial 3- velocity, Ua, are approximately equal. This is not true 
of the contravariant components as the Kasner indices, Pi, in the metric elements 
used to raise indices are not equal in these solutions. The contravariant component 
that diverges the fastest, and dominates the others in the t —>■ velocity-dominated 
limit, is therefore = U^t'^^^, as ps is the largest of the pi. If the 4-velocity is 
normalised, so that UaU"" = —1, and the contravariant 3-velocity component 
diverges the fastest as t — 0, then we must have UqU'^ ~ U^U^ = {U^Yt'"^^^ in that 
limit. The integrated conservation equations above can now be solved approximately 
in this limit to give 

p ^-7(Pl+P2)/(2-7) 

^ ^(pi+P2)(7-l)/(2-7) 

as t ^ 0. 

It is now possible to calculate the leading-order contributions to the energy- 
momentum tensor, T"fe = (p + p)U°'Uh + 

Tl p ~ i-7(P-P3)/{2-7) 

for 7 < 4/3. The component which diverges the fastest here as t —> is T| ~ t"^"^^. 
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for general 7. For the case 7 = | where the 4- velocity of the fluid is comoving, 
Ui = 6f, all the components of T^^ diverge as p ~ ^-2P+2p3_ 

For the generalisation of the Kasner metric in the C = i?" theory, ()5.2.3|) . the 
vacuum terms diverge as t"^". We therefore require that 2n > P + ^3, or p3 < 1, in 
order for the vacuum terms to dominate over the matter terms for a non-comoving 
perfect fluid, with general 7, in the limit t 0. For the 7 = |, comoving perfect fluid 
(f/i = 5°), the condit ion for vacuum domination is 2n > 2P — 2p^, or > n — 1. 
Both of these conditions are ensured by the boundary values on given above. 
In the velocity-dominated limit, the vacuum solutions given here are therefore the 
appropriate asymptotic solutions on approach to the singularity. 

For the generalised Kasner metric in the C = {RabR""^)^ theory, ()5.2.7|1 . the 
vacuum terms diverge as For vacuum domination as t in this solution 

we therefore require An > P + p^, or p^ < — 1 + 8n — 4n^, for the general non- 
comoving 7 fluid; and 4n > 2P — 2p3, or > 1 — 6n + 4n^, for the 7 = 4/3 
comoving fluid. From figure 5.2 it can be seen that solutions of this class can 
be in one of two regions, rii ^ n ^ n2 or 1/2 ^ n ^ n-^. The validity of the 
vacuum solution as t — > is different in these two different regions. For the region 
1/2 ^ n ^ the boundary conditions on the index p^ mean that the conditions 
above are automatically satisfied. For the region ni ^ n ^ n2 these inequalities are 
not always satisfied. For general 7, the condition for vacuum domination is only 
met if n lies in the narrow range 1/6 < n < 722; for any other value of n in this 
region the index p^ can be such that the fluid diverges faster than the vacuum terms. 
Similarly, for the 7 = 4/3 comoving fluid the condition for vacuum domination is 
only satisfied if n lies in the narrow range < n < n2, where is the real root of 
the cubic 

4n^ - 12^2 + lOn - 1 = 0. 

To understand the evolution of an anisotropic solution of the form 1)5.2.11] in the 

theory C = {RabR"''^)^, where rii ^ n ^ 77-2, it is therefore necessary to take into 
account the relativistic motions of any fluid that is present (except in the narrow 
ranges of n identified above). This range of n can, however, be regarded as not 
belonging to physically interesting theories on other grounds. In order to agree with 
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weak-field experiments it is necessary for a gravitational theory to contain a term 
which approximates the Einstein-Hilbert action in the weak-field limit (see ref. |7H]). 
Any theory containing such a limit must therefore have a term in its Lagrangian 
that diverges as t"^. In considering the behaviour of alternative theories close to 
the singularity it is necessary for the extra terms in the Lagrangian to diverge faster 
than this if they are to be influential as t — 0. For the theory C = {RabR"'^)"' this 
requires n > 1/2. For n < 1/2 such a term will diverge slower than the Einstein- 
Hilbert term which will then dominate and display the standard Kasner behaviour 
of GR. For the anisotropic solution ()5.2.7|1 we are therefore only interested in the 
region lying in the range 1 /2 ^ n ^ na , which has vacuum terms that dominate the 
fluid stresses in the velocity-dominated approximation. 

It remains to investigate the effects of fluids with equations of state, 7 > 4/3. 
For sufficiently stiff fluids the velocity-dominated approximation is not valid and 
U°'Ua — * as t ^ due to the very high inertia of the fluid producing a slow down 
under contraction jEJl- In such a limit Uq ^ 1 and the conservation equations can 
be solved to give 

In this approximation we can write UaU" ~ (t^a)^^"^^^ ~ ^2(7-i)p-2p3_ ^^^^ ^^^^ 
be seen that this behaviour occurs when 

P3 ^(7-1)^- (5.2.10) 

In this limit it is only required that p diverges more slowly than the vacuum terms, 
as pUaU^ « P when t ^ 0. 

For the solution ()5.2.3|1 to the theory C = i?" it can be seen from the condition 
fl5.2.1()|l that it is necessary for 7 ^ 4/3 in order for the velocity-dominated approxi- 
mation to break down (this is derived using the upper limit on n and the lower limit 
on P3). The condition that the fluid stresses diverge more slowly than the vacuum 
terms is now 

2n 

^ < (*r3T) (".11) 
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where it has been assumed n > 1/2. A sufficient condition to satisfy this for all 
allowed values of n is 7 < 5/3. For fluids with a stiffer equation of state, 5/3 ^ 7 < 2, 
the vacuum terms dominate over the fluid motions provided that n < 1. For n < 1, 
however, the Einstein- Hilbert term will be the leading one in the gravitational action, 
as discussed above. 

Similarly, for the solution ()5.2.7j) to the theory C = {RabR"'^)"' it can be seen 
from the condition ()5.2.10|) that it is again necessary to have 7 ^ 4/3 in order for 
the velocity-dominated approximation to break down. In the region 1/2 ^ n ^ 
the condition for the vacuum solution to be unperturbed by the stiff fluid is always 
satisfied for any fluid with equation of state 7 < 2. However, in the region rii < 
n < n2 the condition for the vacuum to diverge faster than the energy density of 
the stiff fluid is never satisfied for a fluid which satisfies the necessary condition for 
the velocity-dominated approximation to break down (except in the narrow range 
(7— -\/33)/8 < n < n2 where there are some values of 7 for which the vacuum terms 
can diverge fastest). 

5.2.4 Discussion 

We have investigated some anisotropic cosmological solutions to higher-order La- 
grangian theories of gravity. Whist the standard general relativistic theory appears 
to be consistent with all weak-field tests, there is less reason to think that it should 
be valid in high curvature regimes, such as in the vicinity of a possible initial cosmo- 
logical curvature singularity. In fact, it is in the high-curvature limit that quantum 
effects should become important and we should expect to see deviations from the 
standard theory. Without knowing the exact form of such deviations, we have ap- 
proached this problem by considering a general class of theories that can be derived 
from an arbitrary analytic function of the three curvature invariants R, RabR"^^ and 
RabcdR"'^'^'^ ■ Expanding this function as a power series in these variables we then 
expect the dominant term in the Lagrangian to be of the form i?", {RabR"'^)^ or 
{.RabcdR°''^'^'^Y the singularity is approached. We have found all of the solutions 
to these theories that can be expressed in terms of the Bianchi type I line-element 
()5.2.1|1 in vacuum. These solutions provide simple testing grounds for the explo- 
ration of quantum cosmological effects in higher-order gravity theory. We have also 
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found the homogeneous, isotropic and spatially flat solutions to these theories in the 
presence of a perfect fluid. 

Exact vacuum anisotropic solutions of the form ()5.2.H) were found for the theories 
£ = -R" and £ = {RabR"'^)"'-, whilst it was shown that no such solutions exist 
for theories defined by £ = {Rabcd.R"''"^'^)"' ■ The properties of these solutions, with 
respect to their relation to the more general Bianchi type VIII and IX cosmological 
behaviour, has been investigated. We have argued that for all of the physically 
relevant new solutions of these theories, the Universe will not experience an infinite 
number of mixmaster oscillations as the singularity is approached. The extent to 
which these vacuum solutions can be considered as realistic in the presence of a 
perfect fluid has also been investigated. It has been shown in the velocity-dominated 
approximation that all the anisotropic vacuum solutions found for plausible theories 
are valid in the limit t ^ 0. The case of stiff fluids that do not satisfy the velocity- 
dominated approximation as t ^ has also been investigated. For the solutions to 
the theories C = R^ it was found that for a fluid with equation of state 7 < 5/3 
the vacuum solutions are good approximations in the vicinity of the singularity. 
For the theories C = {RahR^^T it 

was found that the vacuum solutions are good 

approximations for all 7 < 2. 
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CHAPTER 6 



Static Solutions 



In this chapter the vacuum static and spherically symmetric solutions to the Brans- 
Dicke and i?" theories will be presented. For these theories there is no equivalent 
of 'Birkhoff's theorem' (in fact, in the next chapter we will present solutions which 
explicitly violate Birkhoff's theorem). The condition of being static is therefore 
added as an extra condition, in contrast to the usual treatment in GR. 
This chapter is based on the work of Clifton and Barrow fTOj . 

6.1 Scalar- Tensor Theories 

The static and spherically symmetric solutions of the Brans-Dicke theory were found 
long ago by Brans himself [SI]. These solutions have been investigated and redis- 
covered in many papers since, and so we will only present a very brief summary of 
the situation here, for completeness of this study. In isotropic coordinates 
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where a = a{r) and (3 = f3{r), there are four solutions of the Brans-Dicke field 
equations ()2.1.8|1 and ()2.1.9j) which are given by 
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and B, ao, (3o and 0o are constants. Whilst solution / is valid for all values of uj, 
solutions //, III and IV are only valid for u < —3/2. It is probably for this reason 
that it is most often solution I that is used in the literature for analyses of the 
static and spherically symmetric situation. Solution I is conformally related to the 
minimally coupled massless scalar field solution of Buchdahl [HB] • 

It can be seen that these solutions are not all independent of each other. By a 
transformation of the form 

1 

r ' 

and some redefinition of constants, solution // can be transformed into the uj < —3/2 
range of solution / jlH] and solution III can be transformed into solution IV |3H]- 
It was further shown in jlH] that the independent solutions / and IV are both 
conformally related to the general solution of the static and spherically symmetric 
case in the Einstein conformal frame, given by Wyman |192j . 



6.2 Fourth-Order Theories 

In the absence of any matter the C = R^~^^ field equations ()2.1.15|) can be written 
as 

f R."^ R'R'^\ ( 1(1 + 25) \ 

^"'^ " I ~R~ ~ ~ ^) '^2 ' j U'^c^fcrf + 2 ^"''^"^ ) ■ (6-2.1) 

We find that an exact static spherically symmetric solution of these field equations 
is given in Schwarzschild coordinates by the line-element 

ds^ = -A(r)dt^ + ^ + r^de^ + sin^ ed(j)^) (6.2.2) 
B[r) 



where 



r,c{l + 25) C 

A{r) = r''^ + 

r (1-*) 



{1-26 + 452)(i - 25(1 + 6)) \ 
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and C = constant. This solution is conformally related to the Q = limit of the 
one found by Chan, Horne and Mann for a static spherically-symmetric space-time 
containing a scalar- field in a Liouville potential IMJ. It reduces to Schwarzschild in 
the hmit of GR: 6^0. 

In order to evaluate whether or not this solution is physically relevant we will 
proceed as follows. A dynamical systems approach will be used to establish the 
asymptotic attractor solutions of the field equations 1)6.2.11) . The field equations will 
then be perturbed around these asymptotic attractor solutions and solved to first 
order in the perturbed quantities. This linearised solution will then be treated as 
the physically relevant static and spherically-symmetric weak-field limit of the field 
equations ()6.2.1|1 and compared with the exact solution ()6.2.2|1 . 

6.2.1 Dynamical System 

The dynamical systems approach has already been applied to a situation of this 
kind by Mignemi and Wiltshire |125j . We present a brief summary of the relevant 
points of their work which is relevent to this study, in the above notation. 

Taking the value of sign(i?) from (jHT^ as sign(-5(l + 6) /{I - 25(1 + 6))) and 
making a suitable choice of flo allows the scalar-field potential ()2.2.14|1 to be written 



In four dimensions Mignemi and Wiltshire's choice of line-element corresponds 

to 



which, after some manipulation, gives the C = R^^^ Einstein frame field equations 



as 




(6.2.3) 



(6.2.4) 



as 



7] 



II 



II 



2c\{\ - bf + 6,5277^2 _ 24,52r7^C - 2(1 - 2(5 - 8,52)^2 

1 - 25 + 452 ^ 
(1 - 25 - 852)(c2(l - 5)2 + 352r/^2 _ YlS'rjC: - (1 - 25 - 8^2)^2) 

352(1-25 + 452) 



(6.2.5) 




(6.2.6) 
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and 

-{1 - 26- 86^)C^ + {1-26 + 452)e2«) (6.2.7) 

where 

C(0 = f^(0+logr(0 

^(0 = -^-^—^^^U{0 + 2logf{0 - ^i^cie + constant. 

Primes denote differentiation with respect to ^ and ci is a constant of integration. 
The variables X, Y and Z are here defined to be 

X = C Y = 7]' Z = e^ 

(not to be confused with the definitions of the scalar curvature invariants used in 
other chapters). The field equations ()6.2.5|) and ()6.2.6p can then be written in terms 
of these new variables as the following set of first-order autonomous differential 
equations 

yj_ 2c\{l-6f + ^6^Y^- 2^6^ XY -2{1- 26- U^)X^ ^ 8) 



Y 



1-26 + 

, {1-26- 852)(c2(i - 6f + 352^2 _ _ (1 _ 25 - 852)x2) 



3(52(1-25 + 4^2) 



(1-25(1 + 5)) , 
+ ^ (6.2.9) 

Z' = XZ. (6.2.10) 

(The reader should note the different definition of Z here to that of Mignemi and 
Wiltshire). As identified by Mignemi and Wiltshire, the only critical points at finite 
values of X, Y and Z are in the plane Z = along the curve defined by 

cl{l - 6^ + 35^2 - 125^X1" - (1 - 25 - 86'^)X^ = 0. 

These curves are shown as bold lines in figure 16.11 together with some sample tra- 
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Figure 6.1: The Z = plane of the phase space defined by X,Y and Z for 5 = 0.1 and 
Ci = 0.5. The bold lines show the critical points in this plane and the diagonal lines 
show the unphysical trajectories confined to this plane. The dotted line isY = 2X 
and separates the critical points where ^ ^ oo from the points where ^ —>■ —oo 

jectories from equations ()fi.2.8j) and ()fi.2.9|l . From the definition above we see that 
the condition Z = is equivalent to fe^ = 0. Whilst we do not consider trajecto- 
ries confined to this plane to be physically relevant, we do consider the plot to be 
instructive as it gives a picture of the behaviour of trajectories close to this surface 
and displays the attractive or repulsive behaviour of the critical points, which can 
be the end points for trajectories which could be considered as physically meaning- 
ful. The dotted line in figure 16.11 corresponds to the line Y = 2X and separates 
two different types of critical points. The critical points with Y > 2X can be seen 
to be repulsive to the trajectories in the Z = plane and correspond to the limit 
^ — > — oo. Conversely, the points with Y < 2X are attractive and correspond to 
the limit ^ — oo. As Z = fe^ , all critical points of this type in the Z = plane 
correspond either to naked singularities, f — 0, or regular horizons, f — >^constant. 
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The two bold lines in figure Ifi.ll are the points at which the surface defined by 



cl{l - 5f + 3(5V^ - m'^XY -{1-25- %5'^)X'^ + {1-25 + A5'^)Z'^ = 



crosses the Z = plane. This surface splits the phase space into three separate 
regions between which trajectories cannot move. These regions are labelled J, // 
and /// in figure 16.11 It can be seen from ()6.2.7|) that trajectories are confined 
to either regions / or II for the potential defined by ()6.2.3p . If we had chosen 
the opposite value of sign(i?) in ()2.2.14|1 then trajectories would be confined to 
region ///. We will show, however, that region III does not contain solutions with 
asymptotic regions in which f — > cxd and so is of limited interest for our purposes. 

In order to find the remaining critical points it is necessary to analyse the sphere 
at infinity. This can be done by making the transformation 



and taking the limit p ^ oo. The set of equations ()fi.2.8j) . ()fi.2.9|) and ()(12.1()j) then 
give 



X = p sin 6 cos 



Y = p sin 6 sin 



Z = pcosO 



de 

dr 



24(52(1-25 + 45) 
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{65^ cos 0(3 - 3(5(2 -95) + {1 -5{2 + 115)) cos 29) 
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-2(18(5^(1 -5(2 + 7(5)) cos 30 



- (1 - 5(4 + 5(9 - 265 + 325^))) sin 30) sin^ 6) 



and 



d0 
d^ 
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(65^(1 - 5(2 + 415) - 5(1 - 5(2 + 55)) cos 2^)cosec^ sin 



+ 2((1 - 5(4 + 5(9 - 265 + 325^))) cos 30 

+ 185^(1 -5(2 + 75)) sin 30) sin^ 
2 cos 0(4(1 - 25(2 - 5(3 - 25 - 452)))cosec^ 

- (7 - 5(28 + 5(15 - 25(43 + 1285)))) sin^)) 
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0.2 - 



Figure 6.2: The surface at infinity of the phase space defined by X,Y and Z for 
6 = 0.1. The shaded areas show where regions I and II . Region III is unshaded. 

where dr = pdS,. These equations can be used to plot the positions of critical points 
and trajectories on the sphere at infinity. The result of this is shown in figure 16.21 
Once again, these trajectories do not correspond to physical solutions in the phase 
space but are illustrative of trajectories at large distances and help to show the 
attractive or repulsive nature of the critical points. The surface at infinity has eight 
critical points, labelled A-H in figure 16.21 Points A and B are the end-points of the 
trajectory that goes through the origin in figure lOl and are located at 

, , / -6(5^ \ , , 

6 = and cpi lo-) = cos , (+7r) 

2' "^ '^ ^ VVl - 4(5 - 1252 + 3253 + 10054 J ^ ' 

or, in terms of the original functions in the metric ()6.2.4|) . 

33^ 33^ 
f ^ - (1-23-83^) and 

where ^1 is a constant of integration. The points A and B therefore both correspond 
to ^ ^ ^1 and hence to f — * 0. 

Points C, E and F are the four end points of the two curves in figure EH] and 
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therefore correspond to ^ — oo or — oo and r — or constant. 
The remaining points, G and H, are located at 



^1,(2) 



TT, , ^ ^ I f 1-2(5+10(52 

— l+TTl and = - cos 



4 ' ' 2 V 3 - 6(5 + 6(52 



or 



and 




1 - 2(5 - 2(52) 1 
-25 + 452) (^-^2 



1 - 2(5 - 252) 1 



;i-25 + 452)(e-6)' ^^-^-^^^ 

where ^2 is another integration constant. The positive branch corresponds to point 
H and the negative branch to point G. These points are, therefore, the asymptotic 
hmit of the exact solution ()6.2.2|1 and correspond to ^ ^ ^2 and hence f — > 00. 

Whilst it may initially appear that trajectories are repelled from the point H, 
this is only the case in terms of the coordinate ^. In terms of the more physically 
relevant quantity f, the point H is an attractor. This can be seen from the first 
equation in ()6. 2.111) . Taking the positive branch here it can be seen that f increases 
as ^ decreases. So, in terms of f the points G and H are both attractors, as f — > 00. 

We can now see that in region / all trajectories appear to start at critical points 
corresponding to either f — > or constant and end at point G where f — > 00. Region 
II appears to share the same features as region / with all trajectories starting at 
either f — or constant and ending at H where f — 00. Region III has no critical 
points corresponding to f ^ 00 and so all trajectories both begin and end on points 
corresponding to f — or constant. 

Therefore solutions with an asymptotic region in which f ^ 00 only exist in 
regions / and // where the potential can be described by equation ()6.2.3j) . Further- 
more, all trajectories in regions / and // are attracted to the solution 



6<S- 



ds" = -fi^de+ ^^ ^^^^^'^^L '^^^^ df^ + f\de'' + sm^ed(t)^), (6.2.12) 

(1-5)4 

which is the asymptotic behaviour of the solution found by Chan, Horne and Mann 
We therefore conclude that all solutions with an asymptotic region in which 
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f — > oo are attracted towards the solution ()fi.2.12|l as f — ^ oo. 

Rescaling the metric back to the original conformal frame, we therefore conclude 
that the generic asymptotic attractor solution to the field equations, ()6.2.H) . is 

as- = '1 - + - - ^'"^ dr- + r\de- + sin-ed^-) (6.2.13) 

(1 - d)' 

as r — > oo, which reduces to Minkowski space in the 5 — > limit of GR. 
6.2.2 Linearised Solution 

We now proceed to find the general solution, to first order in perturbations, around 
the background described by ()(i.2.1Hj) . Writing the perturbed line-element as 

o r,r (1 + 26) „ 

ds^ = -r^^lT^(l + V{r))dt^ 

+ + + w{r))dr^ + r\de' + sin^ ed<f)') (6.2.14) 

{i-6y 

and making no assumptions about the order of R, the field equations ()6.2.H) become, 
up to first order in V and W, 

5{l + 26){l + 26^) (1 + 26^) V 5{l - 25) W 

(1 - 5)2r2 ^ (1-5) 7" ~ 2(1-5) V ^ T 

_ 5{l + 25) R'^ _ 5{l + 25) R" _ 35 R' 
~ 2 2(1-5) "R 4(1-5)^ 

5{l + 25){2 + 5) R' 5{1 + 25) R' 
- 2(1 -5)V R+W^R ' ^ ^ 



5(l + 25)(l-25-252) 5(1 + 25) r (2-5 + 25^)1^' 



(1 -5)2r2 



35 i? 



/2 



+ 



;i-5) r 
35 i?" 



2 2(1 - 5) i? 



+ 



2(1-5) r 
5(1 + 25) (2 -5 + 252) i?' 



+ 



5(1 



2(1 -5)2r 
25) i?' 



R 

35 R^ 

4(1-5) R' 4(1-5)^ 



ly' (6.2.16) 
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and 



_ 2^(3 - 4(5 + 2(5' + 86"^) ^ 2(1 - 2(5 + 4(5')(1 - 2^ - 2(5') _Y^^^ 
(1 — Syr"^ (1 — 6yr'^ r r 
o.N^" (5(1 + 2(5) i?" (5(4 - (5 + 2(5' + 4^3) 
= -5{l + 26)— + -j^^ + ^Y^r R 

6(1 + 26) R 6(1 + 26)^ (6.2.17) 
2(1-6) R 2(1-6) R ^ ^ 

Expanding i? to first order in V and W gives 

^--(1 - 26- 26^)7^^^' ^^-^-^^^ 



where 

2(l + 5 + (5') ly 2(l-(5)(l + 2(5') r 



(1 - 2(5 - 25') r' (1 - 2(5 - 2(52)(1 - 25 + 45') r 

(1 -5) (2 -5 + 2(5') ly' 
(1 - 25 - 2(5')(1 -26 + 452) V 

(1-5)' 



+ 



(1 -25- 25') (1 -25 + 45') 



V". (6.2.19) 



Substituting (l(i.2.18|) into the field equations (16.2. 15|) . (16.2. 16jl and (I6.2.17|l and 
ehminating Ri using ()6.2.19j) leaves 

(1 + 5 + 5')(5 - 125 + 125' + 45^) W (16 - 475 + 765' - 345^ - 165^ + 325^) , 
3(l-5)'(l + 5) 6(1 -5')(1 - 25 + 45') 

(1 + 5 + 75' - 1953 + 445^ + 205^) Y (8 - 155 + 185' + 165^) 



3(l-5')(l- 25 + 45') r 6(1 + 5)(1 - 25 + 45') 

(1 - 25 - 25') (5 - 125 + 125' + 45^) ^ (1 - 25 - 25') 



12(l-5')(l + 5) r 4(1-5') 
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(1 + 26){l + 6 + 6^){3 - 45 + 45^) w _ {1 + 26){2 - 6 + 26^){3 - A6 + 45^) 
3(1-5)2(1 + 5) ~ 6(l-52)(i_25 + 452) 

(3 - 25 + 1752 - 453 + 4052) Y (6-5 + 252 + 205^) 
^ 3(l-52)(l-25 + 452) 7^ 6(l + 5)(l- 25 + 452) 
_ (1 + 25)(1 - 25 - 252)(3 - 45 + 452) y; (1 + 25)(1 - 25 - 252) , 
~ 12(1-5)2(1 + 5) 7 ^ 12(1-52) ■ 



and 



2(8 - 85 + 352 + 105=^ - 285^ - 125^) ly (13 - 225 + 1252 + 265^-565^) 
3(l-52)(l + 5)(l- 25 - 252) 7^3(1 - 52)(1 - 25 - 252)(1 - 25 + 452) 
2(4 + 952 + 853-125^) Y (5-45-452 + 1253) 

^ 3(1 + 5)(1 - 25 - 252)(1 - 25 + 452) 7 ^ 3(1 + 5)(1 - 25 - 252)(1 - 25 + 452) 

_ (5 - 45 - 452 + 1253) ip (1 + 25) , 
~ 6(1-5)2(1 + 5) 7 ^ 6(1-52)"^ 

where y = rV and = r^i?'^, subject to the constraint ()6.2.19|1 . 

For —-^^-^1^-^ < 5 < — -^^-f^^ the general solution to this first order set of 
coupled equations is given, in terms of V and W, by 

V{r) = CiVi{r) + C2V2ir) + c^Vsir) + constant (6.2.20) 
W{r) = -ciVi{r) + C2W2{r) + csW^ir) (6.2.21) 

where 

^2 = _ 3^ ; ^ ^g^3^ ((1 + 25)^ sin(Alogr) + 2A{1 - 5) cos(Alogr)) 

(1-26+4,6^) 

W2 = r sin(Alogr) 

(l-2i5+4f^) 

n 4- 25V 2(i-«) 

V^3 = 2(2 - 35 + 1252 + 1653) + 2^)' cos(Alogr) - 2A{1 - 5) sin(Alogr)) 
= r 2(1-5) COS (A log r) 
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and 

V7 - 28(5 + 36^2 - 16(53 - 80F 

which is a real number for the any 5 in the range specified above. The extra con- 
stant in ()6.2.20|) is from the integration of Y and can be trivially absorbed into the 
definition of the time coordinate. The above solution satisfies the constraint ()6.2.19|) 
without imposing any conditions upon the arbitrary constants Ci, C2 and C3. 

It can be seen by direct comparison that the constant ci is linearly related to 
the constant C in flfi.2.2j) by a factor that is a function of 5 only. The constants C2 
and C3 correspond to two new oscillating modes. 
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CHAPTER 7 



Inhomogeneous Cosmological Solutions 



We will present in this chapter inhomogeneous cosmological solutions to the scalar- 
tensor and fourth-order field equations. Considerations will be limited to space-times 
admitting spherically symmetric three dimensional space-like hypersurfaces. The 
existence of vacuum solutions with this symmetry shows explicitly that Birkhoff's 
theorem cannot be formulated in these theories. 

In the scalar-tensor case we will also present a study of how the spherical collapse 
model works in these theories. Closed universes are matched to flat universes at a 
boundary and allowed to evolve seperately. Following the evolution of the scalar 
fleld in each region then allows for consideration of how spatial inhomogeneity in 
Newton's gravitational constant G can arise in a realistic universe. 

This chapter is based on the work by Clifton, Mota and Barrow [ZSj and Clifton 

7.1 Scalar- Tensor Theories 

Extensive studies have been made of cosmological solutions of scalar-tensor gravity 
theories [HHl, although they are hmited in two respects. First, they usually focus on 
the simplest case of isotropic expansion with zero spatial curvature, where simple 
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exact solutions exist. Second, they are almost exclusively concerned with spatially 
homogeneous cosmologies. The latter restriction means that the value of G and its 
rate of change in time, G, are required to be the same everywhere in the universe. 
We want to understand how G and G are expected to vary in space in a realistic 
inhomogeneous universe. Since, even on the scale of a typical galaxy, the amplitude 
of visible density inhomogeneities are of order 10^, we need to go beyond linear 
perturbation theory for such an analysis. 

Previous studies on inhomogenieties in physical constants have focussed on the 
formation of black holes (gravitational memory) and perturbative expansions (see 
[IHl 122 EHl inni IIZ3 inni inZI). in this work we confront this problem by tracing 
the evolution of G, first in exact inhomogeneous solutions and then in a simple, 
but not unrealistic, inhomogeneous universe in which a zero-curvature Brans-Dicke 
FRW background universe is populated by spherical overdensities which are mod- 
elled by positive curvature Brans-Dicke FRW universes in the dust-dominated era of 
the Universe's history. This will enable us to track the different evolution followed 
by G{t) in the background universe and in the overdense regions, which eventually 
separate off from the background universe and start to contract to high-density sep- 
arate closed universes. This process can produce significant differences between G 
and G in the background and in the overdensities. Eventually, the collapse of the 
spherical overdensities will be stopped by pressure and a complicated sequence of 
dissipative and relaxation processes will lead to virialisation. This state will provide 
the gravitational environment out of which which stars and planetary systems like 
our own will form, directly reflecting the local value of G{t) inherited from their 
virialised protogalaxy or its parent protocluster. The simple model we use for inho- 
mogeneities in density and in G has many obvious limitations, notably in its neglect 
of pressure, deviations from spherical symmetry, accretion, and interactions between 
inhomogeneities. Nonetheless, we expect that it will be indicative of the importance 
of taking spatial inhomogeneity into account. It provides the first step in a clear 
path towards improved realism in the modelling of inhomogeneities that mirrors the 
route followed in standard cosmological studies of galaxy formation with constant 
G. 
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7.1.1 Exact Solutions 



Starting with a solution of Einstein's field equations with a scalar field we can apply 
a conformal transformation to arrive at a solution of the Brans-Dicke field equations 
in a vacuum. A spherically symmetric exact solution for the collapse of a minimally- 
coupled scalar field, ip, in GR is known and is given by |l()5j 



ds' = -{qt + b){fir)dt' - f-'{r)dr') + R'{r,t){de' + sin' edr), 



(7.1.1) 



where 



fW = (i--)^ 

r 

R\r,t) = (gt + 6)r2(l- —)!-", 

r 



a = ± 



The evolution of the minimally coupled scalar field, tp, in the Einstein frame, is 
given by 



i^{r, t) 



1 



In 



2c\ V3 



r 



d[l (qt + b) 



V3 



(7.1.2) 



Here, q, b, c and d are constants. Now under the transformation to the Jordan 



frame, ()2.2.6|) . where = exp 



we obtain 



ds 



7,2 



[A{rYde - A(r)-°rfr2] 



A(r)'""^^S(t)i-^/'3r^ ,,2 2 



and 



where A{r) 



t) = [dA(r)"/^5(t)^ 



(7.1.3) 



(7.1.4) 



2c 



B{t) = qt + b and [3 = ±y/2u + 3. We now assume that g 7^ 



3 V3 



[i.e. the metric is not static) and define the new time coordinate t = (qt + b)'^ 2/3 
In the limit that c — > the r-dependence of the metric is removed and the space 
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becomes homogeneous. In this case we expect ()7.1.Hj) to reduce to the FRW Brans- 
Dicke metric given in the last section. We see from the form of ()7.1.3j) that the 
metric should reduce to that of a flat FRW Brans-Dicke universe. Insisting on this 
limit requires us to set j3 = y/2uj + 3, q = 3^^^^ and d = 1. This leaves the metric: 

+ A{ry"^^^h^^ [dr^ + A{ry{de^ + sin^ Odcf)^)] . (7.1.5) 
Rewriting ()7.1.4j) with these coordinates and constants gives 

0(r,t)= (1-7) ''?/(v^/5-i). (7.1.6) 

A comparison of ()7.1.5|1 with the FRW solution ()3.2.12jl shows that ()7.1.5j) does 
indeed reduce to a flat vacuum FRW metric in the limit c —>■ (an inhomogeneous 
universe requires c 7^ 0). The metric ()7.1.5|1 is asymptotically flat and has singu- 
larities at t = and r = 2c; the coordinates r and t therefore cover the ranges 
< t < 00 and 2 < - < 00. 

— — c 

The equations ()7.1.39j) and ()7.1.6j) can now be used to construct a plot of G{r, t); 
this is done in figure mi which was constructed by choosing the minus sign in ()7.1.6j) . 
From the form of G{r, t) we see that this choice corresponds to an overdensity in the 
mass distribution (identified by comparison with the inhomogeneous Brans-Dicke 
solution with matter, found below). In this figure u was set equal to 100 and c to 
1. This plot shows how G can vary in space and time in an inhomogeneous universe 
which consists of a static Schwarzschild-like mass sitting at r = in an expanding 
universe. As r ^ 00 it can be seen that this solution approaches the behaviour of a 
Brans-Dicke FRW universe. 

An Inhomogeneous Brans-Dicke Solution With Matter 

We now seek a solution of the Brans-Dicke field equations ()2.1.8j) with the form 

ds^ = -e^de + e^a^{dr^ + r^dd^ + sin^ dd<p^) (7.1.7) 
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G (r, t) 
G (2.1, 5) 



100 



Figure 7.1: This graph illustrates the possible space and time variations that can arise 
in G in inhomogeneous solutions to the Brans-Dicke field equations with uj = 100, 
normalised at r = 2.1 and t = 5. 

where e" = e^(^), = e''^^) and a = a{t). Substituting (ITTTjl into gives, for 

a perfect fluid 



Stt 

—P = - 



a \a J 2 \ (p I acj) (p 



+ 



, u' 2 



e-^a-^ (7.1.8) 
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+ 



\ h h - 

2 4 2r 2 



+ — + 



u' 1 
2" ^ r / 



(7.1.9) 



and 



»7r 



3(^ 

a 



2\(j)] act) 



+ 



„ u!^ 2ll' uo f d)'^ 

+ ^ + — + 77 ' 
4 r 2 



+ — + 



z/' 1 
2"^ r 



e 



(7.1.10) 



as the T^^, T^^ and equations, respectively. The propagation equation ()2.1.9p 
now becomes 



87r(p-3p) 
{2uo + 3)0 













a 







~'~ r 



(7.1.11) 



and the only other non-trivial field equation is that for T^'': 

Svr 



,a 

V ijJ 

a 



(p am V 
— + -— H - 

a 2 



(p + p)ve'' 



(7.1.12) 



where v is the three velocity of the perfect fluid. We now assume is of the form 
0(r, t) = 0(r)0(t) and look for solutions to the set of equations 



/ ^ /iV /iV _ c.(r) /0'(r)y 0'(r) / z.' 2\ 

T + ^ + ^"^*^0(O; 0WV 2" + ry'' ^^•^•^''^ 



/fV; z/ ^ /0-(r)y fir^ _ ( l\ ^ 

2 4 2r 2V0(r); 0(r) V 2 r y' 0(r) ' ^ ^ 

. I / I /0-(r)y 0"(r) /z.- 1\ 0-(r) 

^^4^r 2\d){r) Mr) V 2 ^ r i 0(r) ' ^ 
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0(r) 2 0(r) r 0(r) ^ ' ^ 

Such solutions are given by |136j as 



1 ^\ 2fc 

2kr 



1 + — 



e 



(7.1.17) 



^ s 4 /I _ _^^\ 2(fc-l){fe+2)/fc 
M - f 1 I JL\ / ^ 2fcr 



and 
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1 -2{fc2_i)/fc 

-"^ 2fcr 
1 + — 



(7.1.19) 



3+2c<; 

become 



where k = J |±|^. For e'' and of this form, equations (|7T8|l . ()7.1.1()jl and ()7.1.njl 

87r(pe'- - Spe-) 0(t) d0(t) 
(2cu + 3)0(t) 0(t) a0(t)' ^ ■ ^ 

Stt „ /a\^ uj f(i)(t)Y a (hit) 

-pe^ = ^i-\ _ p^lil +3_:^^ (7.1.21) 



and 



0(t)'^ \aj 2 \ (pit) / a0(t)^ 



Stt {?,ujpe- + ?>pe- + ujpe-) _-d a 0(t) ^ u; / 0(t^ \ ' ^ ^2) 



30(t) (2a; + 3) a a(f){t) 3 I 0(t) 

Here (|7.1.22jl was obtained by substituting (fTTKIl and (|7.1.1()jl into (|7.1.njl and 
discarding the terms involving r derivatives, as these now sum to zero. We see that 
^rrm . (frr^ and rrrm are slmply dHH), and (im|l with k = 0, 

Vfrw = V^" and Pfrw = P^^ i where subscript FRW denotes a quantity derived 
from the field equations using the FRW metric. We also have, from T"^;^ = 0, that 

j^{pe'') + ?,H{pe' +pe'') = i). (7.1.23) 

Looking for solutions of the form a oc t^ and (f){t) oc t^ gives, on substitution into 
()7.1.20|) . the relation y = 2 — 3x7 (assuming an equation of state for the Universe 
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of the form p = (7 — l)p). Using ()7.1.2()|1 . (j7.1.21|l . and this relation then gives the 
solutions 

2t^(2 — 1)+2 
f t \ 3"7(2-7)+4 

<t) = «o - (7.1.24) 



and 



to 



2(4-37) 

f \ 3^7(2-7)+4 

0(t)=0o - • (7.1.25) 



M 

These are exactly the same as would be expected for the scale factor and the Brans- 
Dicke field in a flat FRW universe |lH5j . The form of p is then given by ()7.1.23j) and 

dunjas 

\ 37 /I c \ -2k 

'^0 \ / -L - 2fe7 



and 0(r, t) is given as 

This separable solution displays the same time dependence as the power-law 
FRW Brans-Dicke universes, ()3.2.2|l - ()3.2.3j) . but with an additional inhomogeneous 
r-dependence created by the matter source at r = 0. Such a distribution of matter 
in space is illustrated by figure 17.21 Here we have chosen, for illustrative purposes, 
7 = 1, c<j = 100 and a background value set by the choice p = Pfrw {Pfrw being the 
matter density that would be expected in the corresponding homogeneous universe). 
The temporal evolution of p is the same as the FRW case. 

We see from figure 17.21 that the matter density is isotropic and asymptotically 
constant as r cx3 with a sharp power-law peak near the origin. Now ()7.1.27|) gives 
us 

G{r,t)=Goi ^ t-3.7{2-7)+4. (7.1.28) 

+ 2fc7/ 

Equation ()7.1.28p is used, with the values u = 100, 7=1 and c = 0.5 to create 
figure 17^ which shows the space-time evolution of G{r,t). 

These results show how G{r, t) can vary in space and time in an asymptotically- 
flat universe with a peak of matter at the origin. Observers located near the mass 
concentration will determine different values of G locally, although they will find 



97 




10 



15 



20 



r 



Figure 7.2: Distribution of p as a function ofr, withuo = 100, from equation ^7.1.2(^ 
and c = 0.5. 



the same values oi G/G everywhere because of the separable nature of the G{r,t) 
evolution in equation ()7.1.2(ij) . This was also the case for solution (jT.l.fij) given in 
subsection 17.1.11 In the next section we shall consider a more realistic model in 
which both G and G/G are different from place to place. Plots like figure [731 can 
be generated for universes dominated by other types of cosmological fluid and with 
different rates of density fall off with r. 

7.1.2 Spherical Collapse Model 

Matching Two Vacuum FRW Brans-Dicke Universes 

We will now consider a simple model of a spherically symmetric cosmological inho- 
mogeneity produced by matching together fiat and positively curved vacuum FRW 
Brans-Dicke universes. This is a well studied technique, first introduced by Lemaitre, 
for studying the non-linear evolution of overdensities in general relativistic FRW uni- 
verses. The overdense region is modelled as a closed universe that at first expands 
more slowly than the background, before reaching an expansion maximum and col- 
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100 



Figure 7.3: Evolution of G{r,t) in space and time in an inhomogeneous matter 
dominated Universe, with uj = 100, from equation \7.1.2t^ and c = 0.5. 

lapsing back to high density, whilst the background continues to expand. In this 
section we consider vacuum universes only, with p = p = 0. 

For flat vacuum FRW universes the Friedmann equation ()3.1.3|1 gives 

\aj acpb 6 ycpb J 

where ' = = (t>b{t) is the Brans-Dicke scalar field and a is the scale factor in the 

fiat background. For a positively curved {k = +1) region the scale factor is taken to 
be S{t), which satisfies the Friedmann equation for the closed vacuum Brans-Dicke 
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universe: 



^V + ^^-^r^V-A (7130) 

where ' = 0p = 0p(t) and r and /c are the proper time and curvature of this 
perturbed region. 

In matching these two regions at t = to = To we must satisfy the boundary 
conditions 



S'(ro) = a(to), 0p(to) = M^o) 



and 



Prom rj to t 



dr J Q \dt J Q \ dr J Q \ dt 



The function T{ri) is now obtained by integrating Sdrj = dr and T{ri) can then be 
used to obtain S{t). We now require a relation between t and r. For this we proceed 
as in reference jHS] and use the equation of relativistic hydrostatic equihbrium |1()()[ 

nrmiTTTj 

^ = (7.1.31) 

or P + P 

where $ is the Newtonian gravitational potential and r is radial distance. Equation 
()7.1.31|1 is derived under the assumptions that the configuration is static and the 
gravitational field is weak, so $ completely determines the metric; then ()7. 1.3111 
is given by the conservation of energy-momentum for a perfect fluid. We now use 
dr = e^dt, and for a scalar field we have an effective equation of state with p = p 
and p = ^0^. Combining these results gives 

±-ktll (7 132) 



100 



Now (f)b oc a 



^ and (p'p QC S ^, and so with the solutions ()3.2.17|1 and ()H.2.12j) this 



gives 



(7.1.33) 







Equation ()7.1.33j) and the relation Sdri = dr allow us to obtain S{t) from ()3.2.18j] . 
This is done numerically. Now fixing the constants of proportionality together with 
k in equations ()3.2.18|) . ()3.2.17|) and ()3.2.12p . in order to satisfy the boundary con- 
ditions, we find equations for the evolution of the scale factors and scalar fields in 
the fiat background and perturbed region. These are matched at a boundary, at 
time to = ^0 = '7o- 



Figure [73] shows the evolution of a{t) and S{t) when the region described by S{t) 
becomes positively curved at initial time to = 1- Iii figure [73] we choose uj = 100, 
for illustrative purposes, and ao = 5*0 = 1 so that the boundary condition for the 



We can now express G = G(t) in the regions of different curvature using equation 
(I7.1.39jl and the flat and closed FRW solutions (13.2. 17jl and (13.2. 12|) . along with the 
appropriate coordinate transformations. This gives flgure 17.51 It is clearly seen 
that the evolution of G{t) is quite different in the two regions, as expected. The 
collapsing overdensity evolves faster than the background and possesses a smaller 
value of G but a larger value of \G/G\ at all times after the matching. 



We now generalise the spherical collapse model described in the last section to the 
more realistic case of a flat universe containing matter and a cosmological constant 
(see. e.g. refs. |149j . |15()j or |ll(ij ). As before, we match a flat Brans-Dicke FRW 
background to a spherically symmetric overdensity at an appropriate boundary and 
allow the two regions to evolve separately. 



Results 








7.1.3 A More Refined Splierical Collapse Model 
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Figure 7.4: Evolution of the scale factor S in the perturbed overdense region (dashed 
line) and in the background a (solid line) with respect to the comoving proper time 
in the fiat background. 

G(t)/G(5) 
7 



1 2 3 4 5 

Figure 7.5: Evolution of G{t) in the overdense perturbed overdense region of positive 
curvature (dashed line) and in the spatially flat background universe (solid line). 
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The background universe 

Again, we consider a fiat {k = 0), liomogeneous and isotropic background universe. 
Since we are interested in the matter-dominated epoch, when structure formation 
starts, we can assume that our universe contains only matter and a vacuum energy 
contribution so that p = Pm + Pa and p = pa = — pa give the total density and 
pressure, respectively. So, for the flat background, the Friedmann equation ()3.1.3|) 
gives for a general uj{(j)) theory, 

and the scalar field equation ()3.1.2|1 gives 

a [Zuj + 6) [Zuj + 6) 

Here, pm oc a(t)^^ and p\ = constant. These equations govern the evolution of (f){t) 
and a{t) in the fiat expanding cosmological background. The contribution of the 
vacuum energy stress {p = — p) to the Friedmann equation ()7.1.34j) in Brans-Dicke 
cosmology differs from that in GR because of the presence of the variable field: it 
is not the same as the addition of a cosmological constant term to the right-hand 
side of fl7.1.34|) . However, with this proviso, we shall continue to refer to ACDM 
models in Brans-Dicke theories in the following sections. 

The overdensity 

Again, we consider a spherical overdense region of radius 5* and model the interior 
space-time as a closed FRW Brans-Dicke universe, ignoring any anisotropic effects of 
gravitational instability or collapse. As usual, we assume there is no shell-crossing; 
this implies mass conservation inside the overdensity |149j . The evolution equations 
can now be written in a form that ignores the spatial dependence of the fields. Now 
P = PcDM + Pa and p = Pa = — Pa, where 'CDM' corresponds to cold dark matter. 
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so that the acceleration equation gives 



3 02 2(2tUe + 3) 



, Stt (pcDAf (3 + ^c) + Pa(3 - 2a;,)) ^^ 

+3^ (2:;-r3^ j- 

while the scalar field equation ()3.1.2|) reduces to 

where S = S{t) is the scale factor and 0c = 0c (^) is the Brans-Dicke scalar field in 
the collapsing region of positive curvature where pcdm oc S{t)~^ and pA = constant. 
These equations give the evolution of 0c(t) and S{t). 

We have assumed that the equation of motion of the field inside the cluster 
overdensity is described by the local space-time geometry. This means that the field 
follows the dark-matter collapse from the beginning of the cluster's formation. We 
do not consider this to be fully realistic since there is expected to be an outfiow 
of energy associated with from the overdensity to the background universe, as 
first noticed by Mota and van de Bruck ( |13nj ). The details of this outflow of 
energy and its effect on the collapse can only be determined by a fully relativistic 
hydrodynamical calculation, which is beyond the scope of this study. Nevertheless, 
at late times during the collapse of the dark matter (and especially when the density 
contrast in the dark matter is very large) the field should no longer feel the effects 
of the expanding background and will decouple from it. We are also neglecting the 
effects of deviations from spherical symmetry, which grow during the collapse in the 
absence of pressure, along with rotation, gravitational tidal interactions between 
different overdensities, and all forms of non-linear hydrodynamical complexity. 

Evolution of the overdensity 

Consider a spherical perturbation in the dark matter fluid with a spatially constant 
internal density. Initially, this perturbation is assumed to have a density amplitude 
(5j > where \5i\ <^ 1. The initial density of dark matter inside the overdensity is 
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therefore pcdm = Pm(l + 

Four characteristic phases of the overdensity's evolution can be identified: 

• Expansion: we employ the initial boundary condition (pc = (p assume that 
at early times the overdensity expands along with the background. 

• Turnaround: for a sufficiently large 6i, gravity prevents the overdensity from 
expanding forever; the spherical overdensity breaks away from the general 
expansion and reaches a maximum radius. Turnaround is defined as the time 
when S = Smax, S = and S < 0. 

• Collapse: the overdensity subsequently collapses {S < 0). If pressure and 
dissipative physics are ignored the overdensity would collapse to a singularity 
where the density of matter would tend to infinity. In reality this singularity 
does not occur; instead, the kinetic energy of collapse is transformed into 
random motions. 

• Virialisation: dynamical equilibrium is reached and the system becomes sta- 
tionary with a fixed radius and constant energy density. 

We require our spherical overdensity to evolve from the linear perturbation 
regime at high redshift until it becomes non-linear, collapses, and virialises. 

Virialisation 

In scalar-tensor theories we expect that the gravitational potential will not be of 
the standard local form. This requires reconsideration of the virial condition. 
According to the virial theorem, equilibrium will be reached when jOHl 



where T is the average total kinetic energy, U is the average total potential energy 
and S here denotes the radius of the spherical overdensity. 

The potential energy for a given component x can be calculated from its general 
form in a spherical region |117j 





X 
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where ptot is the total energy density and $2,. is the gravitational potential due to 
the density component px- 

The gravitational potential can be obtained from the weak-field limit of the 
field equations ()2.1.8p . This results in a Poisson equation where the terms associated 
to the scalar field can be absorbed into the definition of the Newtonian constant as 
[Ml 

G.= i±^l. (7,1,39) 
This results in the usual form for the Newtonian potential 

{s) = -27rGeP.(37. - 2) (^S^ - ^ 

where Gc is given by equation ()7.1.39|) and — 1 is Px/ Px for the fiuid component 
with density px and pressure px (appearing due to the relativistic correction to 
Poisson's equation: A$ = AnG (p + 3p)). 

In ACDM models of structure formation it is entirely plausible to set 7^. = 1 
as the energy density of the cosmological constant is negligible on the virialised 
scales we are considering |188t ll3Uj . The potential energy associated with a given 
component x inside the overdensity is now given by 

If) 71-2 

Ux = -^G.ptotPxS''. (7.1.40) 
Therefore, the virial theorem will be satisfied when 

T- - -S- ('^ 



dS 



where 



dU 167r2 



dS 15 



^5 ^ 9 Ptot ^5 ^ dpx q5 . ^ p- ri4 

PtotPxiJ + <-Jc r^o PxlJ + '^cPtot-^^'J + '^cPtotPx^^ 



(7.1.41) 

and we have used Utot = Ucdm + Ua + U^^, ptot = Pcdm + Pa + P</.,, P4>c = ^c4>^l<P 
and 

dG, G, ^ J, 3 + 2uj, , 2071(0,) 



= -G, 0c— ^ I a 



dS S ^ A + 2uj, \ (3 + 2cjc)2 
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The other components of equation ()7.1.4H1 are obtained from equations ()7.1.Hfj|l and 



Using equation ()7.1.38|) . together with energy conservation at turnaround and 
viriahsation, we obtain an equihbrium condition in terms of potential energies only 



where is the redshift at viriahsation and is the redshift of the over-density at 
its turnaround radius. The behaviour of G during the evolution of an overdensity 
can now be obtained by numerically evolving the background equations ()7.1.H4|1 and 



fl7.1.H5|l and the overdensity equations ()7.1.H(j|l and ()7.1.H7|1 until the virial condition 
flTTi^ holds. 



We point out here an inconsistency when one makes use of equation ()7.1.42|) to- 
gether with the assumption that energy is conserved. This inconsistency is removed 
by assuming a negligible outflow of from the overdensity, in which case we regain 
energy conservation within the system and so retain self-consistency. 

Overdensities vs Background 
Brans-Dicke theory 

The Brans-Dicke coupling parameter uj is constant and constrained by a variety of 
local gravitational tests (see |18(ij for a review). The strongest constraint to date 
is derived from observations of the Shapiro time delay of signals from the Cassini 
space craft as it passes behind the Sun. These considerations led Bertotti, less and 
Tortora after a complicated data analysis process, to claim that uj must have 
a value greater than 40000, to 2a. This limit on uj must be satisfied at all times in 
all parts of the universe, and leads to the conclusion that Brans-Dicke theory must 
be phenomenologically very similar to GR throughout most of the history of the 
universe. However, we do still expect a cosmological evolution of the Brans-Dicke 
field which determines the value of Newton's G\ and we expect this evolution to be 
different in regions that collapse to form the structure probed by Cassini compared 
to that in the idealised expanding cosmological background, as described above. 
Hence we expect the measurable value of G to be different in these two distinct 




(7.1.42) 
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Figure 7.6: Plots of G against ln(l + z) in the background (dashed-line) and in 

an overdensity (solid-line), for different values of uj. Initial conditions are chosen 
in both cases so as to give G = Go, the present value of the Newton constant, at 
virialisation. We note that increasing uj decreases the difference in G between the 
overdensity and the background. 
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regions with different histories. 

The plots in figure EEl were constructed using the representative values u = 40000 
and u = 500 and the boundary condition 0^0 — Gq^, so that the value of G measured 
inside the overdensity at present is equal to the value of Newton's G, as measured 
locally. The evolution of the background was determined by matching 0c to (p at 
the time when the overdensity decouples from the background, ti. We see a clear 
difference in the evolution of G in the two regions, as expected. This example 
shows that we expect different values of G and G/G inside and outside virialised 
overdensities. The present value of G and G/G depends on the history of the region 
where it was sampled, as well as on the Brans-Dicke coupling parameter, u. 

It can be seen from the plots in figure I7.fil that increasing u has the effect of 
decreasing the difference in G between the background universe and the overdensity. 
The size of this inhomogeneity is found to be of order l/u and, correspondingly, 
reduces to zero as ^ oo. This is an important consistency check for the methods 
used as we expect Brans-Dicke theory to reduce to GR, with a constant G, in this 
limit. 



Scalar-tensor theory with 2uj + 3 = 2A 



1 



Next, we consider a scalar-tensor theory with a variable ui{(j)). We investigate the 
class of theories defined by the choice of coupling function 



2u{(j)) + 3 = 2A 



1 - 



where A, 0oo, and p are positive definite constants. We refer to this as Theory 1. 
Such a choice of coupling was considered by Barrow and Parsons and was solved 
exactly for the case of a fiat FRW universe containing a perfect fluid 

Setting the constants as 2A = (0oo//5)^ and p = 2 gives us the scalar-tensor the- 
ory considered by Damour and Pichon |77] and by Santiago, Kalligas and Wagoner 
|167j . This choice of uj{(f)) corresponds to setting lnA(0) = lnA(0o) + \P{<Poo — 
0)^, where A'^{(j)) is the conformal factor 1/0 from equation ()2.2.fij) . Damour and 
Nordvedt (TH] consider this function as a potential and therefore justify its choice 
in relation to the generic parabolic form near a potential minimum. Expecting the 
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(a) 2tt;(<^)+3 for p=2 and A=l, 2 and 5 (solid, dashed and dotted 
lines, respectively). 
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(b) 2Lj{(f)) + 3 for A=l and p=1.5, 2 and 2.5 (solid, dashed and 
dotted lines, respectively). 



Figure 7.7: In these gravity theories there is fast approach to general relativity at 
late times when —>■ 4>^, but significantly different behaviour at early times. 
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function to be close to zero (i.e. GR), Santiago, Kalligas and Wagoner justify its ex- 
pression as a perturbative expansion. This choice of uj{(f)) with p > 1/2 corresponds 
to a general two-parameter class of scalar-tensor theories that are close to GR and 
will be drawn ever closer to it with u oo and u' /u^ as the universe expands 
and (j) (poo- We therefore consider it as a representative example of a wide family 
of plausible varying-G theories that generalise Brans-Dicke theory. 

The evolution of this form of u;{(f)) is shown graphically in figure FTTI for different 
values of A and p. Clearly the evolution of uj{(j)) is sensitive to both A and p and so 
the choice of these parameters is important for the form of the underlying theory. 
For illustrative purposes we choose here the values p = 1.5, 2 and 5 and A = 1,2 
and 5. 

In a similar way to the Brans-Dicke case we now create an evolution of 0c that 
virialises at 2; = to give the value Gco = 6.673 x 10^^^, as observed experimentally. 
The corresponding evolution for (pb is calculated as before by matching it to the 
value of (pc at the time the overdensity decouples from the background and begins 
to collapse. In creating these plots we have used the conservative parameter values 
p = 2, uJco = 1-2 X 10^ and A = 6 x 10~^ which are consistent with observation and 
allow structure formation to occur in a similar way to GR. The results of this are 
plotted in figure TTM 

Again, we note the different evolution of G{t) in the two regions, and the differ- 
ence in the asymptotic values of G. We note that experimental measurements of G 
on Earth have a significant uncertainty with the 1998 CODATA value carrying an 
uncertainty 12 times greater than the standard value adopted in 1987. The 1998 
value is given as |127^ Il(i9j 

Gi998 = 6.673 ± 0.010 X 10~^ cm^gm'h"^, 

while the 2002 CODATA pre-publication announcement reverts to the earlier higher 
accuracy consensus with |138j 

G2002 = 6.6742 ± 0.0010 x IQ-^cm^gm-^s-^ 

We could re-run the above analysis with different values of A and p, but expect 
that the results would look qualitatively similar. From figure 17.71 we see that in- 
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Figure 7.8: 2uj + 3 — 2 A 1 — ^ A graph of the variation G ag ainst ln(l + z) for 

the background universe (dashed-line) and inside an overdensity (solid-line) which 
give G — Go,the presently observed terrestrial value, at virialisation. 
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creasing (decreasing) the values of A and p will increase (decrease) the value of uj{(j)) 
for a given 0, thereby making the theory more (less) like GR. We therefore expect 
an analysis with a higher (lower) value of A and/or p to look very similar to the 
analysis presented above with a less (more) rapid evolution of G(t). For the sake of 
brevity we omit such an analysis here. 

Space and Time variations of G 

We now calculate how time and space variations of G evolve with redshift and 
depend on the dark matter density contrast, Ag. In order to do this, we make the 
definitions 

G Gq G Gi)(t) 

and 

A _ Pcdm{Zv) 



Pm ) 

where Gc and Gb correspond to G as measured in the overdensity and in the back- 
ground universe, respectively. 

The results of our numerical calculations, for a cluster which virialises at ^^, = 0, 
are presented in figures 17.91 17.101 and 17.111 respectively. These plots display the 
evolution of G/G, AG/G, and 6G/G with redshift for Brans-Dicke theory, the theory 
of subsection 17. 1.3| and some other choices of 0^(0) that are specified in the captions. 
The parameters used in generating these plots are B = 0.4, G = 10^^^, D = 80, 
A = 6 X 10"'', p = 2 and u = 4x 10^. These values were chosen so as to agree with 
observation and so that structure formation is not significantly different from that 
which occurs in GR. 

It is clear from the plots that different scalar-tensor theories lead to different 
variations of G. The predictions of these models can be quite diverse. While some 
models produce higher values of G inside the overdensity, others produce a lower 
one. A feature common to all models is that the value of G and G/G inside an 
overdensity is different from G and G/G in the background universe. The reason for 
these differences is that in the non-linear regime, when the overdensity decouples 
from the background expansion at turnaround, the field that drives variations 
in the Newtonian gravitational 'constant' stops feeling the background expansion. 
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Figure 7.9: Evolution of AG/G as a function of log {1 + z) for overdensities which 
virialise at z = in a ACDM model. Upper plot: evolution inside the overdensity. 
Lower plot: evolution in the background universe. Thick solid line 2uj + 3 — Ax 10^, 
thin solid line 2uj + 3 — B'^cj), dashed-line 2uj + 3 — 2A\1 — ^|~^; dash-dotted line 

2a; + 3 = C| ln(^)|~^, dotted line 2uj+3 — D\l — (-^y\~^ . Each model is normalised 
in order to have Gq — Gc{z — 0) inside the overdensities. 
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Figure 7.10: Evolution of G/G as a function of log {1 + z) for overdensities which 
virialise at z = in a ACDM model. Upper plot: evolution inside the overdensity. 
Lower plot: evolution in the background universe. Thick solid line 2u! + 3 — Ax 10^, 
thin solid line 2a; + 3 = B'^cf), dashed-line 2uj + 3 — 2A\1 — ^|~^; dash-dotted line 

2a;-|-3 = C| ln(^)|~^, dotted line 2a;+3 = D\l — (-^y\~^ . Each model is normalised 
in order to have Gq — Gc{z — 0) inside the overdensities. 
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log (1+z) 

Figure 7.11: Evolution of 5G/G as a function of log{l + z) for overdensities which 
virialise at z = in a ACDM model. Thick solid line 2^; + 3 = 4 x 10®, thin 
solid line 2u; + 3 = ^^^^{p+i)^ dashed-line 2lu + 3 = 2A\1 - ^|"^ dash-dotted line 

2ct; + 3 = C| ln(^)|~'', dotted line 2uj + 3 = D\l — {-^y\~^ . Each model is normalised 
in order to have Gq = Gc{z = 0) inside the overdensities. 



116 




Figure 7.12: Values of AG/G as a function of log {1 + z^) at virialisation. Upper 
plot: evolution inside the overdensity. Lower plot: evolution in the background 
universe. Thick solid line 2a; + 3 = 4 x 10®; thin solid line 2uj + 3 = B'^cf); dashed- 
line 2a; + 3 = 2^411 — ^|~^; dash-dotted line 2a; + 3 = C|ln(^)|~^; dotted line 

2a; + 3 = D|l — (^)^|~^ • Each model is normalised to have Gq — Gc{z — 0) inside 
the overdensities. 
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After turnaround, the field inside the overdensity, 0c, deviates from the field, 0, in 
the background universe, leading to spatial variations in G. In reality, such spatial 
inhomogeneities in the value of G are small: SG/G ~ 10^^, figure ITTTl The time 
variations of G are even smaller than the spatial inhomogeneities but with a marked 
difference between the inside and outside rates of change. We find Gc/Gc < 10^^*^s^^ 
inside the clusters and G/G < lO^^^s^^ in the background, figure mill 

Figures TTTA TTm and Ell represent the values of AG/G, 5G/G and G/G at 
the redshift of virialisation (z^). Once again the differences between the different 
scalar-tensor theories and between Gc{zy) and G{z^) are evident. 

Figure mHl shows how the dark matter contrast, Ac = Pcdm{zv) / Pm{zv) , affects 
the difference between the value of G inside an overdensity and in the cosmological 
background. (Recall that, in the Einstein-de Sitter model Ac = PcDAiizy)/ pm{zy) ~ 
147 at virialisation, and Ac = PcDuizy)/ pm{zc) ~ 187 at collapse). It is interesting 
to see that different scalar-tensor theories produce a different dependence. 

From the numerical simulations, we see that a variation of the Newtonian gravi- 
tational constant, of the order presented here, does not affect the predictions of the 
structure formation models. For instance, the virialisation radius and the density 
contrast Ac of the virialised clusters are very similar to the ACDM and standard 
CDM models with constant G . Nevertheless, one should point out that different 
scalar-tensor theories lead to slightly different results. Besides that, there is also a 
dependence on the initial conditions, as with many other cosmological scalar fields 
(see e.g. |13m 112!^ I129p . Different initial conditions will lead to slightly different 
allowed values of u and to different cosmological behaviours of G and Gc- 

7.1.4 Discussion 

We have studied the inhomogeneous cosmological evolution of the Newtonian grav- 
itational 'constant' G within the framework of relativistic scalar-tensor theories of 
gravity, of which Brans-Dicke theory is the simplest and best known case. We be- 
gan by first exploiting the conformal equivalence between these theories and GR to 
transform an existing solution of Einstein's field equations to a new exact spherically 
symmetric inhomogeneous vacuum cosmological solution of the Brans-Dicke field 
equations. We then present a second spherically symmetric perfect-fiuid solution to 
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Figure 7.13: Values of G/G as a function of log(l + z^) at virialisation. Upper 
plot: evolution inside the overdensity. Lower plot: evolution in the background 
universe. Thick solid line 2a; + 3 = 4 x 10®; thin solid line 2uj + 3 = B'^cf); dashed- 
line 2a; + 3 = 2^411 — ^|~^; dash-dotted line 2a; + 3 = C|ln(^)|~^; dotted line 

2a; + 3 = D|l — (^)^|~'^- Each model is normalised to have Gq — Gc{z — 0) inside 
the over densities. 
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Figure 7.14: Values of 5G/G as a function oflog{l + Zy) at virialisation. Thick solid 
line 2u;+3 = 4x10^; thin solid line 2a;+3 = 5^0; dashed-line 2oj +3 = 2A\l--^\-'^; 

dash-dotted line 2a; + 3 = C| ln(^)|'^; dotted line 2u; + 3 = ^ll - (^)T^. Each 
model is normalised to have Gq — Gc{z — 0) inside the overdensities. 
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Figure 7.15: Values of 6G/G as a function of for a ACDM model. Thick-solid 
line 2a; + 3 = 4 X 10^- solid line 2a; + 3 = B'^cp; dashed-line 2a; + 3 = 2A|1 - 

dash-dotted line 2a; + 3 = C| ln(^)|'^; dotted line 2a; + 3 = D]! - (^)T^. Each 
model is normalised to have Gq — Gc{z — 0) inside the overdensities. 
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the Brans-Dicke field equations which corresponds to an overdense inhomogeneity 
in an asymptotically homogenous and isotropic expanding cosmological background. 
These exact solutions have simple enough form to allow the G{r,t) evolution to be 
investigated directly and can be used to model the presence of a Schwarzschild-like 
mass in an expanding Brans-Dicke universe that approaches an idealised FRW uni- 
verse asymptotically as r — > oo. It is noted that far from the mass G varies very 
slowly with r and that as r — > oo the variation of G with r is removed altogether 
and only the cosmological evolution of G with time is relevant. Close to the mass, 
the r variation of G becomes more significant and we see that G — > at small r. 
We also note that in the limit a; — > oo all space and time variation in G is removed 
from these solutions and GR is recovered. 

Next, we increased the complexity of the model by matching together two vac- 
uum FRW Brans-Dicke universes of zero and positive curvature on a spacelike time 
slice. This results in a simple model for a spherically symmetric perturbation in 
the density and curvature that is followed into the non-linear regime. The different 
evolutions of G in the two regions were determined and a comparison was made of 
the different evolutions of G on spacelike slices of differing time. We see from this 
toy model that the value of G does indeed have a different value in regions that have 
decoupled from the expanding cosmological background and collapsed, compared to 
its value in the background itself. We were then able to repeat this construction 
for a more realistic matching of two FRW ACDM universes of different curvature in 
an arbitrary scalar-tensor gravity theory. We followed the development of spherical 
overdensitics through their expansion, separation from the background, turnaround, 
and subsequent collapse. Applying a simple approximation to virialise the collaps- 
ing overdensitics, wc were then able to study the differences in G{x^^t) between 
the non-expanding ovcrdcnsity and the expanding background universe. Wc high- 
lighted as a special case the simple Brans-Dicke theory but also presented results 
for other scalar-tensor gravity theories, specified by their defining coupling function 
uj{(j)). Although each theory predicts a different detailed cosmological evolution of 
G, a feature common to all of them is the difference between the value and time 
evolution of G inside bound overdensitics and in the background universe. These 
differences will produce spatial inhomogeneities in G with a value which depends 
on the scalar-tensor theory used. While some models produce higher values of G 
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inside the overdensity, others produce lower ones. In spite of these differences, such 
spatial inhomogeneities are small: SG/G ~ 10^^. The differences in the time vari- 
ations of G were found, with typically Gc/Gc < lO^^^s^^ inside the clusters and 
G/G < lO^^^s^^ in the background universe in the range of theories studied. Such 
variations in G do not significantly alter the virialisation radius and the density con- 
trast Ac of the virialised clusters from those in standard ACDM and CDM models 
with constant G. Nevertheless, different scalar-tensor theories lead to slightly dif- 
ferent results. There is also a dependence on the initial conditions. Different initial 
conditions will lead to a different value of u and to different cosmological behaviours 
of G and Gc- Taken as a whole, these analyses show that local observational limits 
on varying G made within our solar system or Galaxy must be used with caution 
when placing constraints upon the allowed cosmological variation of G on extra- 
galactic scales or in the early universe. The universe is not spatially homogeneous 
and, in cosmological models where it can vary, nor is G. 

7.2 Fourth-Order Theories 

It has been known for some time that Birkhoff's theorem is not ensured in the 
vacuum solutions of generalised fourth-order theories of gravity, as they are in GR 
(see e.g. |151U18l] ). We will now present an exact solution that shows this behaviour 
explicitly; a non-static, vacuum solution admitting spherically symmetric three- 
dimensional space-like hypersurfaces. A perturbation expansion will be performed 
about the background (r — > oo limit) of the exact solution. 

7.2.1 Exact Solution 

A solution to the spherically symmetric vacuum field equations of i?" theory 1)6.2.111 
is given by 

ds'' = -A2{r)dt^ + a'{t)B2{r){dr^ + r^dn^) (7.2.1) 
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where 




5+2(5-1 



a(t) 




1-25 + 45^ 



and 6*2 is a constant. This solution is conformally related to one found by Fonarev 
jHE] and reduces to the Schwarzschild solution in the limit 5 — > 0. It shows explicitly 
the lack of validity of Birkhoff 's theorem in these theories. It becomes r-independent 
in the limit r — > oo, but still displays strong curvature in this limit as the metric 
reduces to the spatially fiat vacuum FRW cosmology found in [70] . 

We will now continue to find the general form of spherically symmetric pertur- 
bations to the background of the above solution. It will be seen that there exist 
extra modes which are not excited in the exact solutions, but that the modes cor- 
responding to the linearised exact solutions above are the ones most interesting for 
performing gravitational experiments in this space-time. 

7.2.2 Linear perturbation analysis 

Perturbative analyses in the literature are often performed about Minkowski space 
or de Sitter space. This is perfectly acceptable practise in GR and fourth-order 
theories in which an Einstein-Hilbert term dominates in the low curvature regime. 
In other fourth-order theories, of the type considered here, in which the Einstein- 
Hilbert term does not dominate the low curvature regime then there is good reason 
to consider perturbing about other backgrounds. We have shown explicitly, with 
exact solutions, the existence of other spherically symmetric space-times. We will 
now proceed to perform a linear perturbation analysis about the background of 
solution ()7.2.H) . The general solution to first order in perturbations will be found. 
Writing the perturbed line-element as 



allows the vacuum field equations to be linearised in P and Q. These linearised field 



ds^ = -(1 + P{r))de + h\t){l + Q{r)){dr^ + rW) 



(7.2.2) 
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equations are 



6S{1 - 26- 26^){5 - U6 - 126^)VP 

- 12(1 - 26 - 26^){1 -Q6 + A6^ + A6^)VQ 



^ 2(1-2(5-2(5^) 

+ (1 + 25)(1 - 6)h (1-*-) 



(2 - 175 + 185^ + 525^ + 86^)P" -2(1-6 + 36^ - 346^ - 325^) — 

r 

+ 2(2 - 86- 96^ + 166^ + 86^)Q" + 2(1 - 45 - 155^ + 206^ + 1Q6^) — 

r 



2(l-2S-2S^) ( tJj" tJj' 

+ 3(1-^)^(^^(^I^ + |L 

and 

V^' = 

where 

il) = VP + 2VQ. 

Primes here denote differentiation with respect to r and V is the Laplacian on the 
three-dimensional subspace. These equations have the general solution 

r (5 - 145 - 1252) 

^ (l-25)c4 , 2 

Q = h ocsr + constant 

r 

where C4 and C5 are constants and the two other constant terms in P and Q are 
independent of each other and can be absorbed into t and s by redefinitions. It can 
be seen that one of the modes, C4, corresponds to the linearised version of the exact 
solution (fTTT^ . 

We immediately see that the form of these linearised solutions are quite different 
to those obtained by expanding around the background of the static spherically 
symmetric solution ()(i.2.2|l . Whist the expansion about ()(i.2.2j) produces damped 
oscillatory modes, as well as the mode corresponding to the linearised exact solution, 
the expansion about ()7.2.1|) produces more familiar looking terms proportional to 
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r^. Aside from the different form of these extra modes, there is also a noticeable 
difference in the terms corresponding to the linearised exact solutions, which both 
go as r^^ in the limit 5 — > 0, but behave differently from each other when 5 7^ 0. 
This shows explicitly the differences that can arise when linearising about different 
backgrounds. Not only are there extra modes which can take different functional 
forms, but even the modes which reduce to the Schwarzschild hmit as 5 ^ are 
different, depending on the background that has been chosen. 
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CHAPTER 8 



Scalar- Tensor Cosmologies with Energy Exchange 



Non-minimally coupled scalar fields arise in a variety of different theories, including 
Kaluza-Klein theory IllOj , string theories jHl] and brane- worlds I157j . The 
same mechanism that creates a scalar field non-minimally coupled to the curvature 
in these theories can also lead to a coupling between the scalar and matter fields. 
This coupling manifests itself through the matter Lagrangian becoming a function of 
0. The possibility of such a coupling is usually neglected in the literature, where the 
matter Lagrangian is a priori assumed to be independent of 0. It is the possibility 
of a coupling between the scalar and matter fields in scalar-tensor theories that will 
be the subject of this chapter. 

The introduction of a coupling between (p and matter greatly enlarges the phe- 
nomenology of the theory. Potentially, this allows greater variability of G in the early 
universe whilst still satisfying the solar system bounds on varying G [47 \. We will 
investigate the extent to which G can vary when energy is exchanged between the 
field and ordinary matter. As well as giving a window into the four- dimensional 
cosmologies associated with higher- dimensional theories, we hope that this direction 
of study might be useful in understanding why the present value of G is so small 
compared to the proton mass scale {Grrip^ ~ 10^^^). The direct exchange of energy 
between and the matter fields offers a non-adiabatic mechanism for G to 'decay' 
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towards its present value from a potentially different initial value. There have been 
a variety of studies which investigate the drain of energy from ordered motion by en- 
tropy generation, due to bulk viscosity [13 HZl lEHI or direct decay |2H1 EZj or energy 
exchange |183j . but few studies of the drain of energy by non-adiabatic processes 
from a scalar field that defines the strength of gravity [HIH EHl IHEj- This creates a 
range of new behaviours in scalar-tensor cosmologies. 

In considering a coupling between and matter we are forced to reconsider the 
equivalence principle. The energy-momentum tensor of perfect-fluid matter fields 
will no longer be covariantly conserved and the trajectories of test-particles will no 
longer follow exact geodesies of the metric. These violations of the experimentally 
well verified weak equivalence principle exclude most possible couplings between 
and matter |191j . Such violations are not necessarily fatal though. We show that 
whilst energy-momentum is not separately conserved by the matter fields there is 
still an exact concept of energy-momentum conservation when the energy density of 
the scalar field is included. Furthermore, the non-geodesic motion of test particles 
is only problematic if the coupling increases above experimentally acceptable levels 
as the Universe ages. The theory we consider is still a geometric one and it remains 
true that at any point on the space-time manifold it is possible to choose normal 
coordinates so that it looks locally flat. 

This chapter is based on the work of Clifton and Barrow [7Hj . 

8.1 Field Equations 

The conservation of T"'', whilst appealing, is not absolutely necessary in deriving a 
theory in which G can vary. There are numerous examples where one might expect 
T"-^ 7^ 0. For example, when considering two fluids the energy-momentum tensor of 
each fluid is not separately conserved unless the fluids are completely non-interacting 
|24j . It is only required that the energy- momentum being lost by one of the fluids 
is equal to the energy-momentum being gained by the other. 

In what follows we will consider the scalar and matter fields as two fluids (or 
more than two fluids if there is more than one matter fluid present) and introduce 
a transfer of energy and momentum between them. Such an interaction can be 
introduced by allowing Cm to be a function of (p and will change the nature of the 
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resulting FRW cosmologies. 

The field equations take their usual form in the Brans-Dicke theory ()2.1.8j) and 
the scalar-field propagation equation and matter energy-momentum conservation 
equations are now 

□^ = r^-7^4 (8.1.1) 
T'^^b = a" (8.1.2) 

where T is the trace of the energy-momentum tensor and a" is an arbitrary vector 
function of the space-time coordinates that determines the rate of transfer of 
energy and momentum between the scalar field and the ordinary matter fields. 
The precise form of a"" depends on the detailed form of the interaction between the 
scalar and matter fields in Cm- For example, a conformal transformation of the form 
Qab — > A^{(t))gab from a frame in which T"'\ = gives 

This particular choice of energy transfer can be interpreted as a space-time variation 
of the rest masses of matter described by Cm- For the moment, we consider the case 
of more general interactions by leaving a" as an arbitrary function. Later, we will 
consider specific forms of a"' that allow direct integration of the field equations. 

Substituting the spatially flat FRW metric into the field equations ()2.1.8j) . ()8.1.1|1 
and ()8.1.2j) gives the generalised Friedmann equations: 

+ 3(3 + 2cu)a- (8.1.3) 



3a((2-7)cu + 1) + 3-ta° (8.1.4) 





pV P(p 

M: = «(4-37)-2^t 
PV pet) 

V_ 

'v 



^^'^): = a(4-37)-2^ta° (8.1.5) 
p + 7^P = c^° (8.1.6) 
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where we have defined a comoving volume V = and a constant a = (3^20;) ' 
energy-momentum tensor is assumed to be a perfect barotropic fiuid with density 
p and pressure p which are hnked by a hnear equation of state j9 = (7 — l)p, and 
over-dots denote differentiation with respect to the comoving proper time, t. It is 
this set of differential equations that we need to solve in order to determine the 
evolution of a{t) and G oc 4>{t)~^ in cosmological models of this type. 



8.2 Transfer of Energy and Entropy 

The conservation of energy and momentum, as well as the second law of thermody- 
namics, are of basic importance to physics. In considering an interaction between a 
gravitational scalar field (p and the matter fields \E', it is therefore necessary to inves- 
tigate the extent to which we can consider energy and momentum to be conserved 
and the second law to be obeyed. 

When we consider the thermodynamics of an exchange of energy between the 
scalar field and matter it is useful to define an effective energy density, p^, for the 
scalar field 0. Defining 

lblTUJ(p 

the scalar-field propagation equation ()8.1.4|) can then be rewritten as 

p<^ + 2^p^ = -A0_^o_ (8 2.2) 

Comparison of this equation with ()8.1.6p shows that acts as a fiuid with equation 
of state 7 = 2 (pfj) = p^) . The two terms on the right hand side of this equation act as 
sources for the energy density p^. The first is the standard Brans-Dicke source term 
for the scalar field and the second, cr^{t), is new and describes the energy exchange 
between and the matter fields. It can be seen that the second term is exactly the 
opposite of the source term in equation ()8.1.3|) . and it is in this sense that the total 
energy is conserved in this theory. 

It is also useful to consider the entropy. Contracting the divergence of the energy- 
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momentum tensor with the comoving four-velocity we obtain 



where, in the last hne, we have used the normalisation — —1 and {U"-Ua)-b — 

U°';bUa + U'^Ua-i, — 0. Defining the particle current by N"' = nU"', where n is the 
number density in a comoving Lorentz frame, this expression can be rewritten as 



P;a - n 



{p + p) 



n 



n 



= -nU" 



P 



where we have used the conservation of particle number, N"'-^a = 0. Recalling the 
first law of thermodynamics, 

QdS = pdV + dE^pd (^^^ + ^ (0 ' 
where © is the temperature and S is the entropy, we now get 



or, making use of our assumption of spatial homogeneity. 

This tells us that, as energy is transferred from (f) to the matter fields, the entropy of 
the matter fields increases, as expected. Conversely, the matter fields can decrease 
their entropy by transferring energy into 0. 

Unfortunately, there is currently no known way of defining the entropy of a non- 
static gravitational field so it is not possible to perform an explicit calculation of 
the entropy changes in and gab- We can only assume that if the Universe can be 
treated as a closed system, and the exchange of energy is an equilibrium process. 
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then the entropy that is lost or gained by the matter through this exchange will be 
gained or lost by the gravitational fields. This direct interaction of the matter with 
then allows an additional mechanism for increasing or decreasing the entropy of 
the matter content of the Universe. 

8.3 General Solutions 

It is convenient to define a new time coordinate r by 

dr = pVdt (8.3.1) 
and to re-parametrise the arbitrary function by 

= pV^A' (8.3.2) 


where a prime denotes differentiation with respect to r and A(r) is a new arbitrary 
function. This re-parameterisation of the interaction is chosen to enable a direct 
integration of the field equations and does not imply any loss of generality, as A is 
an arbitrary function. The field equations ()8.1.4j) and ()8.1.5|1 can now be integrated 
to 

p^pVV' = 3«((2 - 7)u; + l)(r - n) + 3«A (8.3.3) 
p\/y = a(4-37)r + ar2-2aA (8.3.4) 

where ri and T2 are constants of integration. We have a freedom in where we define 
the origin of r and can, therefore, absorb the constant ri into r and the definition 
of T2 by the transformations r ^ r + ri and r2 ^ r2 — (4 — 37) ri. It can now 
be seen from fl8.3.4|) that 0' is sourced by three terms. The first corresponds to 
the source term in ()8.1.2|1 and can be seen to disappear for 7 = 4/3, as expected 
for black-body radiation. The second term is constant and is the contribution of 
the free scalar to the evolution of 0; it is this term which distinguishes the general 
spatially-fiat Brans-Dicke FRW solutions [HZ] from the power-law late-time attractor 
solutions |137j . The third term is new and gives the effect of the energy transfer 
on the evolution of (p. This term is dependent on the arbitrary function A, which 



132 



specifies the interaction between and the matter fields. 

The problem is now reduced to solving the coupled set of first-order ordinary 
differential equations ()8.3.3|) and ()8.3.4|) with the constraint equation ()8.1.6|) . The 
remaining equation ()8.1.3p is rewritten in terms of r and A as 

-+ly=^^. (8.3.5) 

and can be solved for p once V and have been found for some A. 

We can decouple the set of equations ()8.3.3j) and ()8.3.4|1 by differentiating ()8.3.4|1 
and substituting for ()8.3.3j) to get the second-order ordinary differential equation 

((4 _ 3^)^ + ^2 _ 2A) ^ - ((4 - 37) - 2A') ^ 

(j) (j) 

= - [((4 - 37)r + r2 - 2A)A' + 3(2 - 7)(((2 - 7)0; + l)r + A)] 
which can be integrated to 



(8.3.6) 



£ _ (4 - 37)r + T2 - 2A 
where 

A = 3-f^uj/2 - 37(1 + 2lo) + (5 + 6uj) 
B = T2 + {4- 37)A 
C = -A^ + Xt2 + D 

and D is a constant of integration. The three source terms for (p' appear in the 
numerator on the right-hand side of equation ()8.3.(j|l . The equations ()8.3.3j) and 
()8.3.4|) can now be combined to give a' /a in terms of 0'/0 as 

a!_ _ ((2-7)a; + l)r + A 
a~ {At^ + Bt + C) 

where A, B, C and D are defined as before. The constant D can be set using the 
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constraint equation ()8.1.(i|l . Using ()8.3.3j) and ()8.3.7|1 we obtain the expression 



pV^cj) = aiAr'^ + Bt + C). 



This can then be substituted into ()8.1.fi|l which, in terms of r and a, gives the 
generahsed Friedmann equation: 



Substituting ()8.3.(i|l and ()8.3.7j) into this, we find that 



8.4 Particular Solutions 

If A is specified in terms of r, we now have a set of two decoupled first-order ordinary 
differential equations for the two variables a and 0. It is the solution of these 
equations, for specific choices of A(r), that we give in this section. 

A(t) = ci + C2T 

A simple form for A that allows direct integration of equations ()8.3.6|) and ()8.3.7p is 
the linear function A = ci + C2T. From equations ()8.3.3|) and ()8.3.4|) it can be seen 
that the constant ci can be absorbed into ri and T2 by simple redefinitions. The 
equations ()8.3.6|) and ()8.3.7p then become 




^ _ (4 - 37 - 2c2)r + T2 

~ (ir2 + BT + C) 

_ ((2 - 7)cu + 1 + C2)r 

a ~ {At^ + Bt + C) 



where A = A — C2, B = B + C2T2 and C = D. The solutions of these equations 
depend upon the sign of the discriminant 



(8.4.1) 
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For the case A = 0, there exist simple exact power-law solutions 

2(2-7)c^+2+2e2 

a(r) OC T 3^^"-6^(l+2")+2(5+6i^)-2c2 (8.4.2) 

2(4-37)-4e2 

(/)(r) OC 7-3^^"-6^(l+2'^) + 2(5+6i^)-2c| _ (8.4.3) 

Substituting these power-law solutions into ()8.3.5|1 . we can obtain the corresponding 
power-law form 

4c2(2-37-C2)-67(l+(2-7)^) 
p ^ J- 372i^-67(l+2ij)+2{5+6ij)-2c2 

The relationship between r and the cosmological time t can now be obtained by 
integrating the definition dr = pa^dt given in eq. ()8.3.1|) . This gives ()8.4.2j) and 
fj8.4.3|) in terms of t time as 

2+2(2 — i)ui+2c2 

a(t) ~ t4+37c^{2-7)-2c2(7-67-C2) (8.4.4) 

2(4-37)-4e2 

(pit) ~ t4+37-(2-7)-2c2(7-67-C2). (8.4.5) 

The condition required for the occurrence of power-law inflation is obtained by 
requiring the power of time in equation ()8.4.5|) to exceed unity (for the case with 
out energy transfer see refs. |122t 1115] ). For oj > —3/2 we always have A ^ 0, and 
the case A > possesses the exact solutions 

B((2-7)t^ + l+C2) 



{2^7)^+i±£2 / 2ir + 5 + VA\ 
2At + B -V^ 



^ n V z,— / jLL/-r i-TLT I z.r-\i —t- ij —t- \/ I \ \ 

{At' + Bt + C) ^ — . ' (8-4-6) 



(4-37-2c2)-B-T2A 

. ^ . . 4-37-2C2 / 2At + B+ \/A\ 

0(r) = 0o(Ar^ + Bt + C)^a^ [ J ^ ^ ^ ^ \ . (8.4.7) 

where Oq and 0o are constants of integration. For uo < —3/2 we have A ^ 0, and 
the case A < has the exact solutions 

- r, - - (2-7)c^ + l + C2 

a{T) = ao{AT' + Bt + C) 2a 



X exp 



((2-7)cj + l + C2)5^ _Jb + 2At 

tan ' 



f8.4.8) 
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X exp 



2T2A - 2(4 - 37 - 2c2)5 . B + 2At 
— ; tan , 



(8.4.9) 



These solutions have the same functional form as those found by Gurevich, Finkel- 
stein and Ruban ^97, for Brans-Dicke theory, in the absence of energy exchange (A = 
constant), and reduce to them in the limit C2 0. The behaviour of these solutions 
at early and late times will be discussed in the next section. 



A(t) = Car-, n / 1 

We now consider forms of A(r) that are more general than a simple linear function 
of T. Making the choice A = Csr", where n ^ 1 and C3 is constant, and setting 
the free scalar component to zero (r2 = 0), we find that ()8.3.6|1 and ()8.3.7|) can be 
integrated exactly. The form of the solutions again depends upon the roots of the 
denominator. For real roots we require 00 > —3/2, for which we find the solutions 



a(r) = aor^^^^ f±2c3r"-' T (4 - 37) ± (2 - 7)^3(3 + 2cu) 



3+3(2-7)c^-V3(3+2a.) 
3K{n-l) 



X 



±2c3r"^' T (4 - 37) T (2 - 7)v/3(3 + 2u;) 



3+3(2-7)aj + Y^3(3+2uj) 
3K,(n—l) 



(8.4.10) 



2(4-37) 



±2c3r"-i T (4 - 37) ± (2 - 7)v/3(3 + 2a;) 



(4-37) + (2-7)V3(3+2a.) 
K(n — 1) 



±2c3r"-i T (4 - 37) T (2 - 7)^3(3 + 2c^) 



(4-37)-(2-7)V3(3+2a;) 
,,{n-l) 



(8.4.11) 



where k = 2(5 — 37) + 3(2 — 7)^0;. For a denominator with imaginary roots we 
require uj < —3/2, for which we find 



l+(2-7)a; 



a(r) = ao [±2(4 - 37)c3r^-" ± /cr^^^"") ^ 2c^] ''(^-^ 

(4 - 37) - 2c3t"-^ 



3fi;l-n 



(2-7)v/-3(3 + 2o;) 



(8.4.12) 
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0(r) = 00 [±2(4 - 37)c3r^-" ± /tr'^'^") 





tan 



/ (4-37)-2c3t- 



) 



(8.4.13) 



with K defined as above. The ± and =F signs here indicate that there are multiple 
solutions that satisfy the field equations. These signs should be chosen consistently 
within each set of square brackets (solutions ()8.4.10|) and ()8.4.1H) can have upper 
or lower branches chosen independently in each set of square brackets, as long as 
a consistent branch is taken within each separate set of square brackets). The 
physical branch should be chosen as the one for which the quantity in brackets 
remains positive as r ^ oo, so ensuring the existence of a positive real root in this 
limit, 

These solutions display interesting new behaviours at both early and late times, 
which will be discussed in the next section. 

8.5 Behaviour of Solutions 

A(t) = Ci + C2T 

At late times, as r — oo, the solutions ()8.4.6p - ()8.4.9|) all approach the exact power- 
law solutions ()8.4.2|1 and ()8.4.3|) . It can be seen that these solutions reduce to the 
usual spatially flat FRW Brans-Dicke power-law solutions |137j in the limit that the 
rate of energy transfer goes to zero, C2 ^ 0. It can also be seen that these solutions 
reduce to the spatially-fiat FRW general relativistic solutions in the limit ^ oo, 
irrespective of any (finite) amount of energy transfer. 

The early-time behaviour of these solutions approaches that of the general Brans- 
Dicke solutions [nZ|, without energy transfer. Generally, we expect an early period 
of free-scalar-field domination except in the case T2 = 0, in which case the power- 
law solutions ()8.4.2|) and ()8.4.3p are valid right up to the initial singularity. For 
LJ > —3/2, the scalar-field dominated phase causes an early period of power-law 
inflation. In this case there is always an initial singularity and the value of the 
scalar field diverges to infinity or zero as it is approached, depending on the sign 
of T2. For uj < —3/2 there is a 'bounce' and the scale factor has a minimum non- 
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zero value. In these universes there is a phase of contraction followed by a phase of 
expansion, with no singularity separating them. Solutions of this type were the focus 
of jSH where the evolution of through the bounce was used to model the variation 
of various physical constants in such situations. The energy exchange term does not 
play a significant role at early times in these models. The asymptotic solutions as 
the singularity (or bounce) is approached are the same as if the energy exchange 
term had been neglected, and are given by |144j , up to the absorption of Ci into T2 
previously described. 

A(t) = n ^ 1 

The behaviour of the solutions ()8.4.10|1 - ()8.4.13|) depends upon the signs of n — 1 
and C3, as well as on the sign of + 3/2. For illustrative purposes we will consider 
the radiation case 7 = 4/3 which is appropriate for realistic universes dominated by 
asymptotically-free interactions at early times. 

For > 1 it can be seen that the late-time attractors of solutions ()8.4.10|1 - 
fl8.4.13|l are a —>■ constant and (p ~^ constant as r ^ 00, for both uj > —3/2 and 
uj < —3/2. At late times these universes are asymptotically static; the evolution 
of the scale-factor ceases as r ^ 00 and both (p and p become constant. Further 
analysis is required to establish whether these static universes are stable or not (we 
expect them to be stable as no tuning of parameters or initial conditions has been 
performed to obtain these solutions). 

The early-time behaviour of solutions with n > 1 depends upon the sign of C3 
as well as whether u is greater or less than —3/2. We will consider first the case 
of cij > —3/2. For C3 > 0, we see that a — > 00 as r"^"*^ t^; for C3 < we see 
that a ^ 00 as r"^^ Tq (where r(+ = ((2 - -f)^3{3 + 2iu) + (4 - 37))/2c3 and 
r- = -((2 - 7)v/3(3 + 2cj)- (4 - 37))/2c3). For n > 1 and uj> -3/2 we therefore 
have the generic behaviour that a cxd at early times and a constant at late 
times. The behaviour of a at intermediate times varies in form depending on the 
sign of C3, as can be seen in figures |H3 and |H21 The asymptotic form of for n > 1 
and uj > —3/2 depends critically on the sign of C3. For C3 > it can be seen that 
^ as t"-^^ — > Tq*", whereas for C3 < it can be seen that ^ 00 as r"~^ Tq . 
The behaviour of in these two cases is illustrated in figures 18.11 and 18.21 
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Figure 8.1: The evolution of a and (j) for n = 2, uj = 10, C3 = 10 and 7 = 4/3. 



Figure 8.2: The evolution of a and for n = 2, u; = 10, C3 = —10 and 7 = 4/3. 

We now consider the early-time behaviour of solutions with n> 1 and uj < —3/2. 
It can be seen that a ^ as r ^ irrespective of the sign of C3, so that we find 
the generic behaviour a ^ at early times and a constant at late times (this 
is in contrast to the standard theory where an initial singularity is avoided when 
u < —3/2). Again, the behaviour of a at intermediate times is dependent on the 
sign of C3, as can be seen from figures and IH7^ As r — we see that has a 
finite non-zero value and is either increasing or decreasing depending on the sign of 
C3 . This behaviour is shown in figures 18.31 and 18.41 

It remains to investigate the nature of the solutions with n < 1. At late times 
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Figure 8.5: The evolution of a and 4> for n = 0, a; = 10, C3 = 10 and 7 = 4/3. 
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we see that 

2+2(2-7)t^ 
T-2{5-37) + 3(2-7)2a; 

2(4-37) 
7-2{5-37)+3{2-7)2uj ^ 

as r ^ oo, irrespective of the sign of C3 or the value of uo. These late-time attractors 
are the flat FRW power-law Brans-Dicke solutions |137j which reduce to the standard 
general relativistic solutions in the limit uj ^ 00. 

The early-time behaviour when n < 1 depends on the sign of C3 and the sign of 
^7 + 3/2. We consider first the case uj > —3/2. In this case it can be seen that a 
as r"^^ Tq or r"^^ Tq , for C3 > or C3 < respectively. The behaviour of a 
for both C3 > and C3 < is shown in figure 18.51 and 18.61 The behaviour of at 
early times depends on the sign of C3 and goes to 00 for C3 > or to for C3 < 0. 
The behaviour of in these cases is shown in figures IH7S1 and IH^ For uj < —3/2 the 
scale factor a contracts to a finite, but non-zero, minimum value and then expands. 
The exact form of the minimum depends on the values of n, u and C3, but it is 
interesting to note that odd values of n produce symmetric bounces and even values 
of n produce asymmetric bounces, as illustrated in figures IHTI and The evolution 
of (f) through these bounces is smooth with a time direction prescribed by the value of 
C3, as shown in figures IHTTI and IH7S1 (increasing for C3 > and decreasing for C3 < 0). 
The effect of changing the sign of C3 is seen to be a mirroring of the evolution of a 
and in the y— axis. 

8.6 Physical Consequences 

The solutions found in the previous sections are of physical interest for a number of 
reasons. The transfer of energy and momentum between a non-minimally coupled 
scalar field and matter fields is a prediction of a number of fundamental theories of 
current interest, including string theories, Kaluza-Klein theories and brane- worlds. 
The cosmologies produced by such an interaction, therefore, should be of direct 
interest in the consideration of these theories. Furthermore, the solutions we have 
found display modified behaviour at both early and late times. The investigation 
of modified theories of gravity at early times is of particular interest as it is in 
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the high-energy hmit that deviations from general relativity are usually expected. 
Modified behaviour at late times is also of interest as it is at these times that we 
can make direct observations which can be used to constrain deviations from the 
standard general relativistic model. We will now summarise the behaviour of the 
solutions found in the previous section, highlighting the physically significant results 
and constraints that can be placed on the theory by observations. 

For the case of A linear in r it was shown that the late-time attractors of the gen- 
eral solutions are no longer the power-law solutions of the Brans-Dicke theory, but 
are given by equations ()8.4.4|1 and ()8.4.5|1 . These attractors are of special interest 
as they have a simple power-law form that reduces to the general relativistic result 
in the limit u —* oo and to the Brans-Dicke result in the limit C2 ^ 0. Observations 
of cosmic microwave background anisotropies and the products of primordial nu- 
cleosynthesis will therefore be able to constraint any potential late-time deviations 
of this kind, and hence the underlying model. The process of primordial nucle- 
osynthesis in scalar-tensor theories has been used by a number of authors to place 
constraints on the coupling parameter a;(0) jHH Ell 11731 1172t 11671 EZl ED]- In these 
studies the different value of G during nucleosynthesis causes the weak interactions 
to freeze out at a different time and hence the proton to neutron ratio at this time is 
different to the standard case. This modification causes different abundances of the 
light elements to be produced, which can be compared with observations to constrain 
the underlying theory. Studies of this kind usually assume G to be constant during 
nucleosynthesis, which will not be the case when energy is allowed to be exchanged 
between and the matter fields. The effects of a non-constant G were studied in jSH ■ 
A similar study would be required to place constraints upon the parameters C2 and 
uj in this theory. The cosmic microwave background power spectrum has also often 
been used to constrain scalar-tensor theories of gravity |1 181 Il33t 11941 11341 1^ . In 
these studies the redshift of matter-radiation equality is different from its usual gen- 
eral relativistic value due to the modified late-time evolution of the Universe. This 
change in the redshift of equality is imprinted on the spectrum of perturbations as it 
is only after equality that sub-horizon scale perturbations are allowed to grow. The 
main effect is seen as a shift in the first peak of the power-spectrum, which can be 
compared with observations to constrain the theory. Again, the late-time evolution 
of the Universe is modified from the usual Brans-Dicke case by the energy exchange 
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that we consider, so that the previous constraints are not directly apphcable. 

For the case of a non-hnear power-law exchange of energy, described by A oc r", 
the late-time evolution of a and can be significantly modified. For n > 1, the 
solutions do not continue to expand eternally, but are attracted towards a static 
state where the time-evolutions of a, and p cease. For n < 1 the generic late- 
time attractor is the power-law solution of a flat FRW Brans-Dicke universe. It 
appears that theories of energy exchange with n > 1 are ruled out immediately by 
observations of an expanding universe whilst the case of n < 1 is subject to the same 
late-time constraints as the standard Brans-Dicke theory jHH ISS I173| I172| I167| EH 

EH EHi ESI cnsi iMi cnii E] • 

It remains to investigate the physical consequences of the early-time behaviour 
of our solutions. For A linear in r these solutions approach those of the standard 
Brans-Dicke theory as either the initial singularity or the minimum of the bounce 
are approached, according to the sign of — 3/2. The physical significance of this 
behaviour has been discussed many times before, usually focusing on the avoidance 
of the initial singularity and the inflation that can result from the presence of the 
free component of the scalar field. 

For A oc r" the early-time behaviour can be significantly changed from that of 
the standard theory. For n > 1 the scale factor a either approaches infinity or zero, 
depending on our choice of u and C3, as previously described. For the more realistic 
case of n < 1 the evolution of a at early times either undergoes a period of rapid 
expansion or a non-singular bounce, depending on whether uo is greater or less than 
—3/2. This behaviour is similar to that of the general solutions of the standard 
theory, but in this case the free scalar-field-dominated epoch has not been invoked 
and there is more freedom as to the exact form of the evolution. For example, with 
a suitable choice of parameters it is possible to create a universe that contracts and 
then is briefly static before 'bouncing' and continuing on to its late-time power-law 
evolution. This is shown in figure IHUl for the case uo = —10, n = —2 and C3 = 10. (It 
is interesting to note that Peter and Pinto- Neto remark that a static period followed 
by a bounce could potentially produce a scale-invariant spectrum of perturbations 
|152j ). For the physically reasonable models with n ^ 1 the evolution of (p, and hence 
of G, can be significantly altered at early times from what is generally assumed to be 
the case in scalar-tensor theories of gravity. For the case cu > —3/2, the value of G 
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Figure 8.9: The evolution of a and (j) for n = —2, uj = —10, C3 = 10 and 7 = 4/3. 

can be made to diverge to infinity or to zero as the initial singularity is approached, 
independent of whether or not there was an early scalar-dominated phase to the 
universe's history. For the case u < —3/2, the value of G evolves smoothly through 
the bounce in the scale-factor, and is again independent of whether or not there was 
a scalar-dominated phase. More complicated evolutions of can also be constructed, 
as can be seen in figure 18.91 

In performing this analysis we have allowed the coupling constant uj to lie in the 
ranges u > —3/2 and u < —3/2. It should be noted that when u < —3/2 the scalar 
field should be considered as having negative energy. Allowing an interaction with 
matter fields may then be considered as problematic as energy could potentially 
be transferred from the scalar field to the matter fields without bound, allowing 
runaway solutions. 

8.7 Discussion 

We have considered spatially flat FRW universes in scalar-tensor theories of gravity 
where energy is allowed to be exchanged between the Brans-Dicke scalar field that 
determines the strength of gravity and any perfect fluid matter fields in the space- 
time. We have presented a prescription for integrating the field equations exactly 
for some unknown function A which describes the rate at which energy is exchanged. 
For the case of A being a linear function of r we have found the general solutions 
to the problem and for the case of A being a non-linear power law function of r 
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we have been able to find a wide class of exact solutions. These solutions display 
behaviours that can deviate substantially from the corresponding solutions in the 
standard case, where the exchange of energy is absent. Depending upon the values 
of the parameters defining the theory and the exchange of energy, deviations in the 
evolution of a and (f) can occur at both early and late times, providing a richer 
phenomenology than is available in the standard theory. 

We have found that the parameter n must be bounded by the inequality n ^ 1 
if the Universe is to be expanding at late times. For n = 1 we have found late-time 
power-law attractor solutions which can be used to constrain the parameters C2 and 
uj. For n < 1 we have seen that the late-time evolution will be the same as in the 
standard Brans-Dicke case, and so is subject to the same observational constraints 
as these theories. The parameters C3 and T2 have been shown to be influential only in 
the vicinity of the initial singularity, or at the minimum of expansion in non-singular 
solutions. These parameters are therefore less accessible to constraint by late-time 
observations (see, however, [ZD] where the influence on primordial nucleosynthesis is 
used to constrain T2). The parameter oj is, as always, subject to the very tight solar 
system constraint uj > 40000 to 2cr [TFj . 

These results could be of interest in attempting to explain why G is so small in 
the present day Universe compared to the proton mass scale {Gm^^ ~ 10^^^). In 
these models the value of G can decay away by a coupling between the scalar field 
(p and the matter fields which allows energy to be transferred. The small value of 
G is then due to the age of the Universe. It remains to see whether or not the late- 
time modifications found above are consistent with observations of the primordial 
abundance of light elements, microwave background formation and other late-time 
physical processes. These studies should be able to be performed in an analogous 
way to the ones that already exist for the standard Brans-Dicke theory. 

Using the late-time solutions that have been found it is possible to comment on 
the case of FRW cosmologies with non-zero spatial curvature. At early times it is 
expected that the effect of any spatial curvature on the evolution of a{t) should be 
negligible. From the solution ()8.4.2|) we can see that spatial curvature will dominate 
the late time evolution if the condition 

2 -f- 2(2 - 7)cu + 2c2 ^ 
4 + 37o;(2 - 7) - 2c2(7 - 67 - ca) ^ 



146 



is satisfied. If tliis condition is not satisfied, then the power-law solution ()8.4.2j) will 
be an attractor as t ^ cxd even in the case of non-zero spatial curvature, offering 
a potential solution to the flatness problem. This behaviour corresponds to power- 
law inflation and it can be seen that the condition for equation ()8.4.2|) to dominate 
over the spatial curvature at late times is, indeed, also the condition that power-law 
inflation should occur. 

In conclusion, we have found that a direct coupling between (p and the matter 
fields in scalar-tensor cosmologies provides a richer frame-work within which one 
can consider variations of G. We have shown that it is possible to construct models 
where the late-time violations of the equivalence principle can be made arbitrarily 
small (for A oc r) or are attracted to zero (for A oc r" where n < 1). This en- 
larged phenomenology is of interest for the consideration of the four-dimensional 
cosmologies associated with higher-dimensional theories as well as for more general 
considerations of the variation of G and its late-time value. This study has been 
limited to scalar-tensor theories with constant coupling parameters, to fiat FRW 
cosmologies, and to special cases of that allow direct integration of the field 
equations. Obvious extensions exist in which these assumptions are partially or 
completely relaxed. 
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CHAPTER 9 



Cosmological Constraints 



Observations of physical processes occuring in an expanding universe can be used 
to constrain the underlying gravitational theory in a number of ways. These ob- 
servations can range from galaxy surveys in the nearby Universe to the results of 
processes occuring in the very early Universe. We will concentrate in this work on 
the processes of primordial nucleosynthesis and microwave background formation. 
Both of these processes occur early in the Universe's history and the results of both 
are well observed by astronomers and astrophysicists. 

This chapter is based on the work of Clifton, Barrow and Scherrer and 
Clifton and Barrow [THj . 

9.1 Scalar- Tensor Theories and Primordial Nu- 
cleosynthesis 

Using the scalar-tensor theories we investigate the earliest well understood physical 
process, primordial nucleosynthesis. Previous studies on this subject have been 
carried out by a number of authors. In particular, the Brans-Dicke theory has been 
especially well studied in this context by, for example, Casas, Garcia-Bellido and 
Quiros jni], |n21 and Serna, Dominguez-Tenreiro and Yepes |173j . The more general 
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class of Scalar- Tensor theories has also been well studied, most notably by Serna 
and Alimi jl72j . Santiago, Kalligas and Wagoner |167j and Damour and Pichon [77] . 
More recently, a number of studies have been performed investigating the effect of 
a nonminimally coupled quintessence field on primordial nucleosynthesis (see e.g. 
[Hnj. jl^Hj ). Although these studies are very detailed, they all make the simplifying 
assumption of a constant G during the radiation-dominated phase of the Universe's 
history (with the exception of jl73j . who numerically investigate the effect of an 
early scalar-dominated phase on the Brans-Dicke theory, and 77j who use the idea 
of a "kick" on the scalar field during electron-positron annihilation). We relax this 
assumption and investigate the effects of entering the radiation-dominated phase 
with a non-constant G. The constraints we impose upon the variation of G during 
primordial nucleosynthesis are then used to constrain the parameters of the theory. 
In carrying out this study we consider the more general class of scalar-tensor theories, 
paying particular attention to the Brans-Dicke theory. 

9.1.1 Modelling the Form of G{t) 

If primordial nucleosynthesis were to occur during the scalar-dominated period then 
the very different expansion rate would have disastrous consequences for the light- 
element abundances (see e.g. |173j ). Therefore we limit our study to times at which 
the scale factor can be approximated by a form that is close to a{t) oc t^, and 
so primordial nucleosynthesis can safely be described as occurring during radiation 
domination. Performing a power-series expansion of the FRW Brans-Dicke solutions 
(ra, and mM in Vi/iv + V2) we find 

a{r]) = ai{r] + r]2 + 3r]i) + 0{7]l) (9.1.1) 

m = <l^i(^-^)+0{vl) (9.1.2) 

V V + V2J 

for both uj > —3/2 and uj < —3/2. We can then set the origin of the 77 coordinate 
such that a(0) = with the choice ri2 = —3rii. The solutions ()9.1.1|) and ()9.1.2|) 
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then become, in terms of the time coordinate t, 



a(t) = ait^ + 0(?7i) 

where 02 = —Srjiai and the origin of t has been chosen to coincide with the origin 

of T]. 

Similarly, expanding the FRW u^cp) solutions ()3.3.6|1 and ()3.3.7|) in rji/^r] + 772) 
we find that, for both u > —3/2 and a; < —3/2, 

a{t) = ait^ + 0{r]l) 



where 



ai = ^/8Tipro/3 exp{p{^i - ipoof/2) 

02 = \/2proPr]i{tpi - tpoo) exp(/?(V'i - V^oo)V2) 



5)1 = exp(-/3(?/'i - V'oo)^) 



and ri2 has been set so that a(0) = 0. 

In the limit — > we can see from equations ()3.1.2p . ()3.1.3p and ()3.1.5p that the 
standard GR Friedmann equations are recovered, with a different value of Newton's 
constant given by 

(p{t) a[t) + 02 

where Gi = I /(pi- We conclude that the solutions found above correspond to a 
situation that can be described using the GR Friedmann equations with a different, 
and adiabatically changing, value of G. This is just the situation considered by 
Bambi, Giannotti and Villante |3] , but we now have an explicit form for the evolution 
of G{t) derived from scalar-tensor gravity theory. 
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(a) 2a bounds for 02 > 




(b) 2a bounds for a2 < 



Figure 9.1: The upper and lower bounds on Gi/Gn as a function 0/02 are shown by 
the solid lines. The shaded regions correspond to the allowed parameter space and 
the dashed lines show Gi = G^. 
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9.1.2 The Effect of a Non-Constant G 



The time at which weak interactions freeze out in the early Universe is determined 
by equahty between the rate of the relevant weak interactions and the Hubble rate. 
When the weak interaction rate is the greater then the ratio of neutrons to protons 
tracks its equilibrium value, n/p = exp(— (m„ — mp)/T) where m„ and nip are 
the neutron and proton masses. If, however, the Hubble rate is greater than the 
weak-interaction rate then the ratio of neutrons to protons is effectively 'frozen- 
in', and 13 decay is the only weak process that still operates with any efficiency. 
This will be the case until the onset of deuterium formation, at which time the 
neutrons become bound and /5-decay ceases. The onset of deuterium formation is 
primarily determined by the photon to baryon ratio, r]^, which inhibits the formation 
of deuterium nuclei until the critical temperature for photodissociation is past. As 
the vast majority of neutrons end up in ^He the primordial abundance of this element 
is influenced most significantly by the number of neutrons at the onset of deuterium 
formation, which is most sensitive to the temperature of weak-interaction freeze- 
out, and hence the Hubble rate, and so G, at this time. Conversely, the primordial 
abundances of the other light elements are most sensitive to the temperature of 
deuterium formation, and hence rj^, when nuclear reactions occur and the light 
elements form. (See jSj for a more detailed discussion of these points). 

Using the simple forms of a{t) and G{t) derived above we use a modified version 
of the Kawano code | 112j to investigate the effect of this variation of G on primordial 
nucleosynthesis directly. We use the deuterium abundance estimated by Kirkman 
et al. [TTI] 



to create a parameter space in (77^, Gi, 02). The three-dimensional 95% confidence 
region is projected into the Gi, 02 plane to give figure IHUl where three species of 
light neutrinos and a neutron mean lifetime of r = 885.7 seconds have been assumed. 
The apparent disfavoring of the value of Gi much greater than G^v in figure 




to Icr, 



(9.1.4) 



and the "^He abundance estimated by Barger et al. [Zj 



Yp = 0.238 ± 0.005 to la. 



(9.1.5) 
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19.11 is due to the observational limits we have adopted for Yp (equation ()9.1.5j) ). 
which suggest a helium abundance that is already uncomfortably low compared to 
the theoretical prediction for Yp derived from the baryon density corresponding to 
equation ()9.1.4|) . We expect these constraints to be updated by future observations 
(see e.g. |il85i| ): this will require a corresponding update of any work, such as this, 
that seeks to use these constraints to impose limits on physical processes occuring 
during primordial nucleosynthesis. 

These limits on the parameters Gi and 02 can be used to construct plots showing 
the explicit evolution of G{t) for various limiting combinations of the two parameters, 
this is done in figure IHSl It is interesting to note that in both of these plots the lines 
corresponding to different values of 02 all appear to cross at approximately the same 
point, log a ~ —9.4. This confirms our earlier discussion, and the results of that 
the ^He abundance is mostly only sensitive to the value of G at the time when the 
weak interactions freeze out. In reality, this freeze-out happens over a finite time 
interval, but from figure 19.21 we see that it is a good approximation to consider it 
happening instantaneously - where the lines cross. To a reasonable accuracy one 
could then take G throughout primordial nucleosynthesis to be its value during the 
freeze-out process. 

9.1.3 Constraining the Theory 

Using the results in the previous section it is possible to constrain the underlying 
scalar-tensor gravitational theory. This is done separately for the Brans-Dicke theory 
and for the more general scalar-tensor theories. For each theory we consider the case 
of a universe containing matter and radiation only and then a universe containing 
matter, radiation and a non-zero vacuum energy, with equation of state p = —p. 

Brans-Dicke theory 

Universe containing matter and radiation 

Using the Brans-Dicke matter-dominated FRW solutions ()H.2.1()|1 we can write the 
ratio of G at matter-radiation equahty, Geq, to its present-day laboratory value, Gjy, 
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Figure 9.2: These plots show the explicit evolution of G for some limiting values 
of Gi, as taken from figure W~l[ for various values of a2. Solid lines correspond to 
\a2\ = 10"^^, dashed lines correspond to \a2\ = 10^^^-^^^ and dotted lines correspond 
to laal = 10"^^ 
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as 



Geq {2u + 3) 0^ _ {2u + 3) / ao ^ '/^"""'^ _ {2u + 3) , ^ 



Gat (2^ + 4)0e, (2a; + 4) \aej " (2cu + 4) ^ ^ ^ ' ' ^ 

where we have used the expression for G in the weak- field hmit |191j : 

for Gat, and G(0) = 1/0 for Ggg, as in ()9.1.3|) . In reahty, the laboratory value of G 
today is not equal to that in the background cosmology [ZH] , although we take it to 
be so here for simplicity as this effect should be small. 

We now proceed by calculating 1 + Z(.q in the Brans-Dicke cosmology, following 
Liddle, Mazumdar and Barrow [118]. As T"''.^ = in the Jordan frame we have that 
Pr oc and pm oc a~^. For a universe containing matter and radiation only, this 
gives 

1 + Zeq= = = 

(^eq PrO Pmeq PrO 

as well as the usual relation T oc a~^, when entropy increase is neglected. 

Photons and neutrinos both contribute to the value of pro, the present-day 
energy-density of radiation. From T^q = 2.728 ± 0.0047^ [85J, we get p^o = 4.66 x 
10~^^g cm~^ and using the well known result T^q = (4/ll)^/^XyO) aiid assuming three 
families of light neutrinos, we also have pi^o = 0.68p^o- This gives the total present- 
day radiation density as pro = 7.84 x 10~^^g cm'^. Now, recalling our assumption 
of spatial flatness and ()3.1.3|) . we can write 

(4 + 3u;)(4 + 2c<j) 

PtotO — PmO + PrO — -^—^ ^7^— 7^ . 

87rGiv 6(1 -I- o;)^ 

For Gn = 6.673 x 10"" Nm^kg-^, Ho = lOO/i kms-^Mpc-\ and the value of pro 
above, we have 

4 o(4 + 3cj)(4-f2o;) 
l + z,q = 2.39 X lO^h^ ^ ^ - 1 

6(1 + u)^ 

-2\ 



2.39 X lO^/i^ (^1 + + 0{u 
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This correction to 1 + z^q has direct observational consequences in the power spec- 
trum of cosmic microwave background perturbations. After 1 + Zeq the subhorizon 
scale perturbations, that were previously effectively frozen, are allowed to grow. 
Changing the value of 1 + Zeq therefore causes a shift in the power-spectrum peaks. 



which is potentially observable (see |118j . [HK] . |194j . |lH4j and [2] for a more 



detailed discussion of the effect of a varying G on microwave background formation). 
For our purposes this modified expression for 1 + Z(,q CcLll then be substituted into 
fj9.1.6|) to give an equation in terms of cu, G^q and h. Assuming a value oi h = 0.7 
and that Geq = Gi, as appears a very good approximation for the models above, we 
can use our bounds on Gi in terms of 02 to create an allowed parameter space in 
the (co',a2) plane. This is shown in figure 

We remind the reader that a large uj corresponds to a slowly varying Machian 
component of 0(1 + z), as can be seen directly from equations ()3.2.10p . ()3.2.11|) 
and the late-time limits of (prTK|l . (jTTTjl . (prTK|l and (jTT^ . The parameter 02 
determines the evolution of the free component of 0, as can be seen from the early 
time limits of (ITTHIl . (ITT7I) . (ITT^ and ^^TT^ . In the limits w ^ 00 and 02 ^ 
the Machian and free components of 0(1 + z) both become constant, respectively, 
resulting in a constant G{1 + z). The constraints imposed upon uo and 02 in figure 
19. 31 therefore correspond to constraints upon the evolution of G{l + z) in this theory, 
valid both during primordial nucleosynthesis and at other cosmological epochs. 

Universe containing matter, radiation and a nonzero vacuum energy 

A more realistic constraint would involve taking into account a late-time period 
of vacuum domination; so as well as pr oc and pm oc a~^, we now also have 
Pa = constant. Hence, 





{2UJ + 3) 00 _ {2UJ + 3) 0eg2 00 
{2U0 + 4) 0egl ~ (2CJ + 4) (peql 0eg2 
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where we have defined the redshift of matter-radiation equahty, Zeqi, and the redshift 
of matter- vacuum equahty, Zeq2, as 

1 + Zeql = and 1 + Zeq2 = 

PrO \ PmO / 

As above, we still have pro = 7.84 x lO"^'^^' cm^^ but now our assumption of spatial 
flatness gives 

(4 + 3cj)(4 + 2o;) 

PtotO — PAO + PmO + Pro — F71~, \2 ' 

or, using Gn = 6.673 x 10"^^ Nm'^kg"'^ and = lOOh kms~^Mpc~^, 

(2.39 X 104/^2(±i|£i^ _ 1' 



1 + z, 



eql 



I _|_ PAO 

PmO 



Taking the value pAo/pmo = 2.7, consistent with WMAP observations [IS], we get 
the results shown in figure 

Dynamically-coupled theories 

Constraints on uj{(f)) at matter-radiation equality 

Recalling that e^^ = 1/0 and = 2c<j + 3 allows us to re-write ()9.1.7|) as 

G = e2r(l + «2), 
the ratio Geq/GN can then be expressed in terms of F and a as 
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-500L 



Figure 9.3: The allowed parameter space in this plot is the region above the higher 
solid line and below the lower solid line, with h = 0.7 to 2(J for a universe containing 
matter and radiation only. The dashed lines show the corresponding result for a 
universe with a non-zero vacuum energy. 

Assuming that Fq and ln(l + ao) are negligible compared to Fgq (i.e. the Universe 
is close to GR today, jUj) we can write 

ln(^^ ^F,, = l/3(V',,-^^)2 
_ 27r 1 

This allows us to constrain Ueq in terms of f3 and 02, as shown in figure Con- 
straints imposed upon u and 02 have the same consequences for the evolution of G 
as previously discussed. The parameter P controls the evolution of u, as can be seen 
from ()3.3.5|) . In the limit /3 — it can be seen that u becomes constant and this 
class of scalar-tensor theories becomes indistinguishable from the Brans-Dicke the- 
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Figure 9.4: The allowed parameter space in this graph is the region above the line 
in the region > and below the line in the region uj < Q, to 2a. The solid lines 
correspond to 02 = 10~^^ the dashed lines to 02 = 10~^^ and the dotted lines to 
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UJ. 



eg 



is the value of u at matter-radiation equality and (5 is a parameter 



of the theory, defined in \3. 3. 



ory. Constraints upon (3 therefore correspond to constraints on the allowed variation 
of uo and hence G. 

The parameter (5 is taken to be small here so that the "kick" on the scalar 
field during electron-positron annihilation can be neglected. These effects have been 
explored by Damour and Pichon in jTTj, for the case 02 = 0, and are expected to 
have the same result in this more general scenario; we will not repeat their analysis 
of this effect here. 

Constraints on ujq for a universe containing matter and radiation 

The scalar field can be evolved forward in time from the time of matter-radiation 
equality to the present, using the solution ()3.3.9p . The two arbitrary constants in 
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this expression can be fixed using our limiting values of Geq derived above and by 
assuming that the evolution of the scalar field has effectively ceased by this time, 
i.e. V'eg = 0) is the case for the models considered here. 

In order to gain a quantitative limit on ipi —ipooi and hence on cuq, it is necessary 
to calculate iVeg. From the definition of N , we can write 

Neq = - ln(l + Zeq) - Teq + Tq 

~ - ln(l + Zeq) - ^PiAq " V'oo)^ 

where Zeq is the redshift at teq and we have assumed, as before, the term In Aq to be 
negligible. 

It now remains to determine 1 + Zeq = Pmo/Pro for the case uj = uj{(j)). If we now 
assume ujq to be moderately large, and recall our assumption of spatial flatness, we 
can write 

PtotO — PmO + PrO 



For Pro = 7.84 x lO'^'^g cm~^ this gives pmo — 1-87 x 10~^^ h'^g cm~^. We have 

^eq 

1 

Evolving ip — ipac from Neg to A^^q = and using 



1 + Zeq ^ 2.4 X 10"^ h'^ and so finally obtain, for h = 0.7, 



Neg- -9.37 --PitPeq-i^oof. 



27r 3 , , 

we obtain the results shown in figure 19.51 

Constraints on ujq for a universe containing matter, radiation and nonzero 
vacuum energy 

We can repeat the previous analysis for the more realistic case of a universe with a 
period of late-time vacuum domination. We now need the time of matter-radiation 
equality, Negi, and the time of matter- vacuum equality, Neg2- 

Assuming spatial flatness, a large Uq, h = 0.7, pAo/pmo = 2.7, and pro as above 
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Figure 9.5: The allowed parameter space in this graph is the region above the line 
in the region u > and below the line in the region u < 0, to 2a for a universe 
containing matter and radiation only. The solid lines correspond to 02 = lO""*^^, the 
dashed lines to 02 = — lO^^"'^'^^^, and the dot-dashed lines to 02 = — 10~^^'^^^. The 
dotted lines show the corresponding results for a universe with a non-zero vacuum 
energy. 

gives 

Negi - -8.06 - ^PiAgi - i^oof, 

and 

N,,2 ^ -0.33 - -(3{^eg2 ' V^oo)'- 

We can now evolve ipi^N) — ipoo from N^qi to Neq2 using ()3.3.9|) and the same 
boundary conditions as before. (In finding Nf.q2 we used an iterative method to 
evaluate ipeq2 — i'oo)- We then use the solution ()3.3.10j) to evolve the field from Neq2 
to No] the constants C and D in ()3.3.1()|1 are set by matching ipi^N) — ip^o and its 
first derivative with the solution ()3.3.9p at Nf.q2- Finally, we can calculate ujq using 



161 



()9.1.8|1 . The results of this procedure are shown in figure 1^31 



9.1.4 Discussion 

Using the framework provided by scalar-tensor theories of gravity we have investi- 
gated the effect of a time-varying G during primordial nucleosynthesis. We deter- 
mined the effect on primordial nucleosynthesis numerically, using a modified version 
of the Kawano code |112j . and constrained the parameters of the underlying the- 
ory using these results. Our results are consistent with the interpretation that the 
abundance of ^He is primarily only sensitive to the value of G at the time when 
weak interactions freeze out, for the class of scalar-tensor models studied. 

Using our numerically determined constraints on the evolution of G, we imposed 
the 2(7 upper and lower bounds shown in figure 19.31 on the Brans-Dicke parameter 
oj. For a constant G (i.e. 02 = 0) we get the bounds uj > 332 or < —37, for a 
universe containing matter and radiation only, and the bounds uj > 277 or u; < —31, 
for a universe containing matter, radiation and a nonzero vacuum energy. As the 
parameter 02 is increased, the strength of the upper bound decreases whilst that of 
the lower bound increases; the opposite behaviour occurs as 02 is decreased. This 
is because the strength of the bound on uj is essentially due to the allowed value of 
G at matter-radiation equality. If primordial nucleosynthesis allows this value to be 
vastly different from Gat then a significant evolution of G during matter domination, 
and hence a low \uj\, is permitted. A value of Geq close to G^ means that only a 
very slow variation of G, and hence large |cc;| is permitted. As 02 > corresponds 
to an increasing G during primordial nucleosynthesis this corresponds to a higher 
value of Geq and hence a looser upper bound on oj and a tighter lower bound (as 
G decreases during matter domination for uj > —3/2 and increases for u < —3/2); 
02 < corresponds to G decreasing during radiation domination, and so has the 
opposite effect. 

The more stringent effect of a nonzero value of 02 on the upper bounds is due 
to the current observational determinations of the ^He abundance [7j disfavouring 
G > Gn during primordial nucleosynthesis (see e.g. pS]). 

The interpretation of the constraints on the more general ST theories is a little 
more complicated, due to the increased complexity of the theories. The bounds on 
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ci;(0) at matter-radiation equality, shown in figure \9A\ are seen to be stronger for 
smaller /3 and weaker for larger for f3 (assuming (3 to be small enough to safely ignore 
the effect of the e~e^ kick analysed by [HI). This should be expected as ~ in 
the models we are studying. We see from the constraints on u at the present day, 
shown in figure 19. 5| that the bounds on uj become tighter as (3 gets very small and 
large, with an apparent minimum in the bound at /5 ~ 0.2. The tight bounds for 
very small j3 are due to the tight bounds on uJeq for small (3. The tight bounds at 
large j3 are due to the attraction towards GR at late times that occurs for this class 
of scalar-tensor theories (see e.g. [ZE]). This attractor mechanism is more efficient 
for larger values of /5, as can be seen from the solution ()3.3.4j) . and at late times so 
u is drawn to a larger value for a larger f3. 

The effect of including a late-time vacuum-dominated stage of the Universe's 
evolution is to weaken slightly the bounds that can be placed on u at the present 
day, as can be seen from figures 19.31 and 19.51 This weakening of the bounds is 
due to a shortening of the matter-dominated period of the Universe's history which 
is essentially the only period, after the effects of the free scalar-dominated phase 
become negligible, during which G evolves. 

We find that the constraints that can be imposed upon the present day value 
of u from primordial nucleosynthesis are, for most of the allowed parameter space, 
considerably weaker than those obtained from observations within the solar system. 
To date, the tightest constraint upon uq are imposed by Bertotti, less and Tortora 
(TF] who find \ujq\ > 40000, to 2a. This constraint is obtained from observations of 
the Shapiro delay of radio signals from the Cassini spacecraft as it passes behind 
the Sun. We consider the constraints imposed upon u here to be complementary 
to these results as they probe different length and time scales, as well as different 
epochs of the Universe's history. 

9.2 Fourth-Order Theories 

The modified cosmological dynamics for fourth-order gravity lead to different pre- 
dictions for the outcomes of primordial nucleosynthesis and microwave background 
formation, compared to the standard general relativistic model. The relevant mod- 
ifications to these physical processes, and the bounds that they can impose upon 
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the theory, will be discussed in this section. We will use the solutions ()4.1.5|1 as 
they have been shown to be the generic attractors as t ^ oo (except for the case 
— J < 6 < when 7 = 1, which has been excluded as physically unrealistic on the 
grounds of structure formation) . 

9.2.1 Primordial Nucleosythesis 

We find that the temperature-time adiabat during radiation domination for the 
solution ()4.1.5p is given by the exact relation 

= A (9.2.1) 



where, as usual (with units h 



where g is the total number of relativistic spin states at temperature T. The con- 
stant A can be determined from the generalised Friedmann equation ()4.1.3p and is 
dependent on the present day value of the Ricci scalar, through equation 1)2.1.211) . 
(This dependence is analogous to the dependence of scalar-tensor theories on the 
evolution of the non- minimally coupled scalar, as may be expected from the rela- 
tionship between these theories |124j ). As a first approximation, we assume the 
universe to have been matter dominated throughout its later history; this allows us 
to write 

_ / 45(l-2^)(l-2^-5^^) \ (2{1 + S)\ 
""-y 32(l-%vr3G )[ 3Ho J ^'•'•'^ 
where Hq is the value of Bubble's constant today and we have used the solution 
f|4.1.5|) to model the evolution of a{t). Adding a recent period of accelerated expan- 
sion will refine the constant A, but in the interests of brevity we exclude this from 
the current analysis. As usual, the weak-interaction time is given by 

1 
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The freeze-out temperature, T/-, for neutron-proton kinetic equilibrium is then de- 
fined by 

t{Tf)=t^k{Tf). 

Hence the freeze-out temperature in this theory, with S ^ 0, is related to that in the 
general relativistic case with 5 = 0, T^^: 



Tf = C{Tp)T^ (9.2.3) 
where 



.5(1+^) 



^ ( (1 - ^) \ '"'^''^ 45 \ ^(3+5*) / 3^0 \ (3+5.) 

The neutron-proton ratio, n/p, is now determined at temperature T when the equi- 
librium holds by 

n ( Am 



— — exp 



P \ T 

where Am is the neutron-proton mass difference. Hence the neutron-proton ratio 
at freeze-out in the R^^^ early universe is given by 



where 




'3 + 55 

The frozen-out n/p ratio in the R^'^^ theory is given by a power of its value in 
the general relativistic case, ^ 1/7, by 



GR 

p \pji 



n In 



GR 

■^GR\e 



It is seen that when C{T^^Y > 1 {5 < 0) there is a smaller frozen-out neutron- 
proton ratio that in the general relativistic case and consequently a lower final 
helium-4 abundance than in the standard general-relativistic early universe contain- 
ing the same number of relativistic spin states. This happens because the freeze-out 
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temperature is lower than in GR. The neutrons remain in equihbrium to a lower 
temperature and their slightly higher mass shifts the number balance more towards 
the protons the longer they are in equilibrium. Note that a reduction in the helium-4 
abundance compared to the standard model of GR is both astrophysically interest- 
ing and difficult to achieve (all other variants like extra particle species |178[ 1180] . 
anisotropies [ESI El El IH, magnetic fields |2H1 ED], gravitational waves El EI], 
or varying G [HHl 1^ 17^ . lead to an increase in the expansion rate and in the final 
helium-4 abundance). Conversely, when C{Tj^^y < 1 (5 > 0) freeze-out occurs at a 
higher temperature than in general relativity and a higher final helium-4 abundance 
fraction results. The final helium-4 mass fraction Y is well approximated by 

Y = (9.2.5) 

(1 + n/pj 

It is now possible to constrain the value of 6 using observational abundances of 
the light elements. In doing this we will use the results of Carroll and Kaplinghat 
jnni who consider nucleosynthesis with a Hubble constant parametrised by 

Our theory can be cast into this form by substituting 

2 

a 



and 

(1+6) 1 2 

Hi = ^ ^ ' A~^(^{lMeV)(^, 

so, taking g = 43/8, G = 6.72 x lO^^^MeV-^ and Hq = 1.51 x IQ-^^MeV 03], this 
can be rewritten as 

(1 + 5) f 7.96 X lO-^^n _ ^) \ 2{iT^ / 2.23 X 10-39 _ 
Carroll and Kaplinghat use the observational abundances inferred by Olive et. al. 
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0.228 ^ Fp ^ 0.248 

2 ^ lO^x— ^ 5 
H 

1 ^ 10^°x^ ^ 3, 
H 



to impose the constraint 



,MeV 

where = 2.6 ± 0.9 x 10~^^MeV at = 0.2MeV for 0.5 ^ r^io ^ 50, or i/^ = 
2.0 ± 0.3 X 10~^^ for 1 ^ rjio ^ 10 and r/io is 10^° times the baryon to photon ratio. 
These resuhs can now be used to impose upon 6 the constraints 



for 0.5 ^ r/io ^ 50, or 



for 1 ^ r/io ^ 10. 



-0.017 ^ (5 ^ 0.0012, 

-0.0064 ^ 5 ^ 0.0012, (9.2.6) 



9.2.2 Microwave Background Formation 

The horizon size at the epoch of matter-radiation equahty is of great observa- 
tional significance. During radiation domination cosmological perturbations on sub- 
horizon scales are effectively frozen. Once matter domination commences, however, 
perturbations on all scales are allowed to grow and structure formation begins. The 
horizon size at matter-radiation equality is therefore frozen into the power spectrum 
of perturbations and is observable. Calculation of the horizon sizes in this theory 
proceeds in a similar way to that in Brans-Dicke theory |118j . 

In making an estimate of the horizon size in theory we will use the gener- 
alised Friedmann equation, ()4.1.3|) . in the form 

^ +^^R- 6(rT5) - 3(1-25) B^P- ^^-^-^^ 
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Again, we assume the form ()4.1.5|1 to model the evolution of the scale factor during 
the epoch of matter domination. This gives 



«(^) = cto — 

rr 2(1+5) 

-no = 



2(1+^) 
t \ ^ 



pm 

m 



3to 

3H^ {1 - 26){2 - 36 - 86^) 4 
16nG (l-(5)(l + 5)2 ^ 
4(1 + 55 + 4^2) 



3t2 

during the matter-dominated era. In order to simplify matters, we assume the 
above solutions to hold exactly from the time of matter-radiation equality up until 
the present (neglecting the small residual radiation effects and the late time accel- 
eration). Substituting them into ()9.2.7j) . along with p^q = 2pmeq at equality, we can 
then solve for Hgq to first order in 6 to find 



aoHi 







1-26 



2^1 + Zeq (1 - 2.6865) + 0(5') (9.2.8) 



where z^q is the redshift at matter radiation equality and H has been treated as an 
independent parameter. The value of 1 + z^q can now be calculated in this theory 
as 

1 + Ze, = (9.2.9) 

PmO 

Taking the present day temperature of the microwave background as T = 2.728 ± 
0.004/^ IHSl gives 

p^o = 3.37 X lO^^^Mel/^ (9.2.10) 

where three families of light neutrinos have been assumed at a temperature lower 
than that of the microwave background by a factor (4/11)^. Using the same values 
for G and ifo as above we than find from the above expression for that 

= 2.03 X 10-»il:i|^l|^l^i^Afel/^ (9.2.11) 

Substituting (j9.2.1();i . (j9.2.njl and (^TT^ into ^^T^ then gives the expression for 



168 



the horizon size at equahty, to first order in 6, as 

~ 155(1 - 195) + 0{6^). (9.2.12) 

This expression shows that the horizon size at matter-radiation equahty wiU be 
shifted by ~ 1% for a value of 5 ~ 0.0005. This shift in horizon size should be 
observable in a shift of the peak of the power spectrum of perturbations, compared 
to its position in the standard general relativistic cosmology. Microwave background 
observations, therefore, allow a potentially tight bound to be derived on the value 
of S. This effect is analogous to the shift of power-spectrum peaks in Brans-Dicke 
theory (see e.g. |118j . [HK]). 

A full analysis of the spectrum of perturbations in this theory requires a knowl- 
edge of the evolution of linearised perturbations as well as a marginalisation over 
other parameters which can mimic this effect (e.g. baryon density). Such a study is 
beyond the scope of the present work. 
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CHAPTER 10 



Weak-Field Constraints 



It is the weak-field limit of gravity that is most readily accessible to us for experiment 
and observation. By performing gravitational experiments in the solar system, and 
comparing the results of these experiments to the predictions from different gravita- 
tional theories, it is possible to place constraints on the theory under consideration. 
To date, it is weak field gravitational experiments that have so far offered the best 
constraints on the theory (though with astrophysical and cosmological observations 
becoming ever more accurate there is a possibility of the situation changing in the 
not too distant future). 

The consideration of gravitational experiments in the weak-field, low veloc- 
ity limit is usually performed within the frame-work of the parameterised post- 
Newtonian (PPN) formalism. In this chapter we will give a brief explanation of 
the PPN formalism, and explain how it can be used to constrain scalar-tensor and 
fourth-order theories of gravity. We will also point out the areas in which this 
formalism cannot be applied so easily, and provide alternative ways to proceed. 

This chapter is based on the work of Clifton and Barrow [ZOI and Clifton jHSj- 
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10.1 Parameterised Post-Newtonian Approach 



We will provide in this section a brief explanation of the PPN approach (for a more 
detailed explanation of this formalism the reader is referred to Will, |191j ). 

The PPN formalism does not deal with exact solutions, but instead uses pertur- 
bative expansions about a Minkowski background. In the weak-field, low velocity 
limit we need only deal with the lowest orders of the perturbative expansion. Solv- 
ing the field equations order by order, we can then find approximate solutions which 
can be compared with observation. 

Following Will we define an "order of smallness" such that 

Ur^V^r^-r^Ur^ 0(2) 
P 

where U is the Newtonian potential, v is the velocity of a body in the space-time, 
p is pressure, p is energy density and 11 is specific energy density (ratio of energy 
density to rest mass density) and 0(2) means second order in "smallness". In the 
solar system we have that U ^ 10^^ is a small perturbation to flatness, and all of 
the other above quantities are less that U and so are also small (see Will |191j for 
more details). It is also assumed that the time evolution of the solar system is slow 

^^0(1). 

\d/dx\ ^ ' 

A perturbative expansion about the Minkowski background can now be performed 
in the order of smallness that has just been defined. 

To obtain the Newtonian limit for particles following time-like geodesies we now 
need to know only the 0(2) term in the goo component of the metric. The Newtonian 
limit of light rays are straight lines, so for null geodesies the Newtonian limit is given 
simply by the Minkowski background. The post-Newtonian limit is then given by 
the next non-zero order in the perturbative expansion of the metric components. 
For time- like particles this is given by goo to 0(4), go^ to 0(3) and gf^i, to 0(2). For 
null particles it is given by goo to 0(2) and gf^,^ to 0(2). 

The usual approach at this point is to construct a metric containing a number 
of the possible functionals of the matter variables, to the required post-Newtonian 
order. Giving the different functionals (often referred to as potentials) arbitrary 
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coefficients it is then possible to construct a generalised metric. The values of the 
coefficients are then given by solving the field equations of the theory being con- 
sidered, order by order in perturbation. To restrict the number of potentials in the 
metric a number of conditions are usually imposed such as asymptotic fiatness and 
simplicity of the functional form. After specifying a particular gauge the generalised 
metric can then be written down explicitly. 

The appeal of this approach is that the generalised metric can be used by experi- 
menters to place constraints on the coefficients of the metric potentials, independent 
of any further considerations of the specific gravitational theory being tested. For 
any particular gravitational theory (that fits into this framework) the values of the 
coefficients can then be expressed in terms of the parameters of the theory. In this 
way the results of gravitational experiments can be applied to any gravitational 
theory that fits into the PPN frame-work. 

A simplified version of the post-Newtonian metric is given by 

^700 = -1 + 2f/ - 2PU' + (2 + 27 + Ci)$i (10.1.1) 

+ 2(1 - 2/5 + 37 + (2)^2 + 2(1 + C3)'^'3 + 2(37 + 3C4)$4 
g^, = (1 + 27f/)(5^, 



90^ = -^(3 + 47 + Ci + ai - a2)V^ - ^(1 - Ci + (y2)W^) 



where P, 7, (i, (2, Cs; c^i and 02 are some of the possible post-Newtonian coef- 
ficients and $1, $2) "^^3) "^^4, Vfi and are post-Newtonian potentials. (See |191j 
for a more comprehensive version of the post-Newtonian metric, and definitions of 
the above post-Newtonian potentials). For an idealised spherical and non- rotating 
massive object at the origin, in an otherwise empty space-time, this metric reduces 
to 

ds^ = - (1 + f3 — ]dt^ + {1 + 7 daudx^dx" 

\ r J \ r J 

which is exactly the metric considered by Eddington [SSI > Robertson |161j and Schiff 

mni. 

The PPN approach has been shown to work very well for a number of grav- 
itational theories, but it does have some drawbacks. From the point of view of 
the present study the principle drawback is the assumption of the Minkowski back- 
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ground. Whilst the Minkowski background is very well justified in GR, it is less well 
justified in some alternative theories. This is principally due to the lack of any equiv- 
alent to Birkhoff's theorem. In fact, in a previous chapter we have shown explicitly 
the existence of solutions which are spherically symmetric, non-static vacuum so- 
lutions. This breakdown of Birkhoff's theorem can be understood quite simply in 
terms of the conformal equivalence of these theories in vacuum to GR in the pres- 
ence of a scalar field. Just as it conceivable that a scalar field in GR could cause 
a cosmological evolution (even in the absence of any other matter sources), so it is 
conceivable here that the vacuum solutions of these modified theories could evolve 
cosmologically. 

These problems are especially evident in the ET' theories. In these theories, 
even in the static case, it has already been shown that at asymptotically large 
spatial distances the general solution is not attracted towards fiatness. This strange 
asymptotic behaviour means that the theories do not fit straightforwardly into the 
PPN frame-work. We develop new ways of dealing with this strange behaviour 
below, and apply the results of gravitational experiments to these theories. 

10.2 Scalar- Tensor Theories 

The post-Newtonian limit of the scalar-tensor theories of gravity has been well stud- 
ied numerous times in the past. We present here a brief illustration of how the 
post-Newtonian coefficients are obtained, and how they can be used to constrain 
the theory. 

In order to fit into the PPN formalism the metric and scalar field must be 
subjected to a perturbative expansion of the form 
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where 



/loo -0(2) + 0(4) 

V ~ 0(3) 

V ~ 0{2) 

ip ~ 0(2) + 0(4). 

These expansions can then be substituted into the scalar-tensor field equations 
()2.1.8|) and ()2.1.9|) . which can be solved order by order in perturbations to obtain 
(after an appropriate gauge choice) |191j 

2 , . /3 + 2cj 



goo = -l + 2U- 2(1 + A)U' + 4 $i (10-2.1) 

/l + 2cj \ 

g,u = f 1 + 2 (^) u] S,^ 



2 + 

l/10 + 7cu\ ^ 1^^ 



where 



^ _ duo/dcj) 



(4 + 2cj)(3 + 2cu)^ 
and Newton's constant has been set to 



G 



4: + 2uj\ 1 
'o' 



3 + 2cu 



A direct comparison of (jlO.l.lll and ()1(). 2.111 allows the PPN coefficients to be read 
off as 

l+UJ 

^ ~ 2 + uj 
(3=1 + K 

Cl = C2 = Cs = C4 = ttl = "2 = 0. 

The Brans-Dicke theory can be particularly well constrained by observations con- 
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straining the PPN coefficients. For the Brans-Dicke theory uj =constant, so A = 0, 
and the constraints on the coefficient 7, 

7 = 1 + (2.1 ±2.3) X 10-^ 

from observations of the Shapiro time delay of radio signals from the Cassini space- 
probe 113, gives the constraint 

u > 40000. 

This constraint is very limiting indeed for the Brans-Dicke theory, but less so for 
theories in which u is allowed to be a function of 0. For these theories, whilst 
the present day value of u in the solar system is tightly constrained, the value of 
u at other points in the Universe's history (particularly in the early Universe) is 
less restricted. To better constrain these theories we need to obtain limits on the 
parameter jS or, even better, use cosmological constraints from the early Universe, 
such as those which were discussed in the last chapter. 

10.3 Fourth-Order Theories 

In order to calculate the classical tests of metric theories of gravity (i.e. bending 
and time-delay of light rays and the perihelion precession of Mercury) we require 
a spherically symmetric solution to the C = R^~^^ field equations ()2.1.15p . Due to 
the complicated form of these equations we are unable to find the general solution; 
instead we propose to use ffist-order solutions around attractors as r — 00 or t — 00. 
This method should be applicable to gravitational experiments performed in the 
solar system as the gravitational field in this region can be considered weak and we 
will be considering experiments performed at large r (in terms of the Schwarzschild 
radius of the massive objects in the system) and at late times. 

10.3.1 Static Space-Time 

In this section we will use the static and spherically symmetric solution found pre- 
viously ()fi.2.2j] . We choose to arbitrarily set the oscillatory parts of the solution 
to zero, and hence ensure that the gravitational force is always attractive. This 
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considerable simplification of the solution also allows a straightforward calculation 
of both null and time-hke geodesies which can be used to compute the outcomes of 
the classical tests in this space-time. 



Solution in isotropic coordinates 

Having removed the oscillatory parts of the solution we are left with the part corre- 
sponding to the exact solution ()(j.2.2|l . Making the coordinate transformation 



\ / (l-2i+442) 

^(l-25+4<52)/{l-rf) = ( 1 _ ^ 1 (i-2*-2*^) 

4fV 



the solution ()f).2.2j) can be transformed into the isotropic coordinate system 

ds^ = -A{f)dt'^ + B{f){df^ + f'^ide'^ + sin^ Odcj)'^)) (10.3.1) 

where 

/ \ 2 / \ 2(1+4S) 

I \ I \ (1-25+442) 

24(1+24) / C \ I C 

A(r) = r V(l-24-242)(l-24+442) ]^ _| I 



(1-24+442) II ; (1-24+442) 

(1-25-242) y y ^-y (1-24-242) 
and 

4(1-4) 



_2+2 / C 

B{f) = f V{l-24-242)(l-24+442) M 



(1-24+442) 



(1-24+442) 
(1-24-242) 



This is, to linear order in C, 

(l-5)(l-25) C 



24(1+24) 



A{f') = f V(l-25-242)(l-24+442) ]^ ^ 



;i - 2(5 + 4(52) / (1-24+442) 



r 



(1-24-242) 



and 



_B(f) = f V{l-25-242){l-24+442) ]^ ^ ^ 



(1 - 2(5 + 4(52) / (1-24+442) 
(1-24-242) 
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Newtonian limit 



We first investigate the Newtonian limit of the geodesic equation in order to set the 
constant C in the solution ()1U.3.1|) above. As usual, we have 



~ 00 



where $ is the Newtonian gravitational potential. Substituting in the isotropic 
metric this gives 

v$ = riO.3.2) 



2B{f) 



6{l + 26)i 



l-2c5-2«2 



v/(l-25-252)(i_25 + 452) 

3(l-2g) 

V(l-24-242)(i_2i+452) 

+ 0{C 



(1 -85 + 4^2)^^ V{l-2*-2*^){l-2.+4*2) 



2v/(l - 2(5 - 252)(i - 25 + 4^2)3 

The second term in the expression goes as ~ f~2+'^('5^) and so corresponds to 
the Newtonian part of the gravitational force. The first term, however, goes as 
~ f-i+oc^^) g^j-^^ Yiss no Newtonian counterpart. In order for the Newtonian part 
to dominate over the non-Newtonian part we must impose upon 5 the requirement 
that it is at most 

\ r 

If S were larger than this then the non-Newtonian part of the potential would dom- 
inate over the Newtonian part, which is clearly unacceptable at scales over which 
the Newtonian potential has been measured and shown to be accurate. 

This requirement upon the order of magnitude of 6 allows p().3.2|l to be written 

v$ = ^Jjic^ - i^^Som + 

where expansions in C have been carried out separately in the coefficients and the 
powers of r of the two terms. 



177 



Comparison of with the Newtonian force law 

allows the value of C to be read off as 



C= -2Gm + Oi5). 



Post-Newtonian limit 

We now wish to calculate, to post-Newtonian order, the equations of motion for test 
particles in the metric ()10.3.1|) . The geodesic equation can be written in its usual 
form 

(fx^ dx^ dx^ 
L p/^. = 

dX^ dX dX ' 

where A can be taken as proper time for a time-like geodesic or as an affine parameter 
for a null geodesic. In terms of coordinate time this can be written 



d'^x^' / _ ^\ dx^dx^ ^ 
df" [ dt dt dt 



We also have the integral 



dx'^ dx^ 
Qij— — — = S 



(10.3.4) 



(10.3.5) 



dt dt 

where S = ~1 for particles and for photons. 

Substituting p0.3.H) into ()10.3.4|) and p0.3.5p gives, to the relevant order, the 
equations of motion 



d^ 

'dF 



Gm 



1 + 



dt 



G'^m? Gm rfx rfx 

4 — + 4^^e 



dt dt 



d'x 



dt 



.Gm 



er - 4— er + 5^er + 0(GW) (10.3.6) 



and 



(ix 



dt 



1-4- 



Gm S 



r28 



2^;^ + 0(GW) 



(10.3.7) 



(In the interests of concision we have excluded the 0(5^) terms from the powers of 
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r; the reader should regard them as being there imphcitly). The first three terms 
in equation ()l().H.fi|l are identical to their general relativistic counterparts. The 
next two terms are completely new and have no counterparts in general relativity. 
The last term in equation ()10.3.6|) can be removed by rescaling the mass term by 
m —>■ m{l + 6); this has no effect on the Newtonian limit of the geodesic equation 
as any term Gm6 is of post-Newtonian order. 

The bending of light and time delay of radio signals 

From equation ()10.3.7|) it can be seen that the solution for null geodesies, to zeroth 
order, is a straight line that can be parametrised by 

X = n(t - to) 

where n ■ n = 1. Considering a small departure from the zeroth order solution we 
can write 

X = n(t - to) + xi 

where xi is small. To first order, the equations of motion ()10.3.6|) and (I10.3.7|l then 
become 

c/^x Gm Gm , ^ , . 

— = -2-^er + 4-^(n.er)n (10.3.8) 

and 

(ix Gm 
„.- = -2— (10.3.9) 

Equations ()l().3.8j) and ()1().3.9|) can be seen to be identical to the first-order 
equations of motions for photons in GR. We therefore conclude that any observations 
involving the motion of photons in a stationary and spherically symmetric weak field 
situation cannot tell any difference between GR and this R^~^^ theory, to first post- 
Newtonian order. This includes the classical light bending and time delay tests, 
which should measure the post-Newtonian parameter 7 to be one in this theory, as 
in GR. 

Perihelion precession 

In calculating the perihelion precession of a test particle in the geometry ()1(). 3.111 
it is convenient to use the standard procedures for computing the perturbations of 
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orbital elements (see |179j and |162j ). In the notation of Robertson and Noonan 
[162j the measured rate of change of the perihelion in geocentric coordinates is given 
by 

— = ; — cos^H ; sm0 (10.3.10) 

at he he 

where p is the semi-latus rectum of the orbit, h is the angular- momentum per unit 
mass, e is the eccentricity and TZ and J are the components of the acceleration in 
radial and normal to radial directions in the orbital plane, respectively. The radial 
coordinate, r, is defined by 

r = - ^ (10.3.11) 

(l + ecos0) ^ ' 

and is the angle measured from the perihelion. We have, as usual, the additional 
relations 

p = a(l — e^) 

and 



h= ^/G^ = r^^. (10.3.12) 
From p0.3.(ij) . the components of the acceleration can be read off as 

7^ = + + 4 — + -v^-A— 10.3.13 

ry A ry A ryij ry ry ry 

and 

CtTTI 

j = ^—^vnvj (10.3.14) 
where we now have the radial and normal-to-radial components of the velocity as 

eh 

vn = — sm 
P 

Vj = —{1 + ecos0) 
p 

and v'^ = v\ + v'j. In writing p().3.13j) . the last term of ()1().3.6|) has been absorbed 
by a rescaling of m, as mentioned above. 

The expressions p0.3.13|) and p0.3.14|) can now be substituted into ()10.3.10|) and 
integrated from = to 27r, using ()10.3.1H) and ()10.3.12|) to write r and dr in terms 
of and dcf). The perihelion precession per orbit is then given, to post-Newtonian 
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accuracy, by the expression 



Aa; 



QtcGui 



2ti5 



{ 



(l + 45)a(l-e2) 



(10.3.15) 



Gm 



The first term in ()10.3.15|) is clearly the standard general relativistic expression. 
The second term is new and contributes to leading order the term 



27ra A -e^X 



Gm V e2 ) 



Comparing the prediction p().3.15j) with observation is a non-trivial matter. The 
above prediction is the highly idealised precession expected for a time-like geodesic 
in the geometry described by ()10.3.H) . If we assume that the geometry p0.3.1|) is 
a good approximation to the weak field for a static Schwarzschild-like mass then 
it is not trivial to assume that the time-like geodesies used to calculate the rate 
of perihelion precession ()10.3.15|) are the paths that material objects will follow. 
Whilst we are assured from the generalised Bianchi identities |124j of the covariant 
conservation of energy-momentum, T"*.fe = 0, and hence of the geodesic motion of an 
ideal fluid of pressureless dust, f/*[/-'j = 0, this does not ensure the geodesic motion 
of extended bodies. This deviation from geodesic motion is known as the Nordvedt 
effect |141j and, whilst being zero for GR, is generally non-zero for extended theories 
of gravity. From the analysis so far it is also not clear how orbiting matter and other 
nearby sources (other than the central mass) will contribute to the geometry ()1().3.H] . 

In order to make a prediction for a physical system such as the solar system, and 
in the interests of brevity, some assumptions must be made. It is firstly assumed that 
the geometry of space-time in the solar system can be considered, to first approxi- 
mation, as static and spherically symmetric. It is then assumed that this geometry 
is determined by the Sun, which can be treated as a point-like Schwarzschild mass 
at the origin, and is isolated from the effects of matter outside the solar system and 
from the background cosmology. It is also assumed that the Nordvedt effect is neg- 
ligible and that extended massive bodies, such as planets, follow the same time-like 
geodesies of the background geometry as neutral test particles. 
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In comparing with observation it is useful to recast in the form 

GTcGm 

where 

This allows for easy comparison with results which have been used to constrain the 
standard post-Newtonian parameters, for which 

A = ^(2 + 27-/3). 

The observational determination of the perihelion precession of Mercury is not clear 
cut and is subject to a number of uncertainties; most notably the quadrupole mo- 
ment of the Sun (see e.g. We choose to use the result of Shapiro et. al. |174j . 

A = 1.003 ±0.005, (10.3.16) 
which for standard values of a, e and m [3^ gives us the constraint 

5 = 2.7±4.5 X 10"^^ (10.3.17) 

In deriving ()10.3.16p the quadrupole moment of the Sun was assumed to correspond 
to uniform rotation. For more modern estimates of the anomalous perihelion advance 
of Mercury see |155j . 

10.3.2 Non-Static Space-Time 

We will now proceed to calculate the equations of motion of particles following 
geodesies of the space-time fl7.2.1|) . to post-Newtonian order. It will be shown that 
not only are the terms due to the linear perturbations different to the static case 
(as should be expected as the perturbations themselves have been shown to be 
background dependant) but that the background itself contributes an extra term to 
the post-Newtonian equation of motion. 
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Equation of motion 

We will now proceed to calculate the geodesies of the non-static solution ()7.2.2j) . In 
doing this we will neglect the contribution of the mode to ()7.2.2|) so that we only 
take into account the modes which go as r~^, in the limit 5-^0. This is the mode 
corresponding to the linearisation of the exact solution ()7.2.H1 . 
The geodesic equation can be written, as usual, in the form 

d'^x^ dx^ dx^ 

L Y'^. . = 

dX^ dX dX 

where A can be taken as proper time for a time- like geodesic, or as an affine parameter 
for a null geodesic. In terms of coordinate time this can be re-written as 

di ( -^ti \ doc dx"^ „ ^ „ „ „\ 

+ n-r°,-jr -^-3r = o. (lo.s.is) 



dt^ \ dt J dt dt 

Substituting the linearised solution into this equation will then give the equations 
of motion for test particles in this space-time. 

Substituting ()7.2.H) into ()10.3.18p gives, to post-Newtonian order, the equation 
of motion 

571 = -—-r + 4(1 - - (1 - 2S)— ^- j 

+ 4(1-^)— e,.^^-H^ I (10.3.19) 



where 



Ida 

n = -—. 

a dr 



The conformal time coordinate r is defined by dt = adr and C4 has been set by the 
appropriate Newtonian limit. The equations of motion for this solution are consid- 
erably simpler than those of the static solution ()6.2.2j) . but still differ from those of 
GR in significant ways. All of the terms except the last in equation 1)10.3. 19|1 have 
GR counterparts and the powers of r in these terms are all the same as in the general 
relativistic case. The premultiplicative factors of these terms are, however, modified 
and can be described adequately within the frame-work of the PPN approach by 
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assigning 

(3 = 1 and 7 = 1 - 25. 

By making this identification, the constraints on 7 from observations of the Shapiro 
time delay of radio signals from the Cassini space probe |17] can be used to impose 
upon 5 the constraint 

(5 = -1.1 ± 1.2 X 10"^ (10.3.20) 

As well as the usual effects associated with 7 — 1 being non-zero there are extra 
effects in this space-time due to the last term in ()10.3.19p . This term is proportional 
to the velocity of the test-particle (when v « c) and is zero for photons. For this 
reason we identify it as a friction term, with the friction coefficient being given by Ti. 
This 'friction' is a purely gravitational effect and is not due to any non-gravitational 
interaction of test particles with any other matter. 

It is clear that the constraints which can be imposed on this theory depend 
crucially on whether or not the space-time of the solar system can be assumed to 
be static. This is a non-trivial matter, as Birkhoff's theorem is no longer valid. 
Taking the usual assumption of the geometry of the solar system being static leads 
us to a solution with a non-trivial asymptotic dependence on r, as r ^ 00. These 
unusual asymptotics can be removed, but this appears to require that the usual 
assumption of the solar system having a static geometry be abandoned. Further 
study is required on this subject to determine which of these solutions is the most 
appropriate for application to the solar system. 
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CHAPTER 11 



Conclusions 



The objective of this work has been to develop the framework of alternative theo- 
ries of gravity, to investigate the solutions to their field equations and to identify 
the principle ways in which they deviate from GR. This framework has then been 
used to study the ways in which GR can be considered to be special or unique, to 
identify new behaviour that was not previously present and to place observational 
constraints on deviations from the standard theory. In achieving these goals a num- 
ber of approaches have been used in both cosmological and weak-field environments. 

The focus of the analysis of cosmological solutions has been on, but not limited to, 
studies of FRW universes. These universes are appealing for a number of reasons. As 
well as having the same symmetries as our own observable Universe (on the largest 
scales) their high degree of symmetry is particularly useful for finding simple exact 
solutions to the complicated field equations of the theory These solutions have been 
used to model a variety of physical processes such as structure formation, primordial 
nucleosynthesis and microwave background formation. All of these processes have 
been seen to behave differently from the corresponding cases in the standard model. 
In this way we have been able to see exphcitly what effect such modifications to GR 
would have on the evolution of the Universe, as well as on the physical processes 
occurring within it. 
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Further to studying FRW universes we have also considered a variety of other 
cosmological solutions. In particular we have found the conditions for the existence 
of Godel, Einstein static and de Sitter universes in a wide variety of theories. For the 
case of the Godel universe the existence of closed time-like curves, and hence time 
travel into the past, has been analysed in detail. Furthermore, anisotropic solutions 
of Bianchi type I have been found and used to study the approach to anisotropic 
singularities in this extended framework. Using these solutions it has been found that 
an infinite sequence of chaotic mixmaster oscillations does not occur in all theories, 
as it does in GR. As well as studying homogeneous cosmologies, investigation has 
also been made of the effects of inhomogeneity on cosmological models within this 
framework. We have used FRW universes to model spherical collapse and have found 
exact inhomogeneous, spherically symmetric and time-dependent vacuum solutions. 
The very existence of these solutions is forbidden in GR by Birkhoff 's theorem, and 
so they are of particular interest in looking for new behaviour which is not available 
in the standard model. 

The tightest constraints currently available on gravitational theory are from ex- 
periments in the solar system and astrophysical tests in the proximity of (almost) 
spherical massive bodies. With this in mind the spherically symmetric vacuum solu- 
tions to the field equations are of particular interest. We have found exact solutions 
for this situation, and analysed their behaviour. These solutions have been shown 
to exhibit entirely new behaviour not previously obtainable in GR, including time 
dependence and non-asymptotic flatness. The weak field limit of these solutions has 
been found and the geodesies calculated to Post-Newtonian order. Comparison of 
these solutions with observational results has allowed tight constraints to be placed 
on the theory. 

In summary, the principle new results in this work are that it has: 

• Investigated the spatially flat FRW solutions of theories with Lagrangian den- 
sity C^BJ". 

• Found the conditions for the existence of Godel, Einstein static and de Sitter 
universes in a wide class of fourth-order theories. 

• Investigated the phenomenon of closed time-like curves in these theories. 
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• Found exact Bianchi type I anisotropic solutions for a wide class of scale in- 
variant fourth-order gravity theories. 

• Investigated the approach to anisotropic singularities in these theories. 

• Found exact static and spherically symmetric solutions in vacuum for i?" the- 
ories and investigated their stability. 

• Found exact inhomogeneous cosmological solutions for scalar-tensor and i?" 
theories. 

• Explicitly shown the lack of validity of Birkhoff 's theorem in scalar-tensor and 
i?" theories. 

• Investigated the spherical collapse model in scalar-tensor theories, and the 
spatial variations in G that can result from it. 

• Investigated the possibility of generalising the usual scalar-tensor theories by 
allowing a transfer of energy between scalar and matter fields. 

• Used primordial nucleosynthesis to impose constraints upon scalar-tensor and 
i?" theories, and their resulting cosmologies. 

• Considered microwave background formation in RP' theories. 

• Imposed constraints on i?" theories by modelling the solar system as either a 
static or non-static spherically symmetric vacuum space-time. 

This work has confirmed the excellent agreement between GR and gravitational 
experiment and observation. It has uncovered new effects of modifying the stan- 
dard theory and it has shown that such modifications are very strongly constrained, 
particularly in the low curvature regime of the solar system. 

There are a variety of directions in which future research on this subject can 
proceed. Firstly, new effects and behaviours that have been found in this work 
can be used to model other physical processes. Particular processes of interest 
may be particle production near anisotropic singularities, modelling of anomalous 
observations (such as the Pioneer anomaly) or effects on structure formation and 
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the dynamics of astrophysical bodies. A second avenue of investigation would be 
to further investigate the effects that can be found to arise in theories of gravity 
which do not allow Birkhoff 's theorem to be formulated. A start has been made on 
this by identifying exact vacuum solutions in Brans-Dicke and i?" theories which are 
spherically symmetric and either static or non-static. More work is required to fully 
understand how these solutions should be interpreted and what their physical effects 
would be. Other future research may focus on the effects of self-interacting scalar 
fields in some of the contexts that have been explored here. The gravitational scalars 
that have been considered in this work have all been massless. Self-interactions are 
of direct interest for problems such as the late-time acceleration of the Universe as 
well as early-time behaviour, such as that which occurs during inflation. Understand 
these theories, and their physical effects, would be of direct interest for considerations 
of these scenarios. 
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